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Abstract 

SERDP  project  MM-1537  entitled  “Handheld  Frequency  Domain  Vector  EMI  Sensing  for  UXO  Dis¬ 
crimination”  is  complete.  The  original  objectives  of  MM1537  consisted  of  both  hardware  and  software 
goals.  On  the  hardware  side,  the  goal  was  to  develop  an  innovative  vector  (multi-axis)  handheld  UWB 
electromagnetic  induction  sensor  in  the  frequency  domain,  with  precise  3-D  positioning,  for  close  inter¬ 
rogation  of  anomalies.  This  new  instrument  would  allow  a  new  higher  level  of  UXO  discrimination  in 
the  vicinity  of  a  noteworthy  magnetic  response.  Data  processing  software  accompanying  the  new  instru¬ 
ment  would  also  be  developed.  The  goal  here  was  to  develop  clutter-tolerant  signal  processing  for  UXO 
discrimination  using  the  data  provided  by  the  new  sensor,  based  on  new,  high-fidelity,  physically  com¬ 
plete  forward  modeling  (the  Standardized  Excitation  Approach),  rigorous  instrument  characterization, 
and  on  new  processing  techniques. 

We  have  constructed  the  GEM-3  along  with  a  positioning  system  with  subcentimeter  accuracy  at 
a  range  of  about  2  meters.  This  new  instrument  has  been  successfully  used  to  acquire  data  both  under  lab 
conditions  and  in  test  plots  with  emplaced  targets.  Results  from  all  blind  tests  (see  Sec.  5)  indicate  the 
GEM-3  can  acquire  diverse  and  accurate  vector  data  with  a  SNR  similar  to  or  better  than  prior  GEM 
models.  Single  target  inversion  results  for  laboratory  blind  tests  were  100%  accurate.  Blind  tests  from 
the  test  plots  at  CRREL,  with  either  one  or  two  targets  in  close  proximity,  were  also  very  encouraging 
though  less  accurate. 

A  novel  “beacon”  positioning,  which  uses  the  primary  field  of  the  GEM-3  itself  to  locate  the 
sensor  head,  was  developed  and  successfully  deployed.  Early  ground  and  stiffness  sensitivities  were 
overcome  with  limited  redesigning  of  the  main  sensor  head.  The  GEM-3.^+  successfully  acquires  three 
axis  data  of  the  secondary  field  near  the  center  of  the  instrument  head  while  remaining  portable  for 
operators  and  reproducable  by  Geophex  Ltd..  Details  of  the  GEM-3 hardware  development  can  be 
found  in  Sec.  3. 

Our  group  has  developed  algorithms  and  models  which  accompanied  the  development  of  the  GEM- 
3^^  hardware.  These  algorithms  include  EMI  models  of  ground  response  including  the  effects  of  a 
non-flat  surface/air  interface,  methods  to  calibrate  EMI  insttuments  based  on  the  magnitude  of  the  sec¬ 
ondary  field  from  canonical  objects,  models  of  the  broadband  EMI  response  of  ellipsoids,  extensions  of 
the  Standardized  Excitations  Approach  (SEA),  and  improvements  to  the  Normalized  Surface  Magnetic 
Source  (NSMS)  model.  Details  of  these  models  and  algorithms  can  be  found  in  Sec.  4. 

The  GEM-3  is  now  a  fully  functioning  EMI  instrument  capable  of  vector  sensing  of  magnetic 
anomalies  while  being  well  located  within  a  limited  range.  For  precise,  queued  interrogation  of  anoma¬ 
lies,  the  GEM-3  provides  diverse,  accurate,  frequency  domain  data  of  the  secondary  EMI  field  suit¬ 
able  for  inversion  and  discrimination  with  high  fidelity,  rigorous  models. 
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1  OBJECTIVES 


1  Objectives 

•  Develop  innovative  veetor  (multi-axis)  handheld  UWB  (ultra-wideband)  electromagnetic  in¬ 
duction  sensor  in  the  frequency  domain,  with  precise  3-D  positioning,  for  close  interrogation 
of  anomalies,  i.e.  for  a  new  higher  level  of  UXO  discrimination  in  the  vicinity  of  a  note¬ 
worthy  magnetic  responses  identified  either  by  the  new  instrument  or  most  likely  by  other 
surveying. 

•  Develop  clutter-tolerant  signal  processing  for  UXO  discrimination  using  the  data  provided 
by  the  new  sensor,  based  on  new,  high-fidelity,  physically  complete  forward  modeling  (the 
Standardized  Excitation  Approach),  rigorous  instrument  characterization,  and  on  new  pro¬ 
cessing  techniques  (magnetic  field  pattern  matching,  general  inversion  for  equivalent  mag¬ 
netic  charge  capacity,  optimization  by  Differential  Evolution,  and  classification  with  Support 
Vector  Machines). 

This  project  builds  on  SERDP  SEED  work  (project  #1353)  in  response  to  the  needs  expressed 
in  UXSON-05-03,  by  developing  innovative  sensors  and  signal  processing.  The  aim  was  and  is 
to  achieve  a  new  level  of  UXO  survey  target  characterization  and  false  alarm  rejection.  The  work 
will  allow  the  Army  and  other  DoD  components  to  carry  out  mandated  UXO  cleanup  more  safely, 
effectively,  and  efficiently.  The  new  fast  processing  and  the  portable  nature  of  the  products  might 
also  allow  application  in  tactical  situations  and  for  active  range  maintenance,  as  well  as  for  cleanup 
away  from  areas  of  conflict. 

The  positioning  and  measurement  systems  here  are  intended  to  link  to  dual  mode  navigation  as 
may  be  developed  under  other  SONs  (Statement  of  Need).  That  is,  these  devices  would  provide 
the  high  density,  fine  mode,  while  positioning  reference  points  are  linked  to  locations  identified  in 
the  coarse  mode.  More  generally,  the  survey  and  processing  systems  contemplated  are  designed 
to  be  most  effective  under  some  degree  of  cuing.  The  default  assumption  is  a  high  level  of  cuing, 
in  which  prior  information  indicates  the  approximate  locale  around  an  individual  anomaly,  on  the 
order  of  meters.  In  both  the  small  and  larger  scale  areal  coverage,  surveying  with  a  low  level  of 
cuing  is  also  possible.  That  is,  the  cuing  would  simply  consist  of  identification  of  some  specific  area 
some  lO’s  to  ;^100  m  across,  within  which  the  new  system  would  accomplish  both  coarse  detection 
and  fine  interrogation.  Reception  of  the  horizontal  magnetic  field  components  would  help  indicate 
the  direction  of  offset  from  a  target  location;  and  the  user  would  then  perform  arbitrarily  detailed 
sweeps  around  the  indicated  anomaly.  Whether  the  system  to  be  developed  is  employed  in  both 
coarse  and  fine  mode  of  surveying  or  only  in  the  latter,  the  primary  focus  of  the  work  is  on  UXO 
discrimination  during  close  interrogation  of  anomalies. 

This  work  is  designed  in  part  to  address  realistic  conditions:  terrain  with  lumps,  bumps,  moguls, 
trenches,  embankments...  and  possibly  sparse  vegetation.  High  density  UWB  EMI  data  will  be 
obtained  from  a  new  man-portable,  vector,  ED  (frequency  domain)  instrument  with  precise  posi¬ 
tioning  (sub-cm  local  precision).  The  new  system  is  designated  below  as  the  GEM-3Dplus  (GEM- 
3^+),  meaning  a  fully  integrated  GEM-3D  plus  positioning.  The  instrument  will  be  “handheld” 
in  the  sense  that  a  reasonably  hardy  surveyor  will  be  able  to  move  the  sensor  head  about  freely 
in  sweeps,  tilts,  and  elevations,  for  perhaps  a  couple  of  hours,  without  relief.  Ergonomic  consid¬ 
erations  are  not  given  top  priority  here,  as  the  result  will  be  primarily  a  research  instrument.  The 
new  data  features  provided  by  the  instrument  will  be  exploited  by  a  new  generation  of  original 
processing  approaches  based  on  the  SEA  and  NSMS,  as  well  as  on  other  approaches  described 
below.  The  potential  benefits  envisioned  in  some  of  the  proposed  processing  approaches  depend 
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on  the  new  level  of  data  eontent.  Also,  eonversely,  the  precision  and  diversity  of  data  contemplated 
only  make  sense  if  a  new  generation  of  processing  can  take  advantage  of  it.  Thus  we  will  develope 
the  instrumentation,  survey  techniques,  and  processing  together. 

Relative  to  the  system  produced  during  the  antecedent  SEED  work,  the  new  GEM-3  will 
incorporate  a  number  of  desirable  developments  as  listed  below: 

1.  New  small  (fii  5  x  8  cm),  stackable  DSP’s  (Digital  Signal  Processor),  each  1  cm  in  thickness, 
so  total  stacked  height  for  3  DSP’s  ~  approx  3  cm.  The  DPS’s  synchronization,  power, 
and  communication  will  all  be  on  a  bus  line  (no  wiring)  with  only  one  communication  line, 
and  a  single  RS485  connection  at  IMBit/sec  to  an  RS485/USB2  converter  as  the  link  to  a 
computer; 

2.  Each  DSP  will  provide  an  additional  RS232  port  for  GPS  (Global  Positioning  System)  or  any 
other  external  device  such  as  the  laser  positioning  system; 

3.  The  frequency  range  will  be  increased  from  48KHz  to  96kHz,  with  low  frequency  limit  of 
30  Hz; 

4.  Earger  diameter  sensor  head  (54  cm)  for  more  uniform  primary  field  and  greater  depth  of 
view,  without  significant  loss  of  portability; 

5.  The  3  DSP’s  will  be  configured  to  increase  the  system  performance  (convolution  is  done  in 
parallel  for  all  three  channels); 

6.  The  transmitter  will  be  a  separate  module  connected  to  the  DSP  via  RS422  digital  audio 
communication; 

7.  New  higher  power  transmitter  with  four  times  the  voltage  of  the  older  system. 

8.  Standardized  interface  and  software  for  the  positioning  system  connection,  allowing  EM 
(ElectroMagnetic)  data  synchronization  with  timing  pulses  from  GPS,  inertial,  or  laser  posi¬ 
tioning; 

9.  Integration  of  positioning  system  of  choice  with  the  EM  system  so  that  EM  and  position  data 
files  are  merged; 

10.  The  device  will  be  easily  reproducible. 

The  minimum  CPU  required  for  system  control  will  be  on  board,  while  wireless  links  to  other 
computers  will  serve  for  data  storage  and  processing.  This  will  allow  quasi-real  time  application 
of  the  processing,  including  possibly  parallelized  processing,  which  can  then  guide  the  testing  in 
a  flexible  way.  In  the  least,  it  will  streamline  design  by  separating  the  control  and  the  storage/pro¬ 
cessing  functions. 

2  Background 

Detection  and  especially  discrimination  of  buried  unexploded  ordnance  (UXO)  from  omni-present 
metallic  clutter  is  a  persistent,  expensive,  and  pressing  problem.  While  it  is  not  yet  clear  whether  or 
which  other  technologies  might  be  best  used  in  conjunction  with  electromagnetic  induction  (EMI) 
for  UXO  discrimination,  it  is  clear  that  EMI  is  currently  a  front  runner  in  the  development  of 
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new,  more  effeetive  approaehes.  Virtually  all  EMI  sensors  transmit  a  ’’primary”  field  and  reeeive 
signals  (’’seeondary”  field)  using  wire  eoil  antennas.  Frequeney  domain  (FD)  EMI  systems  use 
sueh  eoils  to  transmit  (Tx)  a  “primary”  magnetie  field  with  seleetion  of  frequeneies  over  a  ehosen 
band  -  possibly  only  a  single  frequeney  -  simultaneously  reeeiving  (Rx)  and  reeording  responses 
to  those  frequeneies  (the  ’’seeondary”  field).  Thus  a  eentral  problem  in  frequeney  domain  sensors 
is  the  isolation  of  the  reeeiving  eoils  from  the  primary  field,  the  latter  being  very  large  relative 
to  the  seeondary  field,  partieularly  when  the  reeeivers  are  near  the  transmitters,  as  eontemplated 
here.  Physieal  separation  of  the  Tx  and  Rx  eoils  is  a  possible  approaeh.  However  this  typieally 
degrades  spatial  resolution  and  handy  manipulation  of  the  instrument.  Therefore,  remaining  with 
eo-loeated  Tx  and  Rx  eoils,  we  implement  a  “bueking”  eoil,  whieh  transmits  a  field  that  opposes 
the  primary  field  at  the  loeation  of  the  reeeiving  eoil.  By  eaneeling  out  the  primary  field  at  the 
reeeiver  as  mueh  as  possible,  one  is  left  (ideally)  with  only  the  seeondary  field  from  the  target. 
Other  approaehes  are  possible  for  finding  the  seeondary  field  amidst  the  mueh  stronger  primary 
field,  e.g.  using  differeneing  sehemes  between  staeked  or  symmetrieally  distributed  eoils,  but  the 
fundamental  problem  is  the  same. 

The  “bueking”  diseussion  above  is  partieularly  pertinent  when  one  attempts  3-D  veetor  defi¬ 
nition  of  the  seeondary  field.  The  geometrieal  arrangement  of  the  main  and  the  bueking  eoils  is 
erueial  to  sueeessful  suppression  of  the  primary  field  at  the  reeeiver.  However,  to  our  knowledge 
no  FD  EMI  instruments  have  been  produeed  to  date  that  are  eapable  of  suppressing  the  primary 
field  in  one  reeeiver  orientation  while  also  suppressing  it  in  other  direetions  as  well.  Most  imag¬ 
inable  arrangements  of  two  bueking  eoils  that  work  well  for  one  reeeiver  orientation  would  in  faet 
exaeerbate  the  problem  for  other  reeeiver  orientations.  This  observation  would  appear  to  favor 
the  development  of  time  domain  (TD)  systems  for  obtaining  veetor  EMI  response.  Most  TD  sys¬ 
tems  operate  by  transmitting  a  steady  signal,  saturating  a  metal  objeet  of  interest,  then  shutting  the 
transmitter  off.  A  reeeiver  then  reeords  the  seeondary  field  from  the  objeet  as  that  field  deeays  in 
response  to  the  sudden  shut  off.  The  reeeiver  only  operates  while  the  transmitter  is  off,  thereby 
dodging  the  primary  field.  While  this  bypasses  the  bueking  problem,  a  number  of  things  motivate 
us  to  proeeed  in  the  FD:  FD  systems  elaim  the  advantage  in  having  superior  eontrol  of  seleetion 
and  power  in  the  frequeney  eontent  of  data  produeed,  and  thus  in  the  equivalent  time  domain  sig¬ 
nal,  whieh  ean  readily  be  obtained  from  the  FD  data.  In  praetiee,  FD  systems  are  less  band  limited 
than  TD  systems  and  ean  therefore  offer  a  greater  equivalent  time  range  of  response  than  the  aetual 
TD  systems.  Consider  Fig.  2.0.1  below. 

While  there  has  been  progress  in  ’’very  early  time”  EMI  deviees,  in  general  TD  deviees  are 
unable  to  shut  off  the  transmitter  and  damp  the  effeets  of  the  transition  quiekly  enough  to  get  very 
early  time  data.  Information  in  sueh  data  is  the  equivalent  of  high  frequeney  FD  data.  The  latter 
is  of  strong  interest  beeause  it  provides  information  on  asymptotie  limits  of  seattering  behavior: 
The  indueed  eurrents  penetrate  the  target  negligibly  and  the  seeondary  fields  they  produee  do 
not  depend  on  the  type  of  metal  eneountered.  This  may  offer  hope  of  avoiding  very  substantial 
eomplieations  in  UXO  diserimination,  given  that  many  if  not  most  UXO  are  eomposites  of  different 
metals.  At  the  other  end  of  the  speetrum,  TD  signals  naturally  fade  as  the  target’s  responses  wind 
down.  This  means  that  very  late  time  data  is  unlikely  to  be  deteetable  above  the  noise.  Thus  the 
equivalent  of  low  frequeney  asymptotie  limits  will  be  poorly  defined  or  absent.  While  FD  systems 
must  also  struggle  with  some  ehallenges  in  the  equivalent  very  low  frequeney  range,  on  the  whole 
it  is  mueh  more  feasible  to  get  data  from  them  in  that  part  of  the  EMI  band.  Many  UXO  produee 
signifieant  response  patterns  in  the  lO’s  of  Hz.  Fow  frequeney  asymptotie  limits  and  fundamental 
signal  features  eontaining  basie  target  shape  information  appear  in  many  eases  only  well  below  30 
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log(t) 


Figure  2.0.1:  General  representation  of  received  time  domain  signal  (R)  vs  log  of  time. 

Hz  [SERDP  Project  1 122,  Final  Report].  To  get  at  the  TD  equivalent  of  these  frequencies  down  to 
1  Hz  using  TD  instruments  would  require  recording  out  to  the  order  of  tenths  of  a  second.  Highly 
unlikely.  At  least  at  the  research  level,  FD  FMI  measurements  have  been  made  successfully  over 
a  continuous  UWB  ranging  from  a  few  Hz  up  to  about  300  kHz,  roughly  the  TD  equivalent  of 
spanning  times  from  a  micro-seconds  to  tenths  of  a  second.  Recent  work  by  Geophex  on  TD 
systems  in  a  GFM  configuration  may  be  an  exception  to  the  limitations  cited  above,  and  we  look 
forward  to  encountering  that  new  device. 

Much  of  the  motivation  for  this  project  derives  from  a  desire  to  develop  further  a  new  generation 
of  FMI  forward  modeling;  to  apply  it  in  innovative  processing;  and  to  develop  instrumentation  that 
will  provide  the  data  most  likely  to  support  this  new  processing  to  best  advantage.  The  new  forward 
modeling  will  generally  be  based  on  the  Standardized  Excitations  Approach  (SEA),  though  all 
processing  will  not  be  restricted  to  SEA  implementations.  The  theory  and  detailed  development 
of  SEA  have  been  presented  in  various  fora  [1-12].  Here  we  only  present  a  general  sketch  of  the 
underlying  ideas  so  that  the  direction  taken  in  some  of  the  proposed  work  may  be  understood. 
The  essence  of  the  method  lies  in  its  ability  to  provide  physically  complete  forward  modeling  - 
including  all  target  heterogeneity,  near  and  far  field,  and  internal  interaction  effects  -  while  being 
fast  enough  for  inversion  and  classification  computations. 

During  antecedent  SERDP  SEED  work,  the  configuration  of  the  existing  Geophex  GEM-3  [4] 
was  used  as  a  starting  point  for  the  new  vector  device.  The  essential  difference  was  that  the  new  in¬ 
strument  would  receive  and  record  separately  three  orthogonal  magnetic  field  components  instead 
of  one.  This  was  accomplished  by  adding  two  receiver  coils  that  are  the  same  size  as  the  original 
single  coil  but  which  are  perpendicular  to  it  and  to  each  other.  These  additional  coils  escape  be¬ 
ing  swamped  by  the  primary  field  in  part  because  they  are  near  the  null  region  engineered  for  the 
original  in-plane  receiver  loop.  For  the  most  part,  however,  they  are  protected  by  their  orientation. 
The  plane  of  each  of  the  new  receiver  coils  is  completely  parallel  to  the  transmitted  field  lines  from 
the  perpendicular  transmitter  coils  (Fig.  2.0.2).  Because  the  receiver  coils  respond  only  to  rate  of 
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change  of  magnetic  flux  that  passes  through  them,  in  principle  the  additional  coils  do  not  register 
the  primary  field.  The  bench  level  instrument  that  was  constructed  and  tested  has  an  expanded 
bandwidth  relative  to  established  versions  of  the  GEM-3,  ranging  from  the  same  lower  limit  of  30 
Hz  up  to  at  least  98  kHz,  with  a  maximum  of  about  15  frequencies  recordable  in  each  channel,  i.e. 
for  each  orthogonal  signal  component  (using  slightly  fewer  is  recommended). 


Figure  2.0.2:  An  EMI  sensor  field  impinges  on  a  UXO,  which  tends  to  respond  most  strongly  along  its  axis, 
producing  a  response  that  strikes  the  receiver  at  an  arbitrary  angle  relative  to  the  vertical. 

The  bench  version  of  the  instrument  that  was  produced  appears  to  work  well  by  both  quan¬ 
titative  and  qualitative  criteria.  Proper  (and  illuminating)  symmetries  and  asymmetries  appear  in 
the  different  received  field  components  as  the  sensor  head  is  swept  over  UXO’s  and  other  objects. 
Measured  data  for  machined  spheroids  of  different  materials  and  shapes  compare  very  well  to  cor¬ 
responding  analytical  solutions  that  account  for  the  details  of  the  OEM’s  primary  field.  The  laser 
positioning  system  that  was  implemented  is  said  (by  the  supplier)  to  be  capable  of  sub-millimeter 
accuracy,  and  was  configured  to  provide  all  tilt  angles  in  a  version  of  the  instrument  dubbed  the 
GEM-3DE.  In  various  benign  tests  the  integration  of  the  system  performed  well. 

Regarding  the  EMI  measurement  system,  the  task  that  remains  is  proper  design  and  imple¬ 
mentation  so  that  it  is  melded  into  an  efficient  and  coherent  unit,  including  direct  integration  with 
the  positioning  system  via  a  conventionally  recognizable  timing  pulses.  Basically,  for  the  sake  of 
showing  that  the  basic  concept  of  ED  vector  EM  measurement  could  work,  the  SEED  work  ef¬ 
fectively  patched  together  in  parallel  three  traditional  GEM-3’s,  one  associated  with  each  receiver 
loop.  Given  what  was  learned  during  the  SEED  exercise  and  also  the  recent  advancement  in  com¬ 
ponents  at  affordable  prices  (e.g.  digital  signal  processors  -  DSP’s),  we  are  now  positioned  to 
construct  and  integrate  the  instrument  as  it  should  be  done.  It  has  also  become  clear  that  the  easiest 
and  most  reliable  way  to  perform  the  EM-positioning  integration  is  to  key  the  EM  measurement 
timing  to  a  pulse  from  the  positioning  system,  as  opposed  to  vice  versa. 

The  GEM-3D  produced  in  the  SEED  effort  used  a  laser  positioning  system:  the  ArcSecond 
laser  positioning  system.  The  laser  positioning  system  provided  accuracies  on  the  order  of  a  few 
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mm  of  the  instrument  head.  The  positioning  data  provided  by  this  system  allows  free  movement  of 
the  sensor  head  at  any  elevation  and  tilt  angle  thus  leading  to  a  rieh  and  diverse  data  set.  Data  from 
multiple  elevations  has  helped  in  distinguishing  target  signals  from  widespread  elutter  responses; 
and  in  the  general  inversion/optimization  ealeulations  we  have  done,  it  signifieantly  redueed  or 
eliminated  ill-eonditioning.  Also,  horizontally  distributed  observation  points  allow  “multi-look” 
data,  i.e.  views  that  stimulate  and  also  register  response  from  all  sides  of  an  objeet  to  be  studied. 
All  this  requires  very  flexible  positioning,  i.e.  both  antenna  movement  and  traeking  thereof.  In  faet, 
without  sueh  positioning  the  new  (and  older  GEM-3)  sensor  will  not  be  useful  for  the  next  gen¬ 
eration  of  proeessing  towards  whieh  we  strive.  However  signifieant  design  and  integration  issues 
remain.  Both  inadequate  ealibration  of  the  laser  system  and  espeeially  timing  disparities  between 
both  the  different  laser  reeeivers  and  between  the  laser  and  the  EM  system  preeluded  eontinuous 
motion  of  the  sensor  in  reeent  field  tests  at  the  Aberdeen  Proving  Ground  UXO  Standardized  Test 
Site.  This  makes  for  impraetieally  slow  surveying,  even  for  testing  purposes. 

Initially,  we  planned  to  use  the  AreSeeond  system  with  the  new  GEM-3^^  in  order  to  aehieve 
the  same  level  of  positioning  aeeuraey.  However,  early  on  in  the  projeet  we  reeonsidered  this  design 
deeision  in  light  of  reeent  developments  in  the  area  of  “beaeon”  positioning  systems  (see  See.  3.3) 
and  beeause  of  faetors  mentioned  above  regarding  the  GEM-3^”*'.  Beaeon  positioning  systems 
rely  on  the  primary  field  of  the  instrument  itself  to  provide  an  aeeurate  position  of  the  sensor  head. 
Using  external  reeeivers  to  deteet  and  sample  the  primary  field  of  the  GEM-3^’'',  the  loeation 
of  the  primary  field  (assumed  dipolar)  souree  ean  be  inferred.  The  design,  implementation,  and 
resulting  aeeuraey  of  this  system  is  detailed  below.  This  positioning  system,  in  theory,  is  simpler 
to  operate  and  maintain  than  the  AreSeeond  system.  Eor  example,  the  ealibration  of  the  system  is 
far  simpler  for  the  beaeon  system,  and  the  setup  time  is  also  greatly  redueed.  On  the  other  hand, 
the  range  of  this  system  is  less  than  the  laser  positioning  system. 

In  simulations  and  lab  measurements,  we  have  seen  distinet  benefits  from  veetor  data  in  the 
task  of  defining  unambiguously  the  response  of  UXO’s  and  of  other  objeets  to  the  basie  modes  of 
exeitation.  Veetor  or  “3-D”  data  means  that  three  orthogonal  eomponents  of  the  seattered  magnetie 
field  are  measured  at  the  reeeiver  loeation,  thereby  eompletely  defining  the  field  there  that  eonsti- 
tutes  the  objeet’s  response.  The  idealized  representation  in  Eig.  2.0.2  shows  how  a  UXO  response 
field  will  generally  strike  the  reeeiver  at  some  arbitrary  angle  relative  to  the  vertieal,  i.e.  relative  to 
the  direetion  in  whieh  the  established  sensors  reeord  the  signal. 

The  rest  of  the  report  proeeeds  as  follows.  Seetions  3  and  4  detail  the  materials  and  meth¬ 
ods  assoeiated  with  the  design  and  implementation  of  the  hardware  of  the  GEM-3  as  well  as 
the  software  analysis  algorithms  developed  in  the  eourse  of  this  projeet.  Seetion  4  ineludes  dis- 
eussions  on  some  eharaeteristies  of  the  GEM-3^'''  (See.  4.1),  the  soil  response  of  the  instrument 
(See.  4.2),  absolute  ealibration  of  the  instrument  (See.  4.4),  the  Generalized  Standardized  Exeita- 
tions  Approaeh  (GSEA,  See.  4.5),  the  NSMS  and  HAP  methods  (See.  4.6),  and  the  response  of 
the  GEM-3^~^  to  ellipsoidal  targets  (See.  4.7).  These  seetions  are  followed  by  blind  test  results 
both  from  lab  data  and  from  test  plot  data  at  CRREE  (See.  5).  The  Appendiees  eontain  supporting 
material  sueh  as  operating  proeedures  and  a  publieation  list  follow  the  eoneluding  seetion  (See.  6). 
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3  Materials  and  Methods:  Hardware  Development  of  GEM-3 Instru¬ 
ment 

3.1  GEM-3^+  Design 

3.1. a  GEM-3^+  Specifications 

The  primary  purpose  of  this  projeet  is  to  advanee  the  existing  GEM-3  EMI  sensor  to  ineorporate 
multiple  axis  viewing  (’’veetor”)  eapability  that  is  expeeted  to  enhanee  the  UXO  deteetion  and 
diserimination.  Relative  to  the  system  produeed  during  previous  CRREE  projeet,  the  new  GEM- 
3^'^  will  ineorporate  a  number  of  advaneements  as  listed  below: 

•  Earger  diameter  sensor  head  (54  em)  for  more  uniform  primary  field  and  greater  depth  of 
view,  without  signifieant  loss  of  portability 

•  The  frequeney  range  will  be  increased  from  48  kHz  to  96  kHz,  with  low  frequency  limit  of 
30  Hz 

•  Use  of  new  small  stackable  DSPs,  each  about  1  cm  in  thickness.  Synchronization,  power, 
and  communication  will  all  be  on  a  bus  line.  An  RS485  connection  at  IMb/sec  to  an 
RS485/USB2  converter  will  serve  as  the  link  to  a  computer 

•  Each  DSP  with  an  additional  RS232  port  for  GPS  or  any  other  external  device  such  as  the 
laser  positioning  system 

•  DSPs  configured  to  increase  the  system  performance  (convolution  is  done  in  parallel  for  all 
four  channels) 

•  A  separate  transmitter  module  connected  to  the  DSP  via  RS422  digital  audio  communication 

•  New  higher  power  transmitter  with  increased  drive  the  voltage 

•  Standardized  interface  and  software  for  the  positioning  system  connection,  allowing  EM  data 
synchronization  with  timing  pulses  from  GPS,  inertial,  or  laser  positioning 

•  Integration  of  positioning  system  of  choice  with  the  EM  system  so  that  EM  and  position  data 
files  are  merged 

•  Easily  reproducible. 

3.1. b  Early  Progress 

This  section  describes  our  early  design  efforts  for  the  GEM-3^+  sensing  head,  transmitter  and 
receiver  coils,  and  structural  assemblage  of  the  coils  and  mechanical  support  structures.  The  GEM- 
3^+  has  two  additional  orthogonal  vertical  coils  in  addition  to  the  existing  GEM-3  geometry. 
Eigure  3.1.5  shows  the  overall  appearance  of  the  GEM-3^”''  sensing  head.  The  design  provides 
the  following  features: 

•  Sensing  head:  1/2”  thick  double-skin  fiberglass  foam  core  synthetic  board  with  an  overall 
diameter:  55  cm; 

•  Outer  transmitter:  54  cm  in  diameter; 
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•  Receiver  coils:  one  circular  horizontal  coil  in  the  magnetic  cavity  plus  two  rectangular  ver¬ 
tical  coils.  The  vertically  offset  rectangular  shapes  reduce  the  mandatory  sensor  ground 
clearance 

•  Maximum  visibility  of  the  target  area  by  embedding  many  see-thru  holes  in  the  disk 

•  Minimum  sensor  head  weight; 

•  The  electronic  console  and  battery  to  be  mounted  at  the  end  of  the  handle  bar,  partly  to  act 
as  a  counter  weight. 

All  three  receiver  coils,  one  circular  and  two  rectangular,  are  designed  to  have  an  identical 
area-turn  product.  The  two  rectangular  coils,  as  shown  in  Fig.  3.1.5,  are  symmetrically  centered 
across  the  sensing  disk.  The  main  advantage  for  centering  is  that  all  three  receiver  coils  are  mathe¬ 
matically  co-located  (this  was  later  changed,  see  below),  which  simplifies  locating  and  interpreting 
targets.  Otherwise,  these  coils  may  move  arbitrarily  upward. 

Figure  3.1.1  shows  the  design  sheet  for  this  GEM-3^+.  Basic  design  parameters  include  the 
diameters,  turns,  and  wire  sizes  for  Txl  (outer  transmitter),  Tx2  (inner  bucking  transmitter),  and  Rx 
(the  horizontal  receiver  coil).  The  two  rectangular  vertical  coils  (RxV  1  &  RxV2  -not  shown  in  this 
design  sheet)  are  made  such  that  their  area-turn  product  is  identical  to  that  of  the  horizontal  coil. 
This  sheet  is  based  on  the  12VDC  supply,  which  is  a  mere  scale  factor  for  the  transmitter  moment. 
Higher  supply  voltages  are  used  in  the  final  design.  The  design  shows  a  healthy  transmitter  moment 
of  about  35  A-m2  at  750  Hz,  which  compares  favorably  with  the  previous  40-cm  diameter  version 
having  about  10  A-m2  also  at  750  Hz.  The  dipole  loss  due  to  Tx2  is  only  5.37%.  The  sheet  also 
shows  the  T  and  Q-responses  of  the  Rx  in  ppm  to  a  known  Q-coil  between  30  Hz-96  kHz.  The 
specs  and  location  of  the  Q-coil  are  also  indicated  in  the  sheet. 

Figures  3. 1.2-3. 1.6  show  some  preliminary  drawings  of  the  GEM-3^’''.  Eater  designs  posi¬ 
tioned  the  transverse  receiver  coils  to  be  more  flush  with  the  bottom  of  the  transmitter  coil  (see 
Eig.  3.2.6). 

3.2  GEM-3  Continuing  Development 

3.2. a  Cord  Sensitivity  Issue 

During  the  initial  stages  of  testing  the  GEM-3^+  at  CRREE,  we  noticed  an  extreme  sensitivity  in 
the  real  part  of  the  transverse  components  associated  with  any  movement  of  the  instrument.  We 
eventually  traced  this  sensitivity  problem  more  specifically  to  issues  of  cord  movement.  Whenever 
the  connecting  cord  from  the  sensor  head  to  the  electronics  box  was  moved  in  the  slightest  amount, 
the  transverse  components  of  the  received  magnetic  field  were  corrupted  by  a  tens  if  not  hundreds 
of  ppm.  Together  with  the  Geophex  engineers,  the  problem  was  traced  back  and  solved  so  that  this 
sensitivity  issue  was  resolved. 

3.2. b  Head  Deformation  Issue 

Prototype  grade  instruments  often  have  several  unforeseen  issues  arise  during  their  development. 
After  the  sensitivity  associated  with  cord  movement  was  resolved,  we  noticed  another  type  of  sen¬ 
sitivity  issue  whenever  we  physically  touched  or  moved  the  instrument.  As  it  turned  out,  there  was 
a  further  sensitivity  associated  with  deformation  of  the  relative  position  of  the  transverse  receiver 
coils  to  each  other.  We  decided  to  investigate  whether  it  was  the  deformation  of  the  entire  head  or 
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TXl 


Figure  3.1.1:  GEM-3^+  basic  geometrical  design.  Table  showing  coil  parameters  and  frequency  character¬ 
istics  omitted:  proprietary. 
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30  deg 


Figure  3.1.2:  GEM-3^+  preliminary  drawing  1. 
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50  deg 


Figure  3.1.3:  GEM-3 preliminary  drawing  2. 
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GEM-3D+  RX  coil  form 

Material:  %"  thick  white  SINTRA 
OD:  1.5  cm 
ID:  12.8cm 


Figure  3.1.4:  GEM-3^+  preliminary  drawing  3. 
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TX1 


30  deg 


Figure  3. 1.5:  GEM-3^^  preliminary  drawing  4.  Final  design  for  receiver  coils  is  shown  in  Fig.  3.2.6. 
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3.3cm  dia  1cm  deep 


Figure  3.1.6:  GEM-3^+  preliminary  drawing  5. 
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Figure  3.2.1:  Pictures  of  the  GEM-3^'*’  in  development  at  an  intermediate  stage  before  finalization  of  the 
sensor  head. 
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just  of  the  receiver  coils  that  was  causing  this  broadband  shift.  So,  we  squeezed  the  two  transverse 
receiver  coils  toward  each  other,  without  touching  the  main  head,  and  sure  enough  we  got  only  in 
phase  broadband  shifts  up  or  down  depending  on  which  direction  we  squeezed  in. 

Any  slight  misalignments  from  perfectly  perpendicular  to  the  primary  field  should  be  consid¬ 
ered  background  and  subtracted  if  it  was  constant.  The  problem  is  that  these  transverse  receiver 
coils  are  not  as  rigidly  fixed  in  their  geometry  with  respect  to  the  primary  field  as  the  Z  receiver  coil 
is.  The  transverse  receiver  coils  are  only  embedded  in  a  1  cm  square  cross  sectional  scaffolding 
that  is  much  less  robust  in  terms  of  torsional  and  linear  deformation.  Similarly  when  I  squeeze  the 
receiver  coils  toward  each  other,  even  if  I  deflected  them  only  by  a  few  microns,  it  has  this  large 
effect  of  up  to  100  or  200  ppm  in  the  data.  The  fact  that  these  deformations  have  a  broadband,  con¬ 
stant  (with  respect  to  frequency)  inphase  shift  leads  me  to  suspect  it  is  the  primary  field  response  I 
am  seeing. 

To  lessen  this  effect  in  the  short-term  we  simply  created  some  handholds  on  the  supporting 
plexiglass  structure  we  now  have  on  the  GEM-3^”*'  so  that  we  can  pick  it  up  from  the  center  bot¬ 
tom  of  the  instrument  without  putting  any  torque  on  the  main  head  (see  Fig.  3.2.2).  We  ended  up 
stiffening  the  head  with  plexiglass  as  shown  in  Fig.  3.2.2.  This  helped  the  problem  considerably, 
but  it  neither  looked  good  nor  was  it  very  portable  due  to  the  increased  weight.  The  engineers  at 


Figure  3.2.2:  Plexiglass  gussets  glued  to  GEM-3^+  sensor  head. 

Geophex  redesigned  the  GEM-3^’^  sensor  head  using  lightweight  foam  core  epoxy  coated  stiff¬ 
ening  members.  A  current  picture  of  the  GEM-3^+  including  these  and  members  is  in  Fig.  3.2.3. 


3.2.b.(l)  Final  Head  Design  Figure  3.2.6  shows  the  final  design  of  the  GEM-3^'’'  sensing  head. 
The  only  difference  from  the  design  above  is  that  the  two  vertical  coils  (Hx  and  Hy)  are  raised 
by  3.5  cm,  mainly  to  reduce  the  ground  clearance  distance.  This  was  done  at  the  request  of  the 
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Figure  3.2.3:  Foam-core  gussets  epoxied  to  GEM-3^+  sensor  head. 


CRREL  project  manager,  Dr.  Ben  Barrowes.  Figures  3.2.4  and  3.2.5  show  the  sensor  head  during 
the  construction  and  initial  wiring  stages. 

3.2.b.(2)  Sensor  Electronics  Figure  3.1.1  is  the  final  design  sheet  for  this  GEM-3^'''.  The  main 
change  from  the  previous  report  includes  the  transmitter  wire  gauge  that  changed  from  14  strands 
of  AWG#30  to  a  Litz  wire  of  110  strands  of  AWG#30,  substantially  reducing  its  resistance  (from 
0.251  to  0.032)  but  adding  about  0.45  kg  (1  pound)  of  copper  weight.  The  combined  transmitter 
(TXl  -I-  TX2)  shows  a  measured  inductance  of  53/iH  (close  to  the  predicted  value)  and  a  self¬ 
resonance  frequency  at  1.44  MHz  (highlighted  yellow  in  the  upper  right  comer  (of  proprietary 
table)),  well  above  the  intended  operating  bandwidth. 

The  wire  gauge  of  the  three  receiver  coils  has  also  been  changed  from  AWG#30  to  AWG#26 
to  reduce  the  resistance.  The  three  coils,  designed  to  have  similar  area-turn  products,  have  the  fol¬ 
lowing  measured  electrical  parameters  (Table  1).  Small  variations  are  within  acceptable  tolerance 
and  will  be  equalized  through  calibration  procedures,  once  the  sensor  has  been  completed. 

We  assembled  the  sensor  electronics  beyond  the  sensor  head  in  year  2  of  the  project.  Basically, 
there  is  a  single  electronic  console  mounted  at  the  end  of  the  handlebar,  partly  acting  as  a  coun¬ 
terweight  for  handheld  operation.  The  console,  including  a  rechargeable  battery,  communicates 
wirelessly  thru  Bluetooth  with  either  a  PDA  or  a  laptop  that  will  display  data  and  graphics  during 
a  survey  as  well  as  archive  the  raw  data.  The  Bluetooth  also  communicates  with  the  GEM  Beacon 
Navigator  (fabricated  in  Year  2)  to  receive  wirelessly  the  sensor  location  data. 

Eigure  3.2.7  shows  the  overall  system  functional  block  diagram  for  the  GEM-3^+.  The  re¬ 
ceiver  diagram,  Eig.  3.2.8,  utilizes  a  custom-designed  DSP  board  that,  among  others,  performs 
discrete  sine-  and  cosine-convolution  operations  at  each  transmitted  frequency  to  obtain  the  ppm 
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Figure  3.2.4:  GEM-3^'^  sensing  head  and  coil  frames 


Figure  3.2.5:  Beginning  stage  of  wiring.  All  three  receiver  coils  are  now  potted,  electrically  shielded  (silver 
coating),  and  epoxy-painted  for  protection.  Notice  a  small  preamp  printed  circuit  board  (PCB)  for  each 
receiver  coil  taped  to  the  round  z-axis  receiver  coil. 


Benjamin  Barrowes,  Kevin  O’Neill 


-18- 


ENGINEERING 
AT  DARTMOUTH 


MM- 15  37,  GEM-3  ^+-  Final  Report 


3  MATERIALS  AND  METHODS:  GEM-3 HARDWARE 
3.2  GEM-3^+  Continuing  Development 


TX1 


Figure  3.2.6:  Final  drawing  of  the  CRREL  GEM-3^+  sensing  head. 


coil 

mean  coil 

dimensions 

offset  from 
TX  plane 

#  of 

turns 

measured 

resistance 

measured 

inductance 

self-resonance 

frequency 

Hz 

7.25cm  dia  circle 

0 

100 

6.20 

2.78mH 

254kHz 

Hx 

10.3cm  X  14.6cm  reel 

-i-3.5cm 

110 

7.20 

3.07mH 

242kHz 

Hy 

10.4cm  X  14.7cm  rect 

-i-3.5cm 

106 

7.00 

2.92mH 

242kHz 

Table  1 :  Receiver  coils  characteristics 
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responses  of  the  inphase  and  quadrature  eomponents.  Figure  3.2.9  shows  the  transmitter  funetional 
diagram. 


TX 


Figure  3.2.7:  GEM-3^+  electronic  block  diagram 


Figure  3.2.8:  Custom-designed  DSP  module  for  the  GEM-3^+. 


3.3  GEM-3  as  Beacon  Positioning  System 

The  GEM-3^”''  aequires  high  quality  data  from  its  three  reeeiver  eoils.  However,  for  a  hand¬ 
held  instrument  sueh  as  this,  aeeurate  positioning  of  the  sensor  head  ineluding  all  three  Euler 
angles  is  very  important.  One  reason  why  position  is  so  important  is  that  the  depth  to  whieh  the 
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Figure  3.2.9:  Transmitter  block  diagram  with  red  lines  showing  power  flow  and  orange  lines  data  flow. 


can  ’’see”  is  limited  to  approximately  70  cm  for  realistic  size  targets.  At  the  shallow 
depths,  a  positioning  error  on  the  order  of  2  em  ean  eorrupt  data  and  subsequent  inversion  sehemes. 
Aeouraey  of  all  three  orientation  angles  is  important  for  translating  loeal  veetor  eomponents  to  an 
invariant  global  system  during  rotations.  As  a  eonsequenee,  we  expended  signifieant  efforts  at 
developing  and  integrating  the  beaeon  positioning  system  into  the  GEM-3 

Most  currently  used  positioning  systems  rely  on  loeating  something  attaehed  to  the  sensor, 
for  example  the  robotie  total  station,  GPS  reeeivers,  and  laser  positioning  systems.  The  beaeon 
positioning  system  we  have  integrated  into  the  GEM-3^+  works  on  a  different  prineiple  than 
these  types  of  positioning  systems.  It  makes  the  assumption  that  the  primary  field  transmitted 
from  the  sensor  head  is  dipolar  in  shape.  Under  this  assumption,  two  triaxial  reeeiver  eoils  in 
a  fixed  configuration  should  theoretically  allow  the  determination  of  the  location  and  two  of  the 
three  Euler  angles  assoeiated  with  the  sensor  head.  Because  the  dipole  field  is  axially  symmetrie, 
the  third  Euler  angle  which  describes  the  rotation  of  the  sensor  head  around  its  axis  of  symmetry 
cannot  be  deduced  simply  from  the  dipolar  primary  field.  In  faet,  for  the  usual  GEM  3  sensor  head, 
this  last  Euler  angle  is  unnecessary.  However,  beeause  the  GEM-3  has  the  additional  transverse 
receiver  coils  whose  responses  depend  on  this  rotation  about  the  axis  of  symmetry,  some  method 
to  find  this  last  Euler  angle  was  needed.  We  decided  to  install  an  electronic  compass  onto  the  main 
sensor  head.  This  eompass  when  eombined  with  the  data  from  the  beacon  positioning  system,  will 
provide  the  full  set  of  Euler  angles  neeessary  to  uniquely  determine  the  location  and  orientation  of 
the  sensor  head  at  all  times. 

The  beaeon  positioning  system  that  Geophex  demonstrated  to  us  in  2007  was  for  a  different 
Geophex  instrument.  This  other  beaeon  positioning  system  had  a  eonstant  eleetronic  gain  built 
into  the  system  itself.  This  gain  is  an  important  faetor  in  determining  the  sensitivity  and  aeeuracy 
of  the  final  positioning  system.  Too  small  of  a  gain  leads  to  a  short  range;  too  large  of  a  gain  will 
lead  to  saturation  when  the  sensor  is  near  the  beaeon’s  triaxial  reeeivers.  Eor  the  GEM-3^’*",  the 
decision  was  made  to  have  a  variable  gain  in  order  to  have  both  a  large  range  from  the  reeeiver 
eoils  and  aeouraey  when  the  sensor  head  is  in  olose  proximity  to  the  beaeon  reeeivers.  This  added 
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Figure  3.2.10:  Detail  of  the  transverse  receiver  coils  on  the  GEM-  3^+. 
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Positioning 
Receiver  cube  #1 


Positioning 
Receiver  cube  #2 


GEM 

Transmitter 

Figure  3.3.1:  Concept  of  GEM-3^^  positioning  system.  The  (assumed)  dipolar  shape  of  the  GEM-3^"^ 
primary  field  is  recorded  by  two  triaxial  receivers  held  in  a  fixed  relationship. 


complexity  to  the  positioning  system  delayed  its  final  integration  with  the  GEM-3 

Coneeptually,  the  GEM-3^+  beaeon  positioning  system  operates  as  depieted  in  Eigs.  3.3.1 
and  3.3.2.  The  data  from  the  beaeon  system  is  fully  ineorporated  into  the  GEM-3^'''  data  stream. 
One  important  note  about  this  positioning  system  is  that  the  EMI  primary  field  from  the  GEM-3 
penetrates  trees  and  mounds  as  easily  as  it  does  the  soil,  so  this  system  ean  operate  in  ehallenging 
environments  with  non-metallie  obstaeles. 

A  sehematie  of  the  eoordinate  systems  involved  in  the  GEM-3^'*'  beaeon  system  is  shown  in 
Eig.  3.3.3. 


3.3.a  Straight  Line  Tests 

As  an  initial  test  of  the  beaeon  positioning  system,  we  set  up  a  2x10  on  the  ground  at  the  Geophex 
faeility  and  moved  the  GEM-3  along  straight  lines  next  to  the  boards  at  different  distanees  from 
the  beam.  Eigure  3.3.4  shows  the  setup  used  in  this  test,  while  Eigs.  3. 3. 5-3. 3. 8  shows  the  relative 
aeeuraey  of  the  system  at  near,  intermediate  and  far  distanees  from  the  beam.  The  proeessing 
system  infers  the  x,  y,  and  z  position  of  the  sensor  head  as  the  instrument  is  moved  around  the 
plank.  The  merged  data,  EMI  and  positional,  are  then  reeorded  simultaneously.  If  desired,  the 
position  of  the  sensor  head  in  real  time  ean  be  displayed  (see  Eig.  3.3.9).  This  display  eould  be 
useful  for  real  time  deteetion  of  anomalies  as  both  the  EMI  data  and  the  sensor  head  position  are 
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Small  Rx  cubes 
at  ends  of  a 
hollow  plastic 
tube 


Figure  3.3.2:  Beam  holding  the  two  triaxial  receivers  of  the  GEM-3^+  positioning  system. 


Figure  3.3.3:  Coordinate  systems  for  the  target,  beam,  and  GEM-  3^+. 
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Figure  3.3.4:  Beacon  positioning  system  and  plank  setup  for  initial  demonstration  at  Geophex  facility. 


Figure  3.3.5:  Beacon  positioning  system  data  along  a  board  for  the  entire  length. 


overlaid. 

The  details  of  the  beacon  positioning  system  involve  the  nonlinear  inversion  of  a  dipolar  source 
based  on  6  measurement  points.  The  instrument  and  the  two  triaxial  receivers  are  shown  concep¬ 
tually  in  Fig.  3.3.10.  An  oscillating  dipole  magnetic  field  is  generated  by  the  EMI  transmitter.  The 
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Y 


(cm) 


X  (cm) 

Figure  3.3.6:  Beacon  positioning  system  data  along  a  board.  Near  range. 


Y 

(cm) 


X  (cm) 

Figure  3.3.7:  Beacon  positioning  system  data  along  a  board.  Intermediate  range. 
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X 

Figure  3.3.8:  Beacon  positioning  system  data  along  a  board.  Large  range. 


Figure  3.3.9:  Beacon  positioning  system  real  time  feedback  of  EMI  data  (not  shown)  and  position. 
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Two  3-Axis 
Receiving  Coils 


(XI,  Yl,  Zl) 


(X2,  Y2,  Z2) 


Figure  3.3.10:  GEM-3^+  beacon  positioning  system  principles. 
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where  R  =  _|_  ^2  _|_  ^2  dipole-field  point  distanee,  /io  is  the  permeability  of  free  spaee,  the 

subseripts  refer  to  field  eomponenEsouree  dipole  eomponent,  and  the  measurement  (field)  point  is 
at  the  origin.  For  the  loeator  system  eomprising  two  sets  of  eoils  at  two  field  positions  {xa,0,0) 
and  (xh,0,0),  two  sets  of  equations  as  above  are  formed  by  replaeing  x  with  {xa,xb).  We  re¬ 
parametrize  these  equations  in  terms  of  a  total  dipole  strength  and  the  sine  of  the  tilt  angels  in 
plaee  of  mx,y,z  A  non-linear  iterative  inversion  algorithm  is  performed  by  a  eomputer  to  solve  For 
the  x,y,z  and  tilts;  at  initialization  we  also  solve  for  the  dipole  strength,  then  freeze  it  for  stability. 

The  aeeuraey  of  the  beaeon  positioning  system  was  determined  in  tests  at  the  Geophex  faeility 
in  2007.  The  eomplete  setup  ean  be  seen  in  Fig.  3.3.11.  The  sensor  head  was  then  moved  on  a 
preeise  grid  while  the  beaeon  positioning  system  was  reeording  data.  The  ealeulated  position  was 
then  eompared  to  the  aetual  position  on  the  grid  with  results  shown  in  Fig.  3.3.12.  The  horizontal 
projeetion  of  loeator  measured  relative  position  (blue  diamonds)  vs.  template  grid  (red  erosses). 
There  are  50  samples  at  eaeh  position.  The  3-dimensional  RMS  error  is  0.77  em.  Note,  one  grid 
point  is  missing  due  to  data  aequisition  dropout.  Figure  3.3.13  repeats  the  test  exeept  the  nearest 
edge  of  the  gridboard  was  1.5m  away  from  the  sensor  head.  Relative  measured  position  (blue) 
plotted  on  the  ideal  grid  (red)  for  the  template  at  1 .5  m  elosest  edge  from  the  loeater  beam.  There 
are  50  samples  superposed  at  eaeh  point.  The  3-D  rms  error  is  0.92  em. 

More  information  on  the  proeessing  method  used  in  analyzing  dynamie  data  from  the  GEM- 
3^+  beaeon  positioning  system  ean  be  found  in  See.  4.1.e.(3).  A  more  detailed  explanation  on  the 
eonstruetion  of  the  beaeon  hardware  ean  be  found  in  Appendix  C.2. 
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Locator  beam  with  two  3-axes  coil  sets,  3  DSP's,  Bluetooth 


Locator  reference  coil  DSP,  Bluetooth,  computer  interface 


Standard  GEM-3 
Console 


Locator  reference  coil 


Standard  GEM-3 
Sensor  head 


Figure  3.3.11:  Complete  setup  of  the  beacon  positioning  system  tests  at  Geophex. 


3.4  GEM-3^+  Complete 

Figure  3.4.1  show  the  completed  GEM-3^'*^  before  the  stiffening  was  added  (see  below).  Also  see 
Fig.  4.1.1  for  picture  of  the  completed  GEM-3^+  instrument  head. 

The  GEM-3^’*^  instrument  is  now  completed,  including  the  receiver  beam  for  the  beacon  po¬ 
sitioning  system.  Software  which  runs  both  the  beacon  software  and  the  GEM-3^+  itself  has 
been  provided  to  CRREE  by  Geophex.  Setup  and  operation  is  straightforward  and  is  detailed  in 
Appendix  C.l.  We  have  acquired  both  lab  and  test  plot  data  and  run  blind  tests  of  these  data  (see 
Sec.  5).  As  with  most  prototype  instruments,  the  GEM-S^”*"  is  not  a  production  ready  instru¬ 
ment.  There  remain  some  issues  with  software  stability  and  hardware  communication  (between 
the  GEM-3 and  the  beam).  These  issues  cause  software  crashes  and  resets  that  restrict  larger 
deployment.  The  head  deformation  issue  described  above  is  still  relevant  when  the  GEM-3^'''  is 
held  in  a  dynamic  data  acquisition  mode.  Essentially,  when  the  GEM-3^+  is  tilted  beyond  a  nomi¬ 
nal  degree  (beyond  ~  10®),  the  primary  field  can  leak  into  the  transverse  receivers  and  contaminate 
the  data. 

Erequency  domain  instruments  such  as  the  GEM-3^’*^  are  desirable  for  several  reasons  cited 
above  (Sec.  2).  Eessons  learned  from  the  GEM-3^''^  can  assist  future  researchers  in  developing 
more  sensitive  and  robust  ED  instruments  and  beacon  positioning  systems  that  can  potentially 
resolve  these  issues. 
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3.4  GEM-3  Complete 


Figure  3.3.12:  Results  for  beacon  positioning  system  test  #1.  The  nearest  edge  of  the  gridboard  is  Im  from 
the  beam.  The  horizontal  projection  of  locator  measured  relative  position  (blue  diamonds)  vs.  template  grid 
(red  crosses).  There  are  50  samples  at  each  position.  The  3-dimensional  rms  error  is  0.77  cm.  Note,  one 
grid  point  is  missing  due  to  data  acquisition  dropout. 
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3.4  GEM-3  Complete 


Figure  3.3.13:  Results  for  beacon  positioning  system  test  #2.  The  nearest  edge  of  the  gridboard  is  E5m 
from  the  beam.  Relative  measured  position  (blue)  plotted  on  the  ideal  grid  (red)  for  the  template  at  1.5  m 
closest  edge  from  the  locater  beam.  There  are  50  samples  superposed  at  each  point.  The  3-D  rms  error  is 
0.92  cm. 
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Figure  3.4.1:  Final  GEM-3^'*"  instrument. 


4  Materials  and  Methods:  Software  Modeling  and  Analysis  of  GEM-3^+ 
Data 

During  the  course  of  MM  1537,  we  developed  models  and  algorithms  which  aid  in  helping  to 
analyze  not  only  the  data  from  the  GEM-S^'*',  but  potentially  from  any  other  EMI  instrument.  This 
section  will  proceed  as  follows.  Some  general  characteristics  of  the  GEM-3 and  the  beacon 
positioning  system  are  presented  as  Sec.  4.1.  Also  in  that  section  is  a  preliminary  description  of  the 
GSEA  (Generalized  Standardized  Excitations  Approach),  while  a  more  detailed  presentation  is  in 
Sec.  4.5.  In  Sec.  4.2,  we  first  discuss  the  response  of  two  different  soils  to  the  GEM-S^"*"  acquired 
at  a  local  gravel  pit.  We  also  discuss  the  form  of  the  permeability  of  most  soils  at  EMI  frequencies 
and  discuss  implications  for  both  ED  and  TD  instruments.  We  also  model  soil  properties  such 
as  surface  roughness  and  discuss  the  estimation  of  soil  susceptibility  in  Sec.  4.3.  Section  4.4 
documents  our  efforts  to  find  an  absolute  scale  factor  of  GEM  ppm  converting  it  to  magnetic 
field  H  values.  This  can  be  useful  for  local  soil  characterization  and  background  subtraction,  in 
particular. 

We  have  adapted  our  normalized  surface  magnetic  source  (NSMS)  and  GSEA  models  to  the 
GEM-3^+  geometry  (see  Sec.  4.6  and  Sec.  4.5).  These  ultra-high  fidelity,  physically  complete 
models  facilitate  discrimination  processing  by  a)  allowing  detailed  inspection  of  expected  signals 
under  any  contemplated  circumstances;  b)  enabling  signal  pattern  matching  to  establish  the  pres¬ 
ence  of  known  or  likely  target  types;  and  c)  providing  fully  realistic  synthetic  input  data  on  a  large 
scale  for  training  statistical  learning  machines.  In  addition,  we  also  developed  a  model  for  EMI 
responses  by  elliptical  metallic  bodies,  applicable  at  higher  frequencies  by  design,  but  also  match¬ 
ing  fairly  well  with  data  at  low  frequencies  if  they  are  highly  permeable  (see  Sec.  4.7).  This  is  the 
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Figure  4.1.1:  GEM-3^+  sensor  and  measurement  grid 


only  non-numerical  solution  available  for  objects  with  three  dimensional  geometry.  We  have  also 
machined  6  ellipsoids  of  varying  aspect  ratios  to  compared  our  models  with  measured  data  [5]. 
These  models  will  allow  us  to  accommodate  different  size  and  shape  targets  whether  they  are  in 
our  library  for  not. 

4.1  A  vector  handheld  frequency-domain  sensor  for  UXO  identification 

4.1. a  Executive  Summary 

This  section  presents  some  methods  and  results  related  to  UXO  identification  using  the  GEM-3^+. 
Our  analyses  exploit  data  provided  by  the  sensor  in  both  grid-based  and  dynamic  measurements  to 
characterize  different  objects,  including  metal  spheres  and  actual  UXO.  For  the  data  analysis  we 
alternate  between  the  dipole  model  and  the  more  rigorous  standardized  excitation  approach.  We 
review  some  ill-conditioning  issues  encountered  with  the  latter  model  and  the  different  approaches 
that  we  use  to  overcome  them.  In  applications,  the  availability  of  horizontal  field  components  in 
the  data  allow  us  to  identify  UXO  vs.  non-UXO  items  while  minimizing  the  nonnegligible  effects 
of  ground  response. 

4.1. b  Introduction 

The  improvement  in  sensing  methods  represented  by  the  GEM-3  must  be  accompanied  by  fast, 
physically  complete,  and  accurate  forward  models  that  can  both  reasonably  pinpoint  the  location 
and  orientation  of  an  object  and  provide  sufficient  information  on  its  electromagnetic  properties  so 
as  to  characterize  it  unambiguously.  Here  we  adapt  to  the  GEM-3  some  forward  models  that 
have  shown  success  with  other  sensors.  Perhaps  the  most  widely  used  model  treats  each  anomaly 
as  a  point  dipole  [6-10];  this  procedure,  while  limited,  is  fast  and  provides  good  starting  guesses 
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for  more  elaborate  methods.  An  example  of  the  latter  is  the  data-derived  Standardized  Exeitation 
Approaeh  (SEA),  whieh  ineorporates  the  finite  size  of  the  seatterer  and  the  nonuniformity  of  the 
primary  field  in  the  relevant  length  seale.  The  approaeh  deeomposes  the  primary  field  striking  the 
target  into  a  set  of  standardized  modes,  whose  relative  importanee  in  the  expansion  is  given  by 
the  partieulars  of  position  and  orientation.  The  objeet  is  replaeed  by  rings  of  magnetie  sourees 
plaeed  on  a  surrounding  spheroid.  The  response  of  the  rings  to  eaeh  mode  of  exeitation  ean  be 
deseribed  by  a  Green  funetion,  and  the  full  response  of  the  target  is  a  superposition  of  these  with 
the  eorresponding  geometrie  eoeffieients  as  weights.  Onee  the  objeet’s  response  to  eaeh  mode  is 
established,  the  response  ean  be  ealeulated  for  any  variation  of  exeitation  and  reeeiver  loeation  by 
just  using  an  appropriately  altered  superposition  of  the  modal  responses.  The  method  was  previ¬ 
ously  adapted  to  the  GEM-3  [11]  and  to  various  time-domain  instruments  [12,  13]  using  eharges  as 
the  responding  sourees;  a  generalization  replaees  the  eharges  by  dipoles  and  uses  Stokes’s  theorem 
to  streamline  the  neeessary  integration  over  the  reeeiver  areas  [14,  15]. 

The  overall  seetion  is  organized  as  follows.  After  deseribing  the  instrument  in  Seetion  4.1.e  we 
present  some  studies  related  to  UXO  identifieation  that  use  data  produeed  by  the  GEM-3^”''in  both 
grid-based  and  dynamie  measurements.  We  employ  both  the  dipole  model  and  the  SEA,  whieh  we 
introduee  in  seetions  4.1.d.(l)  and  4.1.d.(2)  respeetively,  to  eharaeterize  different  objeets,  ineluding 
metal  spheres  and  aetual  UXO.  We  look  at  the  possibility  of  using  the  horizontal  eomponents  of  the 
seeondary  field  to  identify  UXO  vs.  non-UXO  items  while  minimizing  the  nonnegligible  effeets 
due  to  the  ground.  We  eonelude  in  Seetion  4. 1  .f. 

4.1.C  The  GEM-3  Sensor 

The  GEM-3^~^  sensor,  designed  and  manufaetured  by  Geophex,  Etd.,  is  a  wideband  frequeney- 
domain  instrument  that  eontinuously  transmits  ehosen  frequeneies  and  outputs  ED  speetra  of  re¬ 
sponse.  The  user  ean  seleet  up  to  15  frequeneies  of  operation  for  the  sensor,  ranging  between 
30  and  90030  Hz.  The  measured  seeondary  field  at  eaeh  frequeney  is  given  in  dimensionless  units 
of  parts  per  million  (PPM),  defined  as  [16,  17] 

Seeondary  magnetie  flux  through  the  reeeiver  eoil 

PPM  =  — - — — ^7 - r - r - q-xl0^  (4.1.1) 

Primary  magnetie  flux  through  a  referenee  eoil 

whieh  serve  to  divide  out  the  frequeney  dependenee  and  the  value  of  the  transmitted  eurrent  and 
aeeount  for  the  faet  that  the  aetual  measured  quantities  are  voltages  proportional  but  not  identieal 
to  magnetie  fluxes. 

The  sensor  head  eontains  two  eireular  eoplanar  eoneentrie  transmitting  eoils  eonneeted  in  se¬ 
ries.  These  transmitters  have  radii  27  em  and  10.84  em  and  are  wound  6  and  2  turns  respeetively, 
making  the  effeetive  eurrent  on  the  outer  eoil  three  times  larger  [18].  This  “transmitter  bueked” 
arrangement  is  designed  to  produee  no  primary  magnetie  flux  through  a  reeeiver  eoil,  eoplanar  and 
eoneentrie  with  the  transmitters  and  of  radius  7.25  em,  to  whieh  we  heneeforth  refer  as  the  H^-coil 
(these  numbers  are  determined  by  using  an  expression  equivalent  to  (4.1.2)  below  but  oast  as  an 
expansion  in  Eegendre  polynomials  [16]).  The  H^-coil  ean  then  measure  the  mueh  smaller  flux  due 
to  the  seeondary  field. 

The  other  two  reeeiving  eoils  on  the  GEM-3  are 

1.  The  “iTx”  coil  with  dimensions  10.3  em  (in  z)  by  14.6  em  (in  y),  offset  3.5  em  from  the  plane 
of  the  transmitter  eoils; 
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Figure  4.1.2:  (Left)  The  field  due  to  the  GEM-3^'*'  slightly  above  the  plane  of  the  transmitter  coils.  The 
central  cavity  where  the  primary  flux  vanishes  is  clearly  visible.  (Right)  The  field  due  to  the  GEM-3^+ 
on  the  plane  of  the  transmitter  coils,  this  time  seen  in  profile. 


2.  The  “Hy”  coil  with  dimensions  10.4  cm  (in  z)  by  14.7  cm  (in  x),  offset  3.5  cm  from  the  plane 
of  the  transmitter  coils. 


The  primary  flux  through  these  vertical  receivers  vanishes  trivially. 

The  primary  field  due  to  a  circular  coil  of  radius  X  in  cylindrical  coordinates  is  [19,  20] 


H 


i-'-kAm-Eik) 


where  the  modulus  k 


4pZ 

{X+p)^+z^' 

(4.1.2) 


where  K{k)  and  E{k)  are  elliptic  functions  [21] 

The  field  slightly  above  the  plane  of  the  transmitter  coils  is  depicted  in  three  dimensions 
and  as  a  function  of  p  in  Figure  4.1.2.  This  figure  also  shows  how  the  total  field  arises  as  the  sum 
of  the  contributions  of  the  two  transmitter  coils. 


4.1.c.(l)  Integrating  over  the  areas  of  the  receivers  Consider  an  upright  dipole  located  zq  ni  below 
the  origin.  With  m  =  mz,  r  =  pp,  and  r'  =  —zqz  we  can  use  the  general  formula 


H 


3  (m  •  ]R)M/i?^  —  m 
4-7tR^ 


m 


to  obtain  iT  =  — 

4;r  (z2  +  p2)5/2 


Integrating  (4.1.3)  over  the  area  of  the  horizontal  receiver  we  find 

A 


m 


d>  =  2;r  /  H^pdp  =  — 


a 


10 


2  (zQ  +  a2)3/2 


(4.1.3) 


(4.1.4) 


On  the  other  hand,  if  we  substitute  the  field  at  the  center  and  multiply  by  the  receiver  area,  we  get 
the  “average”  flux 


ml  2 
- - ^Tta 


m 


(4.1.5) 
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We  see  that  the  average  flux  is  larger  than  the  integrated  eounterpart  when  we  are  right  above  the 
target: 


4>' 


1  + 


3/2 


(4.1.6) 


For  fixed  a  this  effeet  grows  as  zq  deereases;  i^,  when  the  target  gets  eloser  to  the  sensor. 

We  ean  also  let  vary  the  horizontal  position  of  the  target.  We  take  r'  =  xqx  —  zqz,  and  we  obtain 
for  the  field 


^  _  m  2zQ-p2-XQ  +  2pxocos^ 

^  {zl  +  +xl-2pxocos(j)yf^ 


(4.1.7) 

co-rp  -taq- 

There  is  probably  a  elosed  expression  for  the  flux,  but  it  is  undoubtedly  very  eomplieated.  An 
intermediate  expression  is 


m 

d>  =  -  — 


z^k^  p  dp 


SnJo  (p|xo|)3/2(p-|.ro|)^+Zo) 


4iiP-\M)^+4) 

(4.1.8) 

where  K{k)  and  E{k)  are  the  eomplete  elliptie  integrals  of  the  first  and  seeond  kinds  and  the 
modulus 


k^  = 


4p\xo\ 


(p  +  hl)"+4' 

Making  xq  =  0  in  (4.1.8)  we  easily  reeover  (4.1.4).  If  we  merely  substitute  the  field  at  the  center 
and  multiply  by  the  receiver  area,  we  obtain 


m  Izl-xl  2 
d>  = - ^  na  = 


0  -^0 


(4.1.10) 


4;r  (z2+;c2)5/2  4  ^^2  ^^2^5/2 

As  part  of  the  Generalized  Standardized  Excitation  Approach  to  be  considered  later,  it  is  pos¬ 
sible  to  use  a  shortcut  in  this  integration  by  employing  Stokes’s  theorem: 


A-dl  = 


c 


^(VxA) 


hds 


(4.1.11) 


Defining  a  “vector  potential”  (which  in  rigor  should  be  proportional  to  the  permeability  of  the 
medium)  by 

m  X  M 

A  = - ^  (4.1.12) 

we  can  easily  verify  that  its  curl  is  given  by  (4.1.3).  For  the  situation  considered  in  this  section  we 
have 

mp0  .  _  j  *  maq) 


A  = 


in  general  and  A  = 


47t{p^  +Zq)^/^  4;r(a2 -)-Zq)^/2 

along  the  circumference  of  the  receiver,  which  trivially  integrates  to  yield  (4.1.4). 


(4.1.13) 


4.1.c.(2)  The  other  receivers  To  verify  that  our  integration  routines  are  correct  we  compare  the 
integrated  flux  to  the  average  one  for  an  upright  dipole.  In  this  case  the  depth  is  zq  =  —35  cm,  and 
we  have  xq  =  4  cm  and  yo=  12  cm.  Interesting  results:  The  plots  are  very  similar;  the  plots 
show  the  same  tendency  we  saw  in  the  last  section — the  averaged  flux  tends  to  be  larger  than  the 
integrated  one;  the  Hy  plots  show  the  opposite  tendency.  See  Figure  4.1.3.  The  dipole  can  also  be 
tilted,  as  Figure  4.1.4  shows. 
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Figure  4.1.3:  Integrated  (blue)  and  averaged  (green)  flux  profiles  on  a  grid  above  an  upright  dipole;  from 
left  to  right  we  have  the  z-,  x-,  and  y-components.  The  red  is  the  difference  between  the  two. 


4,l.c.(3)  Beacon  positioning  One  of  the  features  that  make  the  GEM-3^”*'  potentially  useful  for 
dependable  UXO  identification  is  its  positioning  system,  which  provides  sub-centimeter  accuracy 
and  can  be  used  in  treed  or  rugged  terrain.  The  beacon  is  also  lightweight  and  unobtrusive,  requir¬ 
ing  almost  no  hardware  to  be  attached  to  the  sensor  itself. 

For  the  data  analysis  it  is  necessary  to  have  the  tilt  information  referred  to  the  same  global 
coordinate  system  in  which  the  sensor  location  is  known.  This  cannot  be  inferred  directly,  however, 
since  the  field  of  a  dipole  is  azimuthahy  symmetric:  the  GEM-3 can  be  rotated  about  its  z-axis 
and  appear  identical  to  the  beacon  system.  To  ameliorate  this,  the  sensor  head  also  has  affixed 
to  it  a  digital  compass  [22]  that  at  every  sensor  location  gives  the  three  Euler  angles,  in  the  yaw- 
pitch-roh  convention  [23],  in  a  coordinate  system  defined  by  magnetic  North  and  the  acceleration 
of  gravity. 

The  connection  between  the  compass  and  beacon  coordinate  systems  can  be  determined  in  a 
least-squares  sense  from  the  fact  that  the  direction  of  the  dipole  moment  (z  in  a  reference  frame 
fixed  to  the  sensor  head)  is  unique.  We  thus  find 

RN^b=[H  H  ]  [4  4  ■■■  4]^  (4.1.14) 


where  the  dagger  denotes  the  pseudoinverse[24].  This  rotation  matrix  can  be  polished  further  (i.e., 
made  “more  orthogonal”)  by  iterating[25] 
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H 


H 


X 


Figure  4.1.4:  Integrated  (blue)  and  averaged  (green)  flux  profiles  on  a  grid  above  an  ^-directed  dipole;  from 
left  to  right  we  have  the  z-,  x-,  and  y-components.  The  red  is  the  difference  between  the  two. 


with  Rq  being  the  matrix  found  from  (4.1.14).  Very  few  iterations  (rarely  more  than  three)  are 
needed. 


4.1.d  Forward  models 


4.1.d.(l)  The  dipole  approximation  A  popular  and  useful  model  for  UXO  discrimination  re¬ 
places  the  buried  target  by  one  or  more  triaxial  point  dipoles,  either  in  the  frequency  domain  [7,  26] 
or  in  the  time  domain  [9].  The  dipole  model  has  been  found  to  fit  measured  data  adequately  when 
bodies  are  small  enough  that  the  primary  field  can  be  assumed  uniform  along  their  extent,  when 
they  are  far  enough  away  from  the  sensor  that  higher  moments  vanish,  and  when  their  composi¬ 
tion  is  homogeneous  enough  that  there  are  no  conflicting  signals  from  different  parts  of  the  ob¬ 
ject  [8,  27]. 

The  signal  picked  up  by  the  /f  th  receiver  is  when  the  sensor  is  centered  at  is 


Sfiirj)  = 


'Ruitj 


3{m-R)R/R^-m 


nil  dsR. 


(4.1.16) 


where  i?  =  r  —  r',  r  is  any  observation  point,  and  r'  is  the  location  of  the  dipole.  The  dipole 
moment  m  =  An  example  fit  to  actual  UXO  data  using  this  model,  as  well  as  the  extracted 
polarizability  elements,  can  be  found  in  Fig.  4.1.5. 

We  can  use  the  results  given  by  this  approximation  as  a  starting  guess  for  the  more  demanding 
SEA  approach,  which  we  discuss  next. 
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Frequency  (Hz) 

Figure  4.1.5:  (Upper)  Quadrature  part  of  the  flux  through  the  three  receiver  coils  for  an  H-83  UXO  with 
nose  tilted  45°  down.  (Lower)  Extracted  polarizability  elements  (dashed  lines:  inphase  part;  solid  lines: 
quadrature  part);  it  is  clearly  seen  that  the  object  exhibits  cylindrical  symmetry. 


4.1.d.(2)  The  Generalized  Standardized  Excitation  Approach  The  Generalized  Standardized  Ex- 
eitation  Approaeh  (GSEA)  [14,  15]  is  simultaneously  a  generalization  and  a  simplifieation  of  the 
usual  SEA[11].  Eor  determining  a  superior  set  of  equivalent  soureesto  produee  the  fundamen¬ 
tal  modal  responses,  its  starting  point  is  the  observation  that  in  the  EMI  regime  the  fields  out¬ 
side  a  metallie  target  are  irrotational  and  ean  henee  be  deseribed  by  a  sealar  potential  that  obeys 
the  Eaplaee  equation,  whose  fundamental  solution  in  a  prolate  spheroidal  system  (— 1  <  t]  <  1, 
l<,^<oo,  0<<^<  In)  eentered  at  the  seatterer  is  the  superposition  [28] 


V/(r) 


-j  CX3  CX3  1 

E  £  where 

m=Qn=m p=0 


oosm0,  p  =  0 
sinm0,  p  =  1 

(4.1.17) 


and  the  E™(-)  are  assoeiated  Eegendre  funetions.  The  interfoeal  distanee  for  a  prolate  spheroid 
with  semiminor  and  semimajor  axes  a  and  b  is  d  =  2\/b^  —  a^;  the  surfaee  of  the  spheroid  is  the 


Benjamin  Barrowes,  Kevin  O’Neill 


-39- 


engineering 

AT  DARTMOUTH 


MM-1537,  GEM-3^+-  Final  Report  4  MATERIALS  AND  METHODS:  GEM-3^+  SOFTWARE 

4.1  A  vector  handheld  frequency-domain  sensor  for  UXO  identification 


set  of  points  ^  =  ^Q  =  2b/d.  We  ehoose  the  seale  faetor  Hq  =  1.  It  is  easier  and  more  efficient  to 
perform  the  decomposition  on  the  normal  magnetic  field  over  the  spheroid,  employing  the  fact  that 


(r)  =  — ^ •  Vt//(r)  =  —  ,  where  the  metric  coefficient  = 

^  oq  ^ 

The  orthogonality  of  the  Legendre  and  trigonometric  functions  [29]  yields 


;o^-i  ■ 
(4.1.18) 


rlTI 


/o 


■'pmn 


V-.T//.2  ^  {n  +  m)\dP-{^) 

2  ^°^2n+l(n-m)! 


^0 


(4.1.19) 


where  a  =  2forp  =  m  =  0  and  a  =  I  otherwise.  To  compute  (4.1.19)  we  use  for  each  transmitter 
coil  the  exact  expression  (4.1.2)  for  the  field  of  a  circular  current  and  evaluate  the  integrals  using 
Gauss-Legendre  quadrature  and  the  periodic  trapezoid  rule[30].  The  signal  in  the  /ith  receiver 
when  the  sensor  is  at  rj  is  given  by 


oo  oo 

(*■;)=  HHipl 

m=Qn=ml=l 


(r  ■) 


(4.1.20) 


where 

^  H  bpmniyj) 


1 

E 

p=o 


/Rq(rj 


rlTt  3  (£L 

ds^p  hp  ■  /  Tp^{(j)') - 

I  Jo 


(4.1.21) 


with  M/  =  r  —  r^I  and  the  position  vector  on  the  Zth  ring 


r'/ =  Z/(*cos())'  +  ^sin())')  +  f  and  Xi  =  ^\/(l  -  -  !)•  (4.1.22) 

In  other  words,  instead  of  the  rings  of  charge  considered  before  [11,  13]  we  use  rings  of  radially 
directed  dipoles.  This  is  a  simplification  of  the  original  SEA  in  that  it  is  not  necessary  to  force  the 
total  charge  to  vanish  in  order  to  comply  with  the  Maxwell  equations.  Equivalently, 


~  H  bpmniTj 


p=0 


d^Rp  • 


i'xRi 

AtzR] 


Xid^'. 


(4.1.23) 


Equation  (4.1.20)  can  be  turned  into  a  linear  system  by  finding  an  appropriate  invertible  map- 


ping  {l,mn)  <-»•  k 

S'  = 

[•SKri)  ■■ 

S'xi’^N,)  Sl(ri) 

■■■ 

(4.1.24) 

and 

r  = 

■n,i(ri) 

rj,^,i(ri) 

rj,,i(r2)  ... 

(4.1.25) 
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Once  we  have  this,  we  minimize  the  diserepaney  to  obtain 

minUrp-S'Yirp-S').  (4.1.26) 

p 

This  is  a  highly  ill-eonditioned  system  that  requires  regularization.  We  have  employed  two  regu¬ 
larization  methods:  the  truncated  Singular  Value  Deeomposition  (TSVD)  [31]  and  the  Tikhonov 
procedure  [24]. 


4.1.e  Results 


We  start  by  stressing  the  neeessity  of  integrating  over  the  reeeivers  when  ealeulating  responses 
for  the  ease  of  monostatie  measurements.  We  eonsider  a  sphere.  It  is  well  known[32]  that  a 
sphere  of  eonduetivity  o,  relative  permeability  and  radius  a  plaeed  in  a  uniform  primary  field 
behaves  like  a  dipole  of  moment  m  =  for  small  spheres  we  ean  assume  HP'"  from  (4.1.2)  to 
be  approximately  uniform  over  the  target  size,  and  the  polarizability  matrix  B  =  /3 1,  where 


/3  =  Ina^ 


(2pr+  l){kacothka  —  1)  —  {kaY 
(/i^  —  l)(fcaeothA:a  —  1) -h  (A:a)2  ’ 


and  k=  {I +  i)/d  =  ^/ipopYYoo . 


(4.1.27) 


Equation  (4.1.6),  whieh  is  valid  only  right  above  a  target,  still  applies  in  this  ease  and  lets  us 
eheek  the  eorreetness  of  our  results.  The  differenees  are  dramatie  when  the  target  is  very  elose  to 
the  sensor,  as  Figure  4.1.6  illustrates.  In  this  ease  the  depth  is  10  em,  and  the  sphere  is  below  the 
origin.  Equation  (4.1.6)  is  satisfied  to  15  figures.  It  is  evident  that  in  the  interest  of  realism  we 
must  integrate  over  the  reeeiver  areas  when  modeling  responses. 

We  now  turn  to  measured  data  taken  over  aetual  spheres  at  three  different  sensor  heights  (16, 
21,  and  26  em  above  the  surfaee  of  the  sphere).  We  use  two  of  the  data  files  to  extraet  the  dipole 
moments  of  the  rings  and  use  these  to  prediet  the  field  in  the  unseen  instanee.  We  have  used 
two  methods  to  regularize  (4.1.26):  1)  We  first  perform  Tikhonov  regularization  [11,  24]  using 
the  E-eurve  method  [31]  to  determine  the  regularization  parameter.  2)  We  also  eomputed  the 
Singular  Value  Deeomposition  (SVD)  of  truneated  it  using  the  eriterion  suggested  by  the 
Pieard  plot  [31],  and  used  the  truneated  version  to  solve  the  normal  equation  (4.1.26). 


4.1.f  Conclusion 

In  this  seetion  we  have  presented  a  handheld  frequeney-domain  sensor  whose  wide  range  of  oper¬ 
ating  frequencies,  three-dimensional  veetor  reeeptivity,  and  innovative  positioning  system  make  it 
a  strong  eandidate  for  providing  robust  and  reliable  diserimination  of  UXO. 

After  introdueing  the  GEM-3^'''  and  seeing  the  importanee  of  obeying  Faraday’s  law  as  ex- 
aetly  as  possible  when  modeling  the  signal  eolleeted  by  eaeh  reeeiver,  we  proeeeded  to  use  the 
point  dipole  model  and  the  Generalized  Standardized  Exeitation  Approaeh  to  model  synthetie  and 
measured  data. 
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Figure  4.1.6:  Integrated  (blue)  and  averaged  (green)  flux  profiles  on  a  grid  above  a  sphere;  from  left  to  right 
we  have  the  z-,  x-,  and  y-components.  The  red  is  the  difference  between  the  two. 
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Figure  4.1.7:  (Upper)  “Training”  field  used  to  extract  the  GSEA  dipole  moments  of  an  aluminum  sphere; 
(Lower)  Picard  plot  showing  the  singular  values  of 
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4.2  Soil  Response  of  the  GEM-3 

Our  April  2007  White  paper  for  the  GEM-3  and  the  MPV  instruments  details  the  theory  behind 
soil  responses  for  both  instruments.  In  this  seetion,  we  will  diseuss  some  of  the  data  we  have  taken 
from  loeal  soils  with  the  GEM-3  and  the  advantages  afforded  by  the  transverse  reeeiver  eoils 
is  their  relative  immunity  to  eontamination  from  geologieal  noise  due  to  ground  response. 

We  arranged  to  test  some  soils  at  a  loeal  gravel  pit,  Eebanon  Crushed  Stone.  The  main  soil 
we  tested  was  known  as  ’’sump  sand”  and  it  exhibited  a  large  magnetie  response  in  the  frequeney 
domain  to  the  GEM-3^+.  We  aequired  height  profiles  at  2  in.  inerements  up  to  36  in.  over 
this  soil.  Our  measurement  set  up  ean  be  seen  in  Eig.  4.2.1.  We  also  took  high  profiles  over  the 
soil  after  burying  a  steel  sphere  northwest  of  the  measurement  position  (see  Eig.  4.2.2).  Height 
profiles  at  different  frequeneies  after  baekground  subtraetion  are  shown  in  Eig.  4.2.3.  As  we  have 
reported  elsewhere,  the  inphase  part  of  the  GEM-3^+  data  displays  a  trend  of  inereasing  signal 
with  deereasing  height  until  a  maximum  at  around  5  to  10  em  above  the  ground.  The  quadrature 
eomponent  of  the  heights  profile,  on  the  other  hand,  displays  only  a  weak  trend  if  at  all. 


Figure  4.2.1:  Setup  of  GEM-3^+  at  Lebanon  Crushed  Stone  over  “sump  sand”. 

Beeause  the  primary  field  is  oriented  downward  in  the  —z  direetion,  the  prineipal  response  from 
the  soil  is  oriented  in  z  with  only  a  seeondary  response  in  the  transverse  reeeivers  (see  Fig.  4.2.4). 
This  effeet  ean  be  seen  by  eomparing  the  data  from  the  height  profiles  with  and  without  this  sphere 
present  (see  Figs.  4.2.5  and  4.2.6).  With  whatever  variations  in  magnitude,  the  isotropie  response 
from  the  sphere  target  should  be  the  same  in  all  veetor  eomponents.  Notiee  in  Fig.  4.2.6  top, 
that  the  apparent  response  is  shifted  downward  (e.f.  bottom  plot)  as  it  also  eontains  the  soil 
response.  These  figures  show  the  strong  broadband  soil  response  in  the  z  reeeiver  eoil  and  the 
almost  negligible  response  in  the  transverse  reeeiver  eoils  (x  and  y). 

Using  the  seale  faetors  eontained  in  Fig.  4.2.7  and  derived  from  the  eomparison  of  analytieal 
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models  with  actual  sphere  data  (see  Sec.  4.4),  we  were  also  able  to  estimate  the  susceptibility, 
of  this  soil.  A  plot  of  the  match  between  model  and  data  is  provided  in  Fig.  4.2.8.  The  estimated 
susceptibility  of  this  soil  is  =  0.0032,  a  very  realistic  value  based  on  the  work  of  Janet  Simms  at 
the  Standardized  Test  Sites  [33-35]. 

Using  this  technique,  we  hope  to  be  able  to  account  for  local  geological  noise  at  remediation 
sites  by  acquiring  height  profile  data  over  local  soil  which  is  not  contaminated  by  metallic  targets; 
and  also  by  inferring  soil  effects  in  the  component  by  analyzing  the  relatively  uncontaminated 
horizontal  components. 
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Inphase  Hz  (Cleaned) 


frequency  =  30  Hz  frequency  =  450  Hz 


frequency  =  13590  Hz  frequency  =  90030  Hz 


Quadrature  Hz  (Cleaned) 


frequency  =  30  Hz  frequency  =  450  Hz 


frequency  =  13590  Hz  frequency  =  90030  Hz 


Figure  4.2.3:  Inphase  and  quadrature  components  of  GEM-3^+  taken  over  “sump  sand”  with  height  de¬ 
creasing  left  to  right. 
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Sensor 

head 


Figure  4.2.5:  GEM-S^"^  setup  comparing  “sump  sand”  data  to  “sump  sand”-i-sphere  data. 
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Hz  =  Vertical  component 


Frequency  (Hz) 


Hy  =  Horizontal  Component 


Figure  4.2.6:  GEM-3^+  inphase  response  showing  the  substantial  soil  response  contaminating  the  re¬ 
ceiver  but  not  the  Hy  receiver. 
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Figure  4.2.7:  Soil  model  match  to  GEM-S^*^  data  using  scale  factor. 
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Figure  4.2.8:  Soil  model  match  to  GEM-S^*^  data  using  scale  factor. 
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See  additional  material  relative  to  soil  response  -  addendum  to  “White  Paper”  for  MM1537  [36] 
that  illustrates  the  use  of  the  eonversion  faetors  for  inferenee  of  geophysieal  eonelusions,  and  for 
eharaeterization  of  instruments  relative  to  UXO  diserimination. 
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The  problem  with  the  view  encouraged  by  Fig.  4.2.8  is  that  it  does  not  explain  the  fact  that 
significant  soil  responses  are  seen  by  time  domain  instruments.  To  understand  the  issue,  consider 
a  FD  expression  for  in  terms  of  magnetic  susceptibility  X  as 

^r{co)  =  \+x{0i)  =  1 +/(«)  + a"  (to)  (4.2.1) 

where  i  is  the  square  root  of  negative  one  and  the  angular  frequency  CO  is  In  times  frequency  in 
Hz.  The  implication  of  spectral  patterns  like  that  in  Fig.  4.2.8  is  that 

X'  (go)  ~  constant,  x”  ( w)  ~  0  (4.2.2) 

Such  a  constant,  real  valued  x  should  produce  no  response  in  time  domain  measurements.  Given 
negligible  electrical  conductivity  in  a  soil,  its  response  would  simply  be  proportional  to^^'times  the 
magnitude  of  the  transmitted  field  at  any  given  point  in  time,  with  no  delay.  All  response  would  be 
instantaneous.  When  the  transmitted  field  from  a  TD  instrument  is  shut  off  during  the  recording  of 
data,  the  soil  response  would  also  shut  off,  producing  no  data  other  than  noise.  However,  many  time 
domain  measurements  do  indeed  show  some  soil  response  in  more  or  less  orderly  decay  patterns. 
Evidently  some  kind  of  relaxation  effect  takes  place,  beyond  the  instantaneous  response. 

Overall,  one  might  hypothesize  a  soil  response  in  the  form  of  the  schematic  in  Fig.  4.2.9,  given 
a  downward  step  input  in  the  primary  field  at  t  =  0. 


Schematic  of  ground  response 
sequence  in  time,  with  Tx  shutoff  at  t  =  0 


Normalized  Time 


Figure  4.2.9:  Hypothesized  relaxation  type  time  response  for  real  soil  permeabilities. 

The  instantaneous  response  of  the  system  is  from  its  initial  value  of  Uo  to  ui .  This  response  per¬ 
sists  for  all  t>0.  It  corresponds  to  the  flat  real-valued  (i.e.  inphase)  impulse  response  in  Fig.  4.2.9. 
On  top  of  that  is  a  relaxation  to  the  ultimate  response  to  u^o- 

The  only  way  to  explain  measurable  TD  instrument  response  to  the  same  soil  as  in  Fig.  4.2.8  is 
to  hypothesize  that  a  relaxation  type  spectrum  (see  MM  1537  white  paper)  lies  within  the  spectrum 
shown  in  the  Fig.  4.2.8,  but  that  it  is  overshadowed  by  the  larger  instantaneous  response,  i.e.  it  is 
below  the  resolution  of  this  data. 
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To  investigate  this,  we  also  eomputed  eonversion  faetors  for  the  TD  MPV  deviee  (see  See.  4.4). 
If  the  eonversion/ealibration  faetors  for  both  the  GEM  and  the  MPV  are  eorreet,  then  when  the 
primary  fields  are  sealed  so  that  a  sphere  is  illuminated  by  the  same  intensity  from  eaeh  deviee, 
eaeh  instrument  will  report  the  same  magnetie  response  at  the  same  distanee.  Figure  4.2.10  below 
shows  that  this  is  the  ease,  when  the  GEM-3  speetrum  is  eonverted  to  a  TD  response  via  Wait’s 
analytieal  solution  [37,  38]. 


Figure  4.2.10:  Measured  GEM-3  (“converted”)  and  MPV  responses  from  a  sphere,  observed  at  the  same 
distance,  when  each  response  has  been  scaled  additionally  by  the  primary  field  from  each  instrument. 


The  implieation  of  all  this  is  apparent  if  one  applies  the  same  kind  of  sealing  as  in  Fig.  4.2.10 
to  measured  soil  responses.  For  a  CRREF  baekyard  soil  with  FD  speetra  like  that  in  Fig.  4.2.8 
above,  we  also  measured  the  MPV  response.  We  then  eonverted  the  GEM-3  soil  response  to  a 
TD  response,  retaining  the  magnitude  indieated  by  our  eonversion  faetor,  with  observation  points 
about  the  same  distanee  from  the  soil.  When  the  GEM-3  and  MPV  soil  responses  are  further  sealed 
so  that  the  illumination  of  the  soil  surface  is  the  same  magnitude  for  both  instruments,  one  obtains 
Fig.  4.2.11  below. 

That  is.  Fig.  4.2. 1 1  shows  approximately  the  same  kind  of  thing  as  Fig.  4.2. 10.  The  latter  shows 
that,  given  our  conversion  factors,  we  can  illustrate  that  a  given  metal  object  (sphere)  produces  the 
same  response  to  the  same  to  the  same  intensity  of  illumination,  when  observed  from  the  same 
distance  by  different  instruments.  The  reason  that  the  similarly  scaled  soil  responses  in  Fig.  4.2. 1 1 
are  different  is  that  the  TD  instrument  only  sees  part  of  the  soil  response,  namely  the  relaxation 
part.  The  instantaneous  part  has  passed  by  the  time  the  first  time  gate  falls  in  the  MPV.  On  the 
other  hand,  the  GEM-3,  and  for  that  matter  any  other  frequency  domain  instrument,  is  dominated 
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Figure  4.2.11:  Comparison  of  measured  responses  of  the  soil  to  the  GEM-3  and  from  the  MPV,  from  about 
the  same  observation  point,  when  the  illumination  of  the  soil  surface  has  been  scaled  to  be  the  same  for  each 
instrument.  The  red  GEM-3  response,  converted  to  TD,  is  flat  before  the  step  shut  off  at  time  zero,  at  which 
point  it  theoretically  goes  to  zero,  instantaneously.  The  blue  MPV  data  start  some  finite  time  after  shutoff 
(~  10^^  s)  and  decay  thereafter. 


by  the  much  larger  (>  3  OM)  instantaneous  response,  that  appears  at  all  frequencies. 

With  all  this  in  mind,  the  most  general  representation  of  soil  susceptibility  may  be  something 
like: 


X{(0)=XoI+Xo2R{C0) 


(4.2.3) 


where  Xoi  and  Xo2  are  constants  and  R{co)  is  something  like  the  commonly  used  expression: 


R{co)  =  1 


_ ^ 

in{%2h\)  \i(OTi  +  lJ 


(4.2.4) 


where  Ti  and  T2  are  chosen  relaxation  time  limits.  The  parameter;t^oi  would  account  for  the  larger 
flat  response  in  the  figs  above.  2^02 would  be  the  separate  magnitude  of  the  relaxation  response.  In 
this  view,  all  the  FD  data  treated  here  reflect  the  fact  that  Xo2 

A  little  more  explanation  in  terms  of  relaxation  times.  The  schematic  below  in  Fig.  4.2.12 
shows  hypothetical  relaxation  responses  by  different  constituents  within  some  soil. 

The  commonly  used  expression  in  (4.2.4)  corresponds  essentially  to  the  collection  of  blue 
curves.  That  model  assumes  that  there  is  a  continuous  distribution  of  relaxation  times,  evenly 
distributed  in  log  space,  and  that  all  the  curves  converge  to  the  same  magnitude  of  contribution. 
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Figure  4.2.12:  Hypothetical  relaxation  curves  in  (log)  time  for  different  constituents  of  some  soil. 


This  has  been  said  to  model  the  TD  responses  effectively  for  many  soils  encountered  in  UXO 
studies  [39-41].  What  our  research  and  conversion  factor  applications  show  is  that  an  additional 
component  with  very  short  relaxation  time  is  also  typically  present,  converging  to  a  much  larger 
magnitude  response  (red  curve/line).  Treating  the  fast  response  as  effectively  instantaneous  simply 
means  adding  the  term  Xoi  in  (4.2.3). 

It  has  been  suggested  that  one  might  cover  all  the  effects  mentioned  here  while  retaining  the 
older  soil  response  model,  i.e.  (4.2.3)-(4.2.4)  with  Xoi  =  by  stretching  the  relaxation  time  limits, 
in  the  very  least  to  encompass  virtually  instantaneous  responses.  However  this  is  likely  to  produce 
a  muddled  formulation,  with  limited  generality.  It  would  imply  a  continuous  distribution  of  relax¬ 
ation  times  between  the  “instantaneous”  one  and  all  the  others,  with  equalization  of  magnitudes. 
Basically,  these  “cures”  of  the  model’s  shortcomings  would  constitute  unphysical  tweakings  of 
the  parameters  or,  in  the  least,  an  unsupported  flight.  While  they  might  succeed  in  matching  an 
individual  FD  spectral  curve,  the  results  would  be  unlikely  to  generalize  successfully.  That  is, 
application  of  the  resultant  X  to  other  instruments  or  circumstances  would  likely  fail. 

CONCLUSIONS/  SIGNIFICANCE: 

1 .  Application  of  the  scaling/calibration/conversion  factors  so  we  can  compare  FD  and  TD  re¬ 
sponses  in  the  same  terms  provides  a  fundamental  insight  into  how  different  kinds  of  in¬ 
struments  respond  to  significantly  different  aspects  of  soil  magnetic  susceptibility.  For  the 
various  backyard  soils  examined  thus  far,  results  suggest  that  the  slow  relaxation  TD  re¬ 
sponse  is  objectively  much  smaller  than  the  FD  response,  which  is  dominated  by  effectively 
instantaneous  response  components. 

2.  The  calibration  factors  that  produced  the  evaluation  of  relative  magnitude  of  FD/TD  re¬ 
sponse,  as  in  Fig.  4.2.10,  can  easily  be  applied  further  to  calculate  the  factors  in  an  expanded 
soil  magnetic  susceptibility  model,  as  per  (4.2.3).  This  information  is  then  potentially  trans¬ 
ferable  to  other  EMI  instruments  of  whatever  type.  As  the  material  above  shows,  existing 
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models  are  not  transferable.  Heretofore,  FD  data  implied  no  TD  response;  TD  data  implied  a 
mueh  weaker,  slower  response  than  was  seen  in  the  FD.  The  expanded  model  eovers  all  these 
bases. 

3.  Item  #1  suggests  that  FD  instruments  will  be  at  a  disadvantage  for  diserimination  in  instanees 
where  signifieant  ground  responses  eannot  readily  be  subtraeted  as  baekground.  The  GEM- 
3D-I-  eounteraets  this  disadvantage  by  reeording  horizontal  field  eomponents  that  are  largely 
immune  to  ground  response  and  that  ean  help  distinguish  the  eontribution  of  target  vs  ground 
in  typieally  stronger  vertieal  eomponents. 

4.3  Estimating  magnetic  susceptibility  from  EMI  data 

4.3. a  Executive  Summary 

Studies  have  showed  that  magnetieally  suseeptible  soils  signifieantly  affeet  EMI  sensor  perfor- 
manee,  whieh  in  return  reduee  the  sensors  diserimination  eapabilities.  In  order  to  improve  EMI 
sensors  deteetion  and  diserimination  performanees  first  soil’s  magnetie  suseeptibility  needs  to  be 
estimated,  and  then  the  soils  EMI  responses  have  to  be  taken  into  aeeount  during  geophysieal  data 
inversion  proeedure.  Until  now  the  soil’s  magnetie  suseeptibility  is  determined  using  a  tiny  amount 
(up  to  15  mg)  of  soil’s  probe.  This  approaeh  in  many  eases  does  not  represent  effeetive  magnetie 
suseeptibility  that  affeets  on  the  EMI  sensors  performanees.  This  paper  presents  an  approaeh  for 
estimating  soil’s  magnetie  suseeptibility  from  low  frequeney  eleetromagnetie  induetion  data  and  it 
is  designed  namely  for  the  GeoPhex  frequeney  domain  GEM-3  sensor.  In  addition,  a  numerieal 
eode  ealled  the  method  auxiliary  sourees  (MAS)  is  employed  for  establishing  relation  between 
magnetieally  suseeptible  soil’s  surfaee  statisties  and  EMI  seattered  field.  Using  the  MAS  eode 
EMI  seatterings  are  studied  for  magnetieally  suseeptible  soils  with  two  types  of  surfaees:  body  of 
revolution  (BOR)  and  3D  rough  surfaee.  To  demonstrate  applieability  of  the  teehnique  first  the 
magnetie  suseeptibility  is  inverted  from  frequeney  domain  data  that  were  eolleeted  at  Cold  Re¬ 
gions  Researeh  and  Engineering  Eaboratory’s  test-stand  site.  Then,  several  numerieal  results  are 
presented  to  demonstrate  the  relation  between  surfaee  roughness  statistie  and  EMI  seattered  fields. 

4.3. b  Introduction 

Active  EMI  sensors  induce  magnetization  via  their  primary  fields.  The  static  magnetic  sensors 
are  sensitive  to  the  presence  of  magnetic  minerals  because  the  earth’s  magnetic  field  creates  an 
induced  magnetization  in  them.  This  anomalous  magnetic  field  produced  by  the  soil  magneti¬ 
zation  will  be  superimposed  on  the  magnetic  anomalies  produced  by  UXO’s  and  other  metallic 
objects  and  will  therefore  introduce  adverse  effects  into  the  discrimination  problem.  There  are 
three  magnetic  effects  that  impact  the  magnetic  and  electromagnetic  characteristics  of  the  subsur¬ 
face:  (4.3.1)  induced  magnetization,  (4.3.2)  viscous  remanent  magnetization  (VRM);  and  (4.3.3) 
permanent/remanent  magnetization 

The  magnetically  susceptible  and  heterogeneous  soils  are  a  major  problem  for  the  near  field 
EMI  sensing  technologies.  Studies  [42-51]  showed  that,  magnetically  susceptible  soils  can  pro¬ 
duce  electromagnetic  anomalies  of  the  same  magnitude  as  buried  metallic  targets.  Even  that  in 
those  studies,  it  is  assumed  that  the  spatial  distribution  of  magnetic  anomalies  are  constant,  similar 
to  a  half  space,  and  that  this  response  can  be  subtracted  from  measured  data,  the  practice  shows  that 
the  magnetic  soils  significantly  degrades  performances  and  applicability  of  the  EMI  technologies 
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for  UXO  discrimination.  Since,  for  discriminating  between  UXO  and  non-UXO  targets  aeeurately 
a  dense  measurement  grid  is  required.  Under  these  eonditions,  variations  in  sensor  height  and  ori¬ 
entation,  as  well  as  small  variations  in  the  surfaee  topography,  ean  produce  anomalies  similar  to 
those  from  UXO  [44].  However,  in  most  existing  approaehes  to  buried  metallie  objeet  diserimina- 
tion,  the  objeet  of  interest  is  assumed  to  be  embedded  in  a  medium  with  the  same  eleetromagnetie 
properties  as  free-spaee.  Any  influenee  of  the  baekground  medium  assumed  to  have  been  removed 
by  pre-proeessing  or  filtering  before  the  data  are  submitted  to  an  inversion  routine.  Until  now, 
teehniques  for  ealeulating  the  baekground  effeets,  whieh  inelude  the  half  layered-spaee  earth  and 
an  equivalent  dipole  layer  models  that  attempts  to  model  the  baekground  geology,  have  almost 
exelusively  negleeted  soils  surfaee  roughness  and  these  teehniques  are  limited  to  the  use  of  one 
dimensional  or  layered  Earth  models.  Sueh  models  are  ineapable  of  eapturing  the  true  proeess 
of  geomagnetie  induetion  that  is  assoeiated  with  the  surfaee  roughness,  as  it  is  the  ease  in  real 
field  environment.  Therefore,  one  must  eharaeterize  soil  response  so  that  it  ean  be  predieted,  or 
so  that  it  ean  at  least  be  aeeounted  for  during  data  processing.  The  characterizing  of  soil’s  EMI 
response  requires  knowledge  of  its  magnetie  suseeptibility.  To  date  the  soil’s  magnetie  suseepti- 
bility  is  measured  using  a  very  small  amount  (up  to  15  mg)  of  soil’s  probe,  whieh  in  many  eases 
does  not  represent  effeetive  magnetic  susceptibility  that  affeets  on  the  state  of  the  art  EMI  sensors 
performanees. 

This  seetion  investigates  EMI  seattering  from  magnetically  susceptible  random  rough  surfaees 
for  establishing  statistieal  relations  between  surfaee  and  signal  eharaeteristies.  The  study  is  done 
for  the  handheld  GEM-3  sensor.  Eirst,  to  understand  interaction  between  magnetieally  susceptible 
soil  and  the  sensor  here,  an  EMI  magneto  quasi-statie  (MQS)  problem  is  solved  is  solved  using 
a  numerieal  approaeh  ealled  the  method  of  auxiliary  sourees  (MAS)  [52-55].  In  MQS  magnetie 
fields  are  irrotational,  and  ean  thus  be  represented  effieiently  using  a  magnetie  sealar  potential. 
In  the  MAS,  boundary  value  problems  are  solved  numerieally  by  representing  the  eleetromagnetie 
fields  in  eaeh  domain  of  the  strueture  under  investigation  by  a  finite  linear  eombination  of  analytieal 
solutions  of  the  relevant  field  equations,  eorresponding  to  sourees  situated  at  some  distanee  away 
from  the  boundaries  of  eaeh  domain.  The  ’’auxiliary  sourees”  produeing  these  analytieal  solutions 
are  ehosen  to  be  elementary  dipoles/eharges  loeated  on  fietitious  auxiliary  surfaoe(s),  usually  eon- 
forming  to  but  offset  slightly  from  the  aetual  surfaee(s)  of  the  strueture.  Enforcement  of  standard 
eleetromagnetie  boundary  eonditions  at  an  array  of  points  over  the  objeet’s  aetual  surface  allows 
us  to  solve  for  the  auxiliary  sources,  from  whieh  we  ean  immediately  express  all  EM  fields  in  the 
problem. 

The  seetion  is  organized  as  follows:  In  seetion  II,  overviews  very  briefly  EMI  problem  and 
the  MAS  Galerkin  teehnique  for  MQS  regime.  Section  III  describes  an  approaeh  for  estimating 
soil’s  magnetie  suseeptibility  from  EMI  data,  seetion  IV  demonstrates  several  numerieal  results, 
and  finally  seetion  V,  presents  eonelusions  and  discusions. 

4.3.C  EMI  Scattering  for  Magentic  Rough  Soil 

To  illustrate  the  MAS  for  MQS  problem,  let  us  eonsider  a  permeable  3D  magnetieally  suseeptible 
rough  surfaee.  The  surfaee  is  illuminated  with  an  eleetromagnetie  field  produeed  by  a  GEM- 
3  sensor  (Eig.  4.3.1).  The  eomputational  region  is  divided  into  two  regions:  the  region  above 
the  surfaee  is  referred  to  as  Region  0  with  free-spaee  eleetromagnetie  parameters,  and  Region 
1  eorresponding  to  the  volume  under  the  rough  surface  with  ei=l  but  with  a  nonzero  magnetie 
suseeptibility.  As  is  well  established  in  the  MQS  regime,  the  displaeement  currents  are  negligible. 
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The  primary  field  induces  magnetic  dipoles  within  the  Region  1,  which  in  turn  produce  secondary 
(scattered)  fields.  The  magnetic  field  inside  each  region  can  be  expressed  as 


Hair)  =  -V'R«(r) 


(4.3.1) 


where 


1  f  Pads 

An^al  \r-ra 

Sa 

X>0 


(4.3.2) 


Here  'Ra(r)  is  the  fundamental  solution  of  the  Poison  equation,  a=0,l,  /i«  is  the  permeability 


GEM-3 


Figure  4.3.1:  3-D  rough  surface  used  for  computational  model. 

of  region  a,  and  pa  is  the  unknown  magnetic  charge  source  coefficient  distributed  on  the  a 
auxiliary  surfaces.  These  sources  reside  on  an  auxiliary  surface  slightly  outside  and  approximately 
conforming  to  the  real  surface.  We  compute  the  total  field  inside  each  a  =  0,1  region  by  assuming 
that  these  sources  radiate  into  an  infinite  homogeneous  medium  with  the  properties  of  the  corre¬ 
sponding  region. 

The  electromagnetic  fields  at  each  interface  between  two  regions  are  connected  by  boundary 
conditions.  Specifically,  the  magnetic  scalar  potential  and  the  normal  component  of  the  magnetic 
flux  density  must  be  continuous  across  the  real  surface: 

iBxW-t”,)]  =  0  (4.3.3) 

=  0  (4.3.4) 

where  h  is  the  unit  normal  vector  on  the  interface  between  the  two  regions.  For  achieving  high 
accuracy  and  saving  computational  resources,  the  boundary  value  problem  is  solved  using  the  MAS 
together  with  the  standard  Galerkin  technique.  To  demonstration  the  MAS-Galerkin  approach,  let 
us  rewrite  the  boundary  value  problem  in  the  following  compact  form 

Lp  =  V/  (4.3.5) 
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Here  L  is  an  integral-differential  operator,  p  is  auxiliary  magnetie  charge  density  and  t//  is  the  right 
hand  side  of  (4.3.2).  Let  a  weighing  function  w  and  (/,  w)  denote  the  inner  product  of  two  scalar 
functions  /and  w,  usually  defined  as 


{f  .w)  ^  j  fwds'  (4.3.6) 

here  ds'  is  a  surface  element.  Then,  by  using  Galerkin’s  method,  one  can  rewrite  equation  (4.3.6) 
as 

{Lp,w)  =  {xj/,w)  (4.3.7) 

upon  discretization  of  equation  (4.3.7)  a  linear  system  of  equations  for  unknown  p  coefficients 
results.  Once  the  amplitudes  of  the  p  coefficients  are  determined,  the  field  at  all  points  in  both 
regions  can  be  readily  computed. 


Figure  4.3.2:  Experimental  setup  of  GEM-3  sensor  at  various  heights  above  the  ground. 


4.3.d  Estimation  of  Soil’s  Magnetic  Susceptibility  from  GEM-3  Measurement  Data 

The  GEM-3  sensor  consists  of  two  transmitter  loops.  Based  on  the  standard  expression  for  the  field 
from  a  wire  loop  in  the  magneto-quasistatic  regime,  the  z  component  of  the  magnetic  field  radiated 
by  the  GEM-3  sensor  along  its  z  rotational  axis  can  be  expressed  as: 


Hz  = 


(af+z2) 


3/2 


+ 


:  +  z2) 


3/2 


(4.3.8) 


According  to  image  theory,  for  a  GEM-3  sensor  placed  at  an  elevation  z  above  a  permeable  half 
space,  the  magnetic  field  produced  at  z  by  the  ground  can  be  written  as 


Hz  = 


l^gr  1  ^0^  1 


1 


1 


l  +  Pgr  2  \(^aj  +  4z^y^  (af+16z2) 


,3/2 


(4.3.9) 
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where  /j  =  Io,andIg  =  has  been  assumed.  So,  if  we  know  the  /„  eurrent  and  the  a\  radius 
then  we  ean  ealeulate  the  permeability  of  ground  jigr.  However,  different  EMI  sensors  reeord 
data  in  different  units.  For  example,  the  GEM  sensor  produees  field  in  “ppm”  units  relative  to  a 
separate  referenee  reeeiver  eoil.  Given  the  known  geometry  and  assumed  relation  of  the  eurrents 
in  the  GEM,  the  aetual  seattering  field  from  any  target  (metal  or  ground)  with  o  and  fi,  ean  be 
written  in  the  following  form 

H  =  lunF  (o,  II, geometry)  =  Hdata/p  (4.3.10) 


lun  is  the  unknown  eurrent  and  p  is  some  unknown  sealing  faetor  (aeeounting  for  amplifieation, 
filtering,  the  size  of  the  referenee  eoil,  ete.)  that  produees  the  GEM’s  readings  in  terms  of  “ppm”. 
Onee  the  sealing  faetor  is  known,  then  the  soil’s  magnetie  suseeptibility  ean  be  estimated  as 


Z  =  Mgr-l 


2y 

1-7 


(4.3.11) 


here 


Fldata  Q  _  ^0^1  I  1  1 

^  \  (aj -f  4^2) (aj -b  16^2)^/^ 


(4.3.12) 


Using  methods  explained  in  See.  4.4  the  sealing  faetor  p  was  determined  in  [56]  for  the  GEM- 
3  sensor  with  =  IQA  eurrent  to  he  p  =  4.32  ■  10^.  This  value  was  determined  from  detailed 
measured  data  for  a  sphere. 


Frequency  [l-|z 


Figure  4.3.3:  Mean  value  of  the  predicted  magnetic  susceptibility  versus  frequency. 
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Figure  4.3.4:  1-D  Rough  surface.  In  meters. 


Figure  4.3.5:  Scattered  magnetic  field  for  1-D  rough  surface. 
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4.3. e  Results 

Experimental  data  were  eolleeted  using  GEM-3  sensor  at  CRREE  faeilities.  To  estimate  soil’s 
magnetie  suseeptibility  the  measurement  were  taken  over  a  partieular  spot  as  a  funetion  of  sensor 
elevation.  Eor  this  the  “blue  rig”  (see  Eig.  4.3.2  left)  was  used.  Eor  the  test  18  shots  were  taken  as 
follows:  three  on  the  top  of  the  rig,  as  in  the  pieture  above  left;  9  deseending  from  the  top  to  within 
10  em  of  the  base;  then  the  rig  was  removed  and  3  shots  were  taken  at  sueeessive  holes  about  10 
em  apart  on  a  board  plaeed  on  the  ground  where  the  rig  had  been  (Eig.  4.3.2  right);  then  the  blue 
rig  was  put  baek  in  position  and  3  shots  were  again  taken  on  the  top.  The  mean  value  (respeet  to 
the  evaluation)  of  predieted  soil 

suseeptibility  real  and  imaginary  parts  for  CRREE  soil  versus  frequeney  is  shown  on  Eig.  4.3.3. 
The  results  show  that  the  soil  has  a  eomplex  magnetie  suseeptibility.  The  amplitude  of  the  real  part 
is  around  0.38  x  10^^,  and  it  is  deereasing  as  frequeney  inereases,  as  expeeted.  The  imaginary  part 
also  shows  very  small  but  frequeney  dependent  trends,  indieating  that  the  present  soil  will  affeet  on 
the  time  domain  EMI  sensors  performanees  as  well.  We  intend  to  deploy  the  time  domain  EM-63 
sensor  on  the  same  soil  in  near  future,  and  outeome  investigations  will  be  reported  in  subsequent 
publieations. 

Onee  the  soil’s  magnetie  suseeptibility  was  determined  then  the  MAS  eode  was  used  for  GEM- 
3  sensor  for  establishing  statistieal  relation  between  surfaee  roughness  and  EMI  seattered  field 
for  magnetie  soil.  Eirst  BOR  rough  surfaees  were  generated  using  Id  eode  [57]  for  generating 
Id  BOR  rough  surfaee.  A  snap  shot  of  the  surfaee  in  a  eylindrieal  eoordinate  system  is  depieted 
on  Eig.  4.3.4.  On  this  figure  vertieal  and  horizontal  axis  eoineide  with  the  z  and  p  axis  in  the 
eylindrieal  eoordinate  system,  respeetively.  ID  line  on  Eig.  4.3.4  was  rotated  around  the  vertieal 
axis  for  generating  the  BOR  surfaee.  The  GEM-3  sensor  that  is  the  exeitation  souree  was  plaeed 
at  15  em  distanee  from  the  origin  (0,0  point  on  Eig.  4.3.4). 

Eor  eaeh  fixed  eorrelation  length  and  RMS  heights  100  realizations  were  simulated.  One  of 
sueh  runs  shown  on  Eig.  4.3.5.  The  standard  deviation  and  mean  are  shown  on  Eig.  4.3.6.  These 
results  show  that  as  eorrelation  length  and  RMS  height  inereases  the  EMI  response  from  magnetie 
soil  inereases  as  expeeted.  The  results  also  show  that  for  high  eorrelation  length  (surfaee  is  beeom- 
ing  more  flat)  and  for  high  RMS  height  (sensor  is  elose  to  the  surfaee),  the  standard  deviation  of 
the  EMI  response  deereases  slightly. 

Einally,  we  studied  EMI  seattering  for  3D  rough  surfaees.  Eirst,  we  eonstruet  3d  rough  surfaees. 
One  of  rough  surfaees  is  shown  on  Eig.  4.3.6.  The  GEM-3  sensor  was  plaeed  25  em  above  the 
surfaee.  Using  the  MAS-Galerkin  teehnique  the  full  EMI  seattering  problem  was  solved.  Eor  eaeh 
fixed  eorrelation  length  and  RMC  height  EMI  seatted  filed  was  eomputed  for  100  realizations. 
Eigure  4.3.8  shows  the  seattered  magnetie  field  versus  realization  number.  The  standard  deviation 
and  mean  values  are  depieted  on  Eig.  4.3.6  for  3d  rough  surfaee,  with  fixed  eorrelation  length.  The 
results  show  that  both  standard  deviation  and  mean  values  inerease  as  RMS  height  of  the  3d  rough 
surfaees  inereases. 

4.3. f  Conclusion 

In  this  seetion,  the  EMI  seeondary  response  from  a  magnetieally  suseeptible  soil  with  rough  sur¬ 
faees  were  investigated  and  analyzed.  The  eombined  MAS  -Galerkin  teehnique  is  employed  for 
investigating  EMI  seattering  from  3d-rough  surfaees.  The  relation  between  surfaee  statisties  and 
EMI  seattered  field  is  investigated  and  presented.  In  future  work  the  relation  between  sensors  el¬ 
evation,  orientation  and  surfaee  roughness  statistie  will  be  studied  and  reported  in  a  subsequent 
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Correlation  length  [cm]  ®  ^  RMS  height  [cm] 


Figure  4.3.6:  Standard  deviation  and  mean  value  for  the  BOR  rough  surfaces  versus  RMS  height  and  corre¬ 
lation  length. 
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Figure  4.3.7:  3-D  random  rough  surface. 


Figure  4.3.8:  Scattered  EMI  field  from  3d-rough  surface. 
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Figure  4.3.9:  Standard  deviation  and  mean  value  for  the  3d  rough  surfaces  versus  RMS  height  at  fixed 
surface  correlation  length. 
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4.4  Absolute  Calibration  of  EMI  Measurements  and  Application  to  Soil  Magnetic  Suscep¬ 
tibility  Inference 

This  section  presents  a  method  for  determining  a  conversion  factor  for  electromagnetic  induction 
(EMI)  survey  measurements  so  that  fields  received  in  an  instrument’s  native  units  can  be  translated 
into  standard  (“absolute”)  EM  units.  The  method  is  non-invasive,  requiring  only  a  rudimentary 
characterization  of  the  instrument,  such  as  loop  geometry  and  the  basic  nature  of  the  output  such 
as  the  magnetic  field,  differencing  of  magnetic  field,  or  time  derivative  of  magnetic  field.  It  does 
not  require  knowledge  of  the  device’s  internal  electrical  engineering  features,  such  as  amplifica¬ 
tion,  filtering,  or  transfer  functions  between  components.  The  technique  proceeds  by  comparing 
data  from  controlled  measurements  to  model  results  and  thus  allows  one  to  do  similar  comparisons 
henceforth.  As  an  example  application  and  reasonableness  check,  we  use  a  particular  frequency 
domain  (ED)  instrument  to  infer  soil  magnetic  susceptibility  value  in-situ  directly  from  survey  data. 
The  same  methodology  is  applied  to  a  new  time  domain  (TD)  instrument.  This  allows  validation 
of  the  basic  scaling  methodology  via  benchmark  cases  and  also  illustrates  its  transferability.  Con¬ 
version  of  data  from  both  instrument  types  into  standard  units,  for  measurements  from  the  same 
soil,  illuminates  the  contrasting  nature  of  the  soil  responses  to  TD  and  ED  sensors.  It  also  points 
to  the  fundamentally  different  magnitudes  of  different  components  of  soil  magnetic  susceptibility 
in  the  soil  studied,  relating  to  instantaneous  response  and  relaxation  response. 

We  acquired  data  over  canonical  targets  and  over  some  different  soil  types  when  we  first  ob¬ 
tained  the  GEM-3^”''.  Using  these  data  we  were  able  to  extract  scaling  factors  for  the  three  receiver 
coils  on  the  GEM-S^"*".  We  were  also  able  to  estimate  the  susceptibility  of  the  soils  we  measured 
directly  from  the  measured  data. 

Data  Collected 

•  Data  taken  over  Spheres  and  other  canonical  targets 

•  Data  collected  over  2  different  soil  types 

4.4.a  Introduction  and  Background 

To  understand  the  needed  for  an  ability  to  convert  a  sensor’s  output  into  standardized  units,  it 
is  helpful  to  view  the  UXO  problem  from  an  inversion  perspective.  Inversion  studies  use  the 
secondary  field  to  infer  the  target’s  properties.  Possibly  the  most  basic  example  is  the  estimation 
of  soil  properties,  particularly  broadband  magnetic  susceptibility.  EMI  responses  caused  by  this 
property  of  soil  can  be  a  major  problem  in  surveys  seeking  to  discriminate  UXO  [45,  48,  50,  58, 
59].  The  ability  to  estimate  susceptibility,  in  situ,  directly  from  sensor  data  would  provide  valuable 
soil  characteristic  data  in  an  extremely  efficient  manner.  A  commonly  used  method  of  inversion 
is  to  match  the  measured  secondary  response  with  an  analytical  model’s  response  of  homogenous 
or  layered  half-space  media.  The  susceptibility  of  the  soil  will  correspond  to  the  values  of  the 
model’s  susceptibility  parameters  that  produce  the  best  match.  In  many  prior  studies  of  this  nature, 
the  secondary  fields  are  reported  in  units  of  parts  per  million  (ppm)  of  the  primary  field  intensity 
[60-65]. 

Examining  the  past  work  reveals  two  key  aspects  which  differentiate  our  work.  Eirst,  all  in¬ 
struments  in  these  prior  studies  use  receiver  and  transmitter  coils  consisting  of  laterally  separated 
horizontal  current  loops.  To  our  knowledge,  no  similar  study  has  been  made  with  more  compact, 
complex  instruments  where  the  transmitting  and  receiving  coils  may  be  coaxial.  Secondly,  more 
significantly,  the  prior  studies  assumed  that  the  results  reported  by  their  instruments  were  truly 


Benjamin  Barrowes,  Kevin  O’Neill 


-66- 


r  AlSSvl  THAYER  SCHOOl  OF 


ENGINEERING 
AT  DARTMOUTH 


MM- 15  37,  GEM-3  ^+-  Final  Report 


4  MATERIALS  AND  METHODS:  GEM-3^+  SOFTWARE 
4.4  Absolute  Calibration  of  EMI  Measurements 


ppm  relative  to  a  speeifie  primary  field  sample.  However,  modern  EMI  deteetors  inelude  filters, 
amplifiers,  ete,  attaehed  to  partieular  eomponents,  the  details  of  whieh  a  general  user  has  no  way 
of  knowing.  Therefore  the  user  has  no  real  way  of  knowing  the  relation  between  aetual  transmitted 
and  reeeived  fields  in  A/m  or  Teslas  and  the  aetual  units  of  measurements. 

This  situation  argues  for  a  simple,  fast,  non-invasive  method  of  both  assessing  the  performanee 
and  ealibrating  the  measurements  of  instruments  that  do  not  report  data  in  known,  universal  units 
to  produee  a  mateh  with  the  seeondary  fields  predieted  by  a  model.  The  latter  is  required  to 
produee  seientifie  inferenees,  sueh  as  estimates  of  ground  magnetie  suseeptibility  using  the  method 
deseribed  earlier.  Thus,  the  ealibration  essentially  amounts  to  translating  reported  ppm  data  to  be 
eonsistent  with  the  inversion  researeher’s  model. 

The  need  for  this  translation  faetor  within  UXO  inversion  studies  ean  be  seen  in  the  researeh 
to  diseriminate  UXO  from  elutter.  Maehine  learning  teehniques  ean  make  use  of  the  modeled 
responses  of  a  wide  range  of  metallie  objeets  to  train  a  systematie  diseriminator  that  separates 
UXO  from  elutter  objeets  on  the  basis  of  their  EMI  responses  [66,  67].  To  elassify  objeets  with 
measured  data,  previous  studies  frequently  normalize  all  measured  and  simulated  data  [67]  to 
ensure  eonsisteney.  This  normalization  entails  a  loss  of  valuable  information  about  the  magnitude 
of  the  response.  Eor  example,  a  very  small  pieee  of  elutter  may  likely  produee  a  weaker  signal 
than  a  large  UXO  if  both  are  equidistant  from  the  sensor.  Normalization  would  erase  sueh  basie 
distinetion.  However,  if  one  were  able  to  eonvert  data  into  units  eonsistent  with  the  standardized 
units  of  the  modeled  data — in  other  words,  to  ensure  that  the  modeled  responses  and  measured 
responses  are  direetly  eomparable — no  information  would  be  lost.  Therefore,  the  reeovery  of  this 
eonversion  relationship  is  a  prerequisite  for  many  types  of  advaneed  inversion  studies. 

4.4.b  Test  Instrumentation 

As  a  handy  example  of  an  ultra-wideband  frequeney  domain  (ED)  EMI  instrument,  here  we  use 
the  GEM-3  sensor,  manufaetured  by  Geophex  [4,  64]  and  used  in  inversion  researeh  [68,  69].  Our 
ehoiee  of  the  GEM-3  is  largely  due  to  availability  and  relative  ease  of  use.  As  mentioned  earlier, 
this  seetion  illustrates  a  general  method  of  ealibration  that  may  be  applied  to  many  EMI  instru¬ 
ments.  This  analysis  of  the  GEM-3  instrument  ean  be  direetly  applied  to  any  non-adaptive  EMI 
instrument  whieh  reports  measurements  proportional  the  transmitted  primary  field  or  whieh  uti¬ 
lizes  a  fixed  but  possibly  unknown  transmitted  field  strength.  Eollowing  a  detailed  treatment  of  the 
method  in  applieation  to  the  GEM-3,  we  then  apply  it  to  a  new  time  domain  (TD)  instrument,  the 
Man-portable  Veetor  Time  Domain  Eleetromagnetie  Sensor  (MPV).  This  deviee  has  just  reeently 
been  developed  in  an  ongoing  researeh  projeet  [70]. 

To  determine,  through  non-invasive  means,  a  eonversion  faetor  that  ean  translate  our  partieular 
GEM-3  instrument’s  reeord  of  the  seeondary  field  into  standard  units,  we  eompare  the  aetual  re¬ 
eeived  fields  from  a  oanonieal  objeet  to  the  fields  determined  by  a  model.  Thereby,  this  analysis 
permits  one  to  aequire  the  ability  to  evaluate  the  sensor’s  data  in  standardized  quantities  in  all  other 
future  applieations.  Eurthermore,  we  examine  the  signifieanee  of  the  diserete  size  of  the  GEM-3 
reeeiver  loop  in  the  isolation  of  the  eonversion  faetor.  Ultimately,  with  proper  ealibration,  data  in 
the  GEM-3 ’s  native  units  from  magnetieally  responsive  soil  is  used  direetly  to  estimate  broadband 
magnetie  suseeptibility  of  the  soil. 

As  an  additional,  eomparative  test  we  utilized  a  new  time  domain  (TD)  instrument  designed  and 
manufaetured  under  the  direetion  of  one  of  the  authors  [70].  This  time  domain  instrument  eonsists 
of  two  eoaxial  eoils  of  37.5om  radius  that  generate  a  step  down  funetion.  Erom  46  mieroseeonds 
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F  H  F 

Tx  Tx->REF 


Figure  4.4. 1 :  Flow  chart  of  instrument’s  process,  to  be  paralleled  by  the  modeling  sequence.  If  normalization 
is  performed  by  the  instrument  then  the  lower,  shaded  set  of  boxes  may  be  deleted  in  the  model  and  an 
arbitrary  (e.g.,  unit)  current  assumed. 


to  about  20  milliseconds  after  shut  off,  the  step  down  response  is  then  measured  at  a  receiver  loop 
located  at  the  center  of  the  transmitter  coils. 


4.4.C  Recovery  of  the  Conversion  Factor 

A  given  transmitter  (Tx)  current  in  any  instrument  will  produce  a  particular  magnetic  field  at  a 
target’s  location  (e.g.  in  A/m),  which  in  turn  will  produce  a  particular  secondary  field  value  at  the 
instrument’s  receiver  (Rx).  In  the  case  of  the  GEM-3,  the  Rx  signal  is  normalized,  frequency  by 
frequency,  by  a  sample  of  the  primary  field  in  a  secondary  Rx  coil.  Thus,  in  effect,  the  instrument 
provides  a  signal  at  each  frequency  that  is  proportional  to  received  field  per  unit  Tx  current,  so  we 
need  not  know  the  specific  Tx  current  to  model  what  occurs.  We  do  require,  however,  knowledge 
of  the  geometry  of  the  transmitter  and  receiver  loops. 

Through  our  analysis,  we  can  confirm  if  the  proportionality  between  the  actual  received  field 
and  the  sensor  measurement  signal  is  truly  fixed.  While  the  signal  may  be  selectively  amplified 
or  filtered,  this  manipulation  is  always  done  in  the  same  way,  without  adaptive  measures  which 
would  alter  that  relation.  Therefore,  the  recovered  conversion  factor  should  have  no  dependence 
on  measurement  position  or  object  type. 

We  can  readily  calculate  the  response  of  standard  objects  to  fields  produced  by  arbitrary  as¬ 
sumed  fixed  Tx  currents  in  the  known  Rx  geometry,  and  thus  determine  the  proportionality  be¬ 
tween  measured  and  actual  fields.  The  standard  objects  consist  here  of  three  spheres  made  of 
aluminum,  brass,  and  steel.  The  underlying  physical  process  by  which  the  instrument  produces 
signals  is  presented  symbolically  in  the  flow  chart  of  Fig.  4.4.1  and  equations  that  follow. 

Production  of  the  instrument’s  signal,  in  the  case  of  the  GEM-3,  is 


G{ppm) 


hx  ■  Ftx  ■  Htx^s  •  S  ■  Hs^rx  ■  Frjc 
hx  ■  Ftx  ■  Htx-*ref  ■  Free 


(4.4.1) 


where  hx  is  the  current  in  the  transmitter  loops.  Fjx  is  what  the  instrument  does  with  this  current 
in  producing  the  Tx  field.  In  the  case  of  the  GEM-3  this  contains  little  besides  the  geometry  of 
the  coils,  as  we  may  simply  assume  that  the  current  is  continuously  on  during  the  measurement 
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portion  of  the  cycle.  Htx->s  is  a  transfer  function  yielding  the  field  produced  by  the  transmitter  at 
the  target  location.  S  is  the  inherent  target  response  to  a  unit  input,  while  Hs->rx  is  the  function 
that  converts  that  response  to  the  field  at  the  Rx  location.  Frx  is  the  process  that  is  performed  on 
that  field  by  the  receiver,  including  any  associated  amplification,  filtering,  etc.  The  quantities  in 
the  denominator  are  explicable  in  similar  terms,  except  that  REF  refers  to  the  reference  coil,  in  a 
fixed  position  relative  to  to  the  Tx  coils,  where  the  Tx  field  is  sampled.  We  do  not  need  to  know 
any  of  its  particulars. 

The  process  producing  the  actual  (modeled)  field  M  at  the  Rx  location  can  be  represented  as 


M{A/m)  =  Im-FTx-HTx^s-S-Hs^Rx-F^, 


(4.4.2) 


where  Im  is  the  current  assumed  in  the  transmitter.  may  include  integration  of  the  secondary 
field  over  the  area  of  the  Rx  coil  but,  unlike  Frx  for  the  actual  instrument,  need  not  contain  any  other 
functions  performed  by  the  sensor’s  electronics  or  processors.  Note  that  both  G  and  M  contain  the 
same  product  Fjx  ■  Htx-*s  ■  S  ■  Hr^rx  ,  which  we  consider  known  (calculable).  The  ratio  of  recorded 
data  to  actual  (calculated)  field  values  is  thus 


(4.4.3) 


The  product  /„  ■  is  known;  we  may  also  normalize  by  taking  /„  to  be  unity  and  adjusting  M 
accordingly.  The  quantity  in  parenthesis  is  a  function  of  the  internals  of  the  instrument  and  is  what 
we  infer  as  a  lumped  factor  from  measurements  of  the  ratio  on  the  left.  When  this  is  done,  the 
conversion  factor  has  been  obtained. 

This  method  to  recover  the  conversion  factor  can  also  be  used  for  any  comparable  instrument 
by  accounting  for  its  particular  loop  geometry.  Note  that,  to  be  fully  comparable,  an  instrument 
must  have  the  signal  normalized  by  Tx  field  or  Tx  current  or  use  a  known,  fixed  Tx  field  or  current. 
Also,  the  GEM-3  data  is  processed  within  the  instrument  so  that  it  is  proportional  to  the  Rx  field 
value,  not  to  the  time  derivative  thereof.  Simple  alterations  of  the  method  can  be  brought  to  bear 
to  account  for  differences  in  any  of  these  particulars.  For  example,  in  our  TD  MPV  instrument,  the 
current  is  shut  off  prior  to  measurement;  the  data  is  not  inherently  normalized  by  the  Tx  current  or 
field;  and  the  output  is  proportional  to  the  time  derivative  of  B  and  not  to  B.  This  simply  means 
that  we  must  heed  these  facts  when  comparing  the  calculated  and  measured  responses  from  our 
standard  objects,  making  sure  that  the  processes  parallel  one  another  by  adjusting,  respectively,  the 
content  of  Frx,  the  denominator  in  equation  (4.4.3) ,  and  the  contents  of  Frx  andF^. 

4.4.c.(l)  Modeling  of  Metallic  Spheres  For  the  frequency  domain,  we  model  metallic  spheres  as 
a  special  case  of  spheroids  in  the  magnetoquasistic  regime.  This  analytical  solution  is  well  vali¬ 
dated  and  produces  EMI  secondary  field  predictions  for  spheroidal  objects  with  any  conductivity 
and  permeability  values  and  in  response  to  arbitrary  excitation  [71-73].  The  solution  provides 
secondary  field  predictions  in  terms  of  A/m.  The  values  can  be  obtained  at  any  observations  point, 
e.g.  at  the  nominal  location  of  an  Rx  coil  or  over  arrays  of  quadrature  points  within  its  area.  The 
three  specific  spheres  that  were  modeled  were  all  3  inches  in  diameter.  In  addition,  we  imple¬ 
mented  the  time  domain  step  down  response  of  a  metallic  sphere  to  a  uniform  field  as  given  by 
Wait  and  Spies,  1969[38].  Given  the  relatively  small  size  of  the  sphere  and  the  distance  between 
the  sphere  and  MPV  transmitter  loops  in  all  measurements,  the  incident  field  is  approximately 
uniform.  Therefore,  Wait  and  Spies ’s  solution  is  a  valid  model. 
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4.4.c.(2)  Modeling  of  the  Primary  field  The  excitation  field  is  generated  by  an  analytical  model 
that  represents  the  GEM-3  instrument  as  a  set  of  idealized  wire  loops.  The  transmitting  sensor 
head  consists  of  two  concentric  current  loops  in  series  with  radii  of  20  cm  and  11.074  cm.  The 
smaller  loop  is  placed  inside  the  larger  loop  and  serves  as  a  bucking  coil,  approximately  nulling  the 
primary  field  at  the  center  of  the  sensor  head,  where  the  receiver  coil  resides.  There  are  8  windings 
in  the  outer  and  4  windings  in  the  inner  transmitter  coil.  We  assume  in  the  calculations  that  the 
wires  in  the  transmitter  coils  carry  10  Amps,  which  is  the  case  at  most  frequencies.  However,  as 
mentioned  earlier,  this  varies  somewhat  by  frequency  but  because  of  the  normalization  by  Tx  field, 
the  results  are  not  affected. 

We  use  the  complete  elliptical  integral  functions  of  the  first  and  second  kind  to  characterize  the 
magnetic  vector  potential  and  vector  field  from  a  circular  current  loop  in  a  cylindrical  coordinate 
system.  This  calculation  method  is  well  established  [19,  74,  75].  For  a  current  loop  of  radius  a  and 
normal  to  the  z  axis  in  the  z=  Zo  plane  and  with  /„  current  flowing  in  the  ^girection,  the  transverse 
and  axial  fields  are. 


where  . 


(4.4.6) 


and  jiQ  is  the  permeability  of  the  surrounding  medium.  The  standard  cylindrical  coordinates  of 
radius,  rotational  angle  around  the  z  axis,  and  elevation  are  expressed  by  (r,  z).  K  and  E  are  the 

complete  elliptical  integral  functions  of  the  first  and  second  kind,  respectively.  These  equations  will 
produce  the  magnetic  field  from  a  current  loop  at  any  observation  position.  The  total  primary  field 
is  created  from  the  superposition  of  two  current  loops  that  correspond  to  the  two  actual  transmitting 
loops  in  the  sensor  head.  Our  implementation  of  the  model  agrees  with  the  alternative  method  of 
using  numerical  integration  and  application  of  Biot-Savart’s  Law  but  has  the  added  benefit  of  being 
more  computationally  efficient. 

For  the  time  domain  MPV  instrument,  the  analytical  TD  solution  for  a  sphere’s  response  to 
a  uniform  step  down  function  only  requires  calculation  of  the  equivalent  of  a  static  uniform  field 
generated  by  a  current  loop  and  observed  at  the  location  of  the  sphere.  This  can  be  done  through 
the  current  loop  equations  provided  above  where  the  radius  of  the  transmitter  loop  now  37.5cm 
and  the  number  of  windings  is  100  to  correspond  to  the  MPV. 

4.4.c.(3)  Modeling  of  the  Receiver  Loop  Given  the  finite  size  of  the  receiver  loop,  a  circular  loop 
with  a  radius  of  6cm  for  the  GEM-3,  the  secondary  field  will  not  necessarily  be  uniform  across 
the  area  of  this  loop.  This  is  especially  true  if  the  object  is  close  to  the  instrument.  Therefore 
all  modeled  secondary  fields  are  numerically  integrated  over  the  loop  area  to  parallel  the  actual 
secondary  field  measurement.  This  integration  is  calculated  by  partitioning  the  loop  area  into 
concentric  rings.  Then  each  ring  is  again  segmented  along  its  length  such  that  the  entire  area  is 
partitioned  into  curved  trapezoids.  In  our  experience,  over  300  divisions  are  needed  for  consistent 
results. 
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Figure  4.4.2:  Block  diagram  illustrating  the  search  algorithm. 


4.4. c.(4)  Measurements  Measurements  were  taken  above  eaeh  metal  sphere  by  plaeing  the  GEM- 
3  instrument  at  points  in  a  vertieal  grid.  This  was  aeeomplished  by  measuring  along  horizontal 
lines  above  one  another  at  different  elevations.  The  horizontal  point  spaeing  was  in  2  em  inere- 
ments  from  -30  em  to  -i-30  em  from  the  position  above  the  sphere.  Vertieal  spaeing  of  the  lines  was 
in  5em  inerements  from  13  em  to  28  em  above  the  surfaee  of  the  sphere.  The  frequeney  range  was 
from  30  Hz  to  47  kHz.  This  wide  range  of  measurement  points  both  in  spaee  and  frequeney  was 
taken  to  ensure  that  the  recovered  scaling  factor  is  consistent:  there  must  not  be  any  variation  due 
to  object  distance,  frequency,  or  choice  of  target  object.  In  other  words,  all  the  data  must  point  to 
a  scaling  factor  of  the  same  value. 

4.4. c.(5)  Matching  Algorithm  We  performed  a  Fibonacci  search  [7 6]  within  the  range  of  possible 
values  for  the  scaling  factor  and  the  sphere’s  material  properties  to  find  the  best  match  between  the 
model’s  prediction  and  the  data.  A  Fibonacci  search  first  partitions  the  search  range  into  two 
halves.  The  median  value  within  each  half  is  given  to  the  model  as  input  to  produce  two  different 
outputs.  Whichever  half  produces  the  more  accurate  output  will  become  the  new  search  space  and 
the  process  repeats.  The  overview  of  this  algorithm  is  illustrated  in  Fig.  4.4.2.  Accuracy  is  gauged 
by  utilizing  the  normalized  mean  square  error  as  an  objective  function.  This  type  of  simple  search 
is  relatively  fast  but  usually  prone  to  converging  on  a  solution  that  is  not  the  best  overall  match. 
In  our  investigation  that  difficulty  was  not  encountered.  The  matching  algorithm  is  independently 
repeated  for  the  frequency  response  at  each  measurement  point  for  each  sphere.  This  independence 
allows  us  to  examine  the  results  for  consistency  over  space  and  frequency. 

The  main  difficulty  in  the  search  algorithm  was  that  the  metallic  spheres’  exact  permeability 
(/i)  and  conductivity  (a)  values  were  not  precisely  known.  These  spheres  were  purchased  off  the 
shelf  from  a  third  party  and  not  manufactured  to  any  specifications  other  than  that  they  were  of  the 
correct  size  and  type  of  metal.  Therefore,  to  match  the  model’s  EMI  prediction  with  measurements, 
both  the  scaling  factor  and  /i  and  a  must  be  simultaneously  recovered.  As  a  simplification,  the 
aluminum  and  brass  spheres  can  be  assumed  to  be  non-permeable,  and  thus  only  their  o  values 
and  the  scaling  factor  need  to  be  recovered.  Approximate  values  of  sigma  for  these  materials  are 
available  from  standard  textbooks,  furnishing  a  reasonableness  check.  For  the  highly  permeable 
steel  sphere,  only  the  ratio  between  a  and  /i  affects  the  secondary  response  [77].  Therefore, 
relative  /i  is  assumed  to  be  100  and  the  normalized  o  value  is  recovered  along  with  the  scaling 
factor.  Adopting  a  higher  relative  /i  value  produces  no  change  in  the  resulting  conversion  factor. 

Fig.  4.4.3  shows  the  match  between  the  modeled  response  and  the  measurement  of  the  steel 
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Data  Fit  for  Steel 


Figure  4.4.3:  Figure  showing  the  matching  between  model  and  measurement  of  a  steel  sphere  at  a  single 
grid  point. 


sphere.  The  match  does  extremely  well  save  for  the  lowest  and  highest  frequency  points.  We  have 
observed  previously  that  the  GEM-3  can  have  difficulties  capturing  consistent  responses  at  very 
low  and  very  high  EMI  frequencies.  Therefore,  the  lowest  and  highest  frequency  data  are  ignored 
during  the  matching  process.  Eurthermore,  we  proceeded  to  recover  the  scaling  factor  of  the  MPV 
instrument  in  the  same  fashion  and  with  similarly  good  results.  The  matching  results  for  only  the 
GEM-3  will  be  presented  to  avoid  redundancy. 

4.4.d  Results 

After  performing  the  Eibonacci  search  for  the  material  parameters  of  the  spheres  and  finding  the 
best  scaling  factor  for  each  GEM-3  measurement  point,  we  then  averaged  the  scaling  factors  for 
all  objects  recovered  at  the  closest  9  measurement  points.  We  arrived  at  an  overall  value  of  4.93  x 
10^.  Division  by  this  scalar  will  convert  all  GEM-3  measurements  into  the  average  magnetic  field 
that  would  pass  through  the  receiver  coil  within  our  model,  given  our  specific  transmitter  loop. 

All  the  data  can  be  converted  into  a  percent  error  from  the  averaged  conversion  factor  of  4.93  x 
10^.  This  variation  of  the  scaling  factor  can  likewise  be  plotted  over  space.  Shown  in  Eig.  4.4.4  is 
the  variation  of  the  scaling  over  space  for  the  steel  sphere.  As  the  figure  demonstrates,  the  variation 
is  minimal  with  most  of  the  error  under  5%  error.  We  see  similarly  good  results  for  the  brass  and 
aluminum  spheres  shown  in  Eigs.  4.4.5  and  4.4.6,  respectively.  Since  the  model  of  the  GEM-3 
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Figure  4.4.4:  Plot  of  the  variation  of  the  retrieved  conversion  factor,  as  a  percent  error  from  4.933  x  10^,  for 
the  steel  sphere  as  a  function  of  position. 


used  a  fixed  eurrent  in  the  transmitter  eoils,  these  results  verify  that  the  ppm  reported  by  GEM-3  is 
eonsistently  normalized  by  some  sample  of  the  primary  transmitted  field.  In  addition  to  providing 
a  eonversion  faetor,  the  method  also  potentially  provides  a  quiek  way  for  inversion  researehers  to 
aseertain  if  their  instruments  are  funetioning  eorreetly  with  respeet  to  any  normalization  produeed 
internally  in  the  instrument. 

If  we  ean  assume  the  seeondary  field  is  relatively  eonstant  over  the  reeeiver  loop,  then  it  would 
be  unneeessary  to  numerieally  integrate  the  fields  within  the  loop  area.  However,  eomparison 
of  the  eonsisteney  of  results  with  and  without  sueh  integration  shows  that  signifieant  errors  are 
introdueed  by  failure  to  integrate  over  the  reeeivers. 


4.4.e  Application  of  the  Conversion  Factor 

As  an  example  applieation  and  reasonableness  eheek,  we  apply  the  eonversion  faetor  to  the  prob¬ 
lem  of  obtaining  an  estimate  of  in  situ  soil  magnetie  suseeptibility.  The  solid  eurves  in  Fig.  4.4.8 
show  the  Inphase  and  Quadrature  response  of  a  partieular  sample  of  soil  queried  with  the  GEM-3 
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Figure  4.4.5:  Plot  of  the  variation  of  the  retrieved  conversion  factor,  as  a  percent  error  from  4.933  x  10^,  for 
the  brass  sphere  as  a  function  of  position. 


at  12.5  cm  above  the  surface.  The  Inphase  curve  shows  virtually  no  ehange  over  frequency,  and 
the  Quadrature  component  is  negligible.  Therefore,  we  implemented  a  model  of  the  homogeneous, 
permeable,  non-condueting  half-space  response  to  the  GEM-3  excitation  by  the  image  method,  as 
per  Wait,  1985  [78].  The  output  of  this  model  was  eonverted  from  model  units  to  GEM-3  ppm  units 
through  multiplieation  of  the  recovered  eonversion  factor.  Then  a  d.e.  soil  suseeptibility  value,  K  = 
70  X  10  ^  SI,  was  found  to  provide  the  best  match  to  data  as  shown  in  Eig.  4.4.8.  This  susceptibility 
value  is  solidly  within  the  range  of  reasonable  soil  values  as  reported  by  others  [33-35]. 

The  example  above  provides  primarily  an  order  of  magnitude  eheek,  whieh  is  often  all  we  may 
require  in  terms  of  knowledge  of  soil  suseeptibility.  It  also  provides  an  illustration  of  ED  soil 
response,  whieh  we  ean  now  eompare  to  TD  response.  By  reeovering  the  eonversion  faetor  for 
multiple  instruments,  one  ean  also  make  meaningful  eomparison  between  data  from  those  instru¬ 
ments,  e.g.  for  soil  response,  illuminating  the  different  aspeets  of  the  suseeptibility  that  dominate 
the  responses  of  the  respeetive  sensor  types.  To  this  end,  the  eonversion  faetor  of  6810  for  the 
TD  MPV  instrument  was  obtained  in  exaetly  the  same  manner  as  for  the  GEM-3.  Results  of  this 
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Figure  4.4.6:  Plot  of  the  variation  of  the  retrieved  conversion  factor,  as  a  percent  error  from  4.933  x  10^,  for 
the  aluminum  sphere  as  a  function  of  position. 


comparison  are  shown  in  Sec.  4.2. 

4.4.f  Conclusion  and  Discussion 

The  formulation  presented  here  shows  how,  by  reeourse  to  simple  modeling,  one  may  obtain  faetors 
to  convert  EMI  instrument  output  to  received  field  values  in  standard  EM  units.  As  all  instruments, 
even  of  the  same  type,  generally  have  at  least  slightly  different  ealibration,  the  method  is  best  ap¬ 
plied  to  each  individual  sensor  to  be  used.  In  any  case,  the  formulation  makes  clear  the  adjustments 
needed  to  treat  different  instrument  types.  Tests  reported  here  reveal  that  the  finite  size  of  reeeiver 
loops  may  have  to  be  taken  into  account  to  obtain  a  precise  scaling  factor,  even  when  those  loops 
are  small. 

An  appeal  of  the  method  is  that  it  does  not  require  any  detailed  knowledge  of  the  internal 
electronics  and  processing  of  the  instrument  under  consideration.  As  long  as  basie  geometry  of  the 
transmitter  coils  is  known  together  with  a  few  other  rudimentary  instrument  characteristics,  one 
can  obtain  actual  received  field  values  from  data  by  one-time  non-invasive  ealeulation. 
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Figure  4.4.7:  Plot  of  the  variation  of  the  retrieved  conversion  factor,  as  a  percent  error  from  4.243  x  10^,  for 
the  aluminum  sphere  as  a  function  of  position.  No  integration  was  performed  over  the  receiver  loop. 


As  an  example  applieation  and  reasonableness  test,  we  used  our  retrieved  sealing  faetor  to 
mateh  the  modeled  response  of  permeable,  suseeptible  soil  to  measurements  of  soil  and  have  found 
reasonable  soil  suseeptibility  values.  This  not  only  supports  the  validity  of  our  sealing  faetor  but 
also  points  to  the  possibility  of  eharaeterizing  soil  from  sensor  measurements  alone,  avoiding  the 
need  to  eonduet  laboratory  analysis  on  eore  samples  outside  of  in  situ  eonditions.  By  this  means, 
one  sensor’s  data  ean  yield  geophysieal  information  direetly  that  is  applieable  to  other  instruments. 
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Figure  4.4.8:  Comparison  between  GEM-3  measured  inphase  (I)  and  quadrature  (Q)  ground  response  and 
prediction  using  half-space  soil  model  with  an  applied  scaling  factor. 


4.5  Generalized  Standardized  Excitations  Approach 

The  generalized  standardized  exeitation  approach  (GSEA)  is  presented  to  enhance  UXO  discrim¬ 
ination  under  realistic  field  conditions.  The  GSEA  is  a  fast,  numerical,  forward  model  for  rep¬ 
resenting  an  object’s  EMI  responses  over  the  entire  frequency  band  from  near  DC  to  100s  of 
kHz,  including  all  physical  relevant  physical  considerations,  such  as  near  and  far  field  phenomena, 
arbitrary  material  and  geometrical  heterogeneity,  and  internal  target  interactions.  It  has  been  de¬ 
veloped  and  tested  in  both  the  frequency  and  time  domains  for  actual  UXOs  placed  in  free  space. 
The  GSEA,  which  uses  magnetic  dipoles  instead  of  magnetic  charges  as  responding  sources,  is 
capable  of  taking  into  account  the  background  medium  surrounding  an  object.  Given  a  modeled 
UWB  frequency  domain  (ED)  response,  the  corresponding  time  domain  (TD)  response  is  easily 
obtained  by  the  inverse  Eourier  transform.  Thus  the  technique  is  applicable  for  any  ED  or  TD 
sensor  configuration  and  can  treat  complex  data  sets:  novel  waveforms,  multi-axis,  vector,  or  ten¬ 
sor  magnetic  or  electromagnetic  induction  data,  or  any  combination  of  magnetic  and  EMI  data. 
Host  media  effects  are  taken  into  account  via  appropriate  types  of  Green’s  function  and  equivalent 
dipole  sources.  Comparisons  between  simulations  and  experimental  data  illustrate  that  the  GSEA 
is  a  unified  approach  for  reproducing  both  TD  and  ED  EMI  signals  for  actual  UXOs.  The  EMI 
response  from  a  soil  that  has  a  frequency-dependent  magnetic  susceptibility  is  studied.  The  EMI 
responses  in  both  ED  and  TD  domains  are  analyzed  for  the  model  of  an  actual  UXO  that  is  buried 
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in  a  magnetically  susceptible  half  space. 

4.5.a  Introduction 

Cleanup  of  unexploded  ordnance  (UXO)  is  greatly  complicated  by  the  difficulties  in  distinguishing 
UXO  from  non-UXO  metallic  debris,  found  at  most  UXO  sites.  The  problem  becomes  much  more 
complicated  when  signals  are  contaminated  by  noise  that  originates  from  magnetically  susceptible 
and  electrically  conductive  soils  [42-48,  50,  51,  79-81].  Until  now,  in  most  existing  approaches 
to  UXO  classification,  the  object  of  interest  is  assumed  to  be  placed  in  a  free  space  [55,  68,  82- 
90].  Any  influence  of  the  host  medium  is  considered  to  be  removed  by  filtering  before  data  are 
submitted  to  an  inversion  algorithm.  Recent  studies  show  that  “geologically  hostile”  sites  cause 
significant  problems  for  magnetometers  and  electromagnetic  induction  (EMI)  sensors,  in  terms 
of  both  decreased  probability  of  detection  and  increased  probability  of  false  alarm.  In  regions 
of  highly  magnetic  soil,  magnetometry  and  electromagnetic  sensors  often  detect  large  anomalies 
that  are  of  geologic,  rather  than  of  metallic,  origin.  For  instance,  [79]  documented  the  problems 
encountered  at  the  former  Naval  Training  Range  on  Kaho’olawe  Island  when  using  EMI  sensors. 
During  production  surveys  at  the  site,  approximately  30%  of  identified  anomalies  were  from  false 
positives  due  to  geology,  attributed  to  the  strong  magnetic  viscosity  exhibited  by  the  basaltic  soils. 

There  is  also  a  need  for  detection  and  discrimination  of  UXO  in  undersea  environments  [91]. 
In  this  case,  unlike  land  surveys  where  the  conductivity  (<  10^^[S/m])  of  non-permeable  soils  can 
be  neglected,  it  is  impossible  to  neglect  EMI  responses  due  to  the  conductivity  of  seawater  and  the 
ocean  bottom  in  marine  surveys.  All  of  this  leads  to  a  high  level  of  false  alarms,  which  translates 
into  an  increased  workload  because  each  detected  anomaly  must  be  treated  as  if  it  were  an  actual 
UXO.  Therefore,  innovative  discrimination  techniques  that  apply  to  any  field  condition  and  that 
reliably,  quickly,  and  accurately  distinguish  between  hazardous  UXO  and  non-hazardous  metallic 
items  are  required.  To  address  this  issue  here  we  present  a  generalized  standardized  excitation 
approach  (GSEA)  that  is  suitable  for  complex  data  sets:  novel  waveforms,  multi-axis,  vector, 
tensor  magnetic  or  electromagnetic  induction  data,  or  any  combination  of  magnetic  and  EMI  data. 

The  SEA  for  objects  placed  in  free  space  is  described  in  detail  in  [55,  88,  90].  This  work 
extends  the  SEA  for  more  general  cases  by  using  magnetic  dipoles  as  responding  sources  instead 
of  magnetic  charges  as  in  [55,  88,  90].  By  using  magnetic  dipoles  with  the  corresponding  dyadic 
Green’s  function  the  GSEA  becomes  applicable  for  objects  placed  in  a  conducting  as  well  as  a 
permeable  host  medium. 

The  GSEA  can  be  briefly  outlined  as  follows.  For  any  given  object  the  amplitudes  of  the  re¬ 
sponding  magnetic  dipoles  are  determined  and  sorted  in  the  universal  library  for  any  number  of 
basic  spheroidal  modes.  Then,  any  primary  field  is  decomposed  into  a  set  of  basis  excitations, 
which  are  then  multiplied  by  appropriate  weights  (e.g.,  spheroidal  modal  decomposition  coeffi¬ 
cients  from  which  we  can  calculate  the  target’s  complete  response  just  by  superposing  responses 
of  each  basic  excitation).  The  key  element  in  the  GSEA  is  to  determine  the  amplitudes  of  the  re¬ 
sponding  sources,  which  are  characteristic  only  of  the  object,  and  are  independent  of  sensor  type, 
object  location  and  orientation,  and  transmitted  waveforms. 

There  are  two  ways  to  determine  the  amplitudes  on  the  responding  magnetic  dipoles:  (1)  using 
measurement  data  and  (2)  solving  the  full  3-D  EMI  problem  in  detail.  The  most  straightforward 
way  to  determine  the  amplitudes  of  the  responding  magnetic  dipole  sources  is  to  solve  a  standard 
inverse  problem  based  on  the  measured  data.  Obviously,  this  process  requires  very  good  experi¬ 
mental  conditions  and  a  sufficient  number  of  independent  measurements  of  an  object  of  interest  in 
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Figure  4.5.1:  Problem  geometry  and  reduced  set  of  sources  distributed  along  rings  on  the  spheroid 
surface. 


order  to  reduce  the  degree  of  ill-posedness.  The  ill-posedness  makes  the  solution  inaccurate.  Re¬ 
cently,  [88]  applied  such  a  data-derived  approach  to  extracting  the  modal  response  coefficients  for 
each  candidate  by  carefully  designing  the  measurements  at  different  distances  and  orientations  in 
free  space.  However,  the  accuracy  and  reliability  of  the  model  parameters  determined  in  this  way 
may  not  always  be  satisfactory  due  to  unavoidable  measurement  noise  and  numerical  difficulties 
arising  from  the  inherent  ill-conditioning  of  the  problem,  although  a  special  treatment  was  applied 
in  [88].  In  addition,  the  model  parameters  were  obtained  from  measurement  data  with  the  given 
sensor.  Currently  all  available  EMI  sensors  have  certain  limitations  in  both  frequency  band  and 
time.  Therefore,  the  amplitudes  of  the  responding  source  that  are  derived  from  these  measured 
data  have  limitations,  they  can’t  cover  all  possible  EMI  and  magnetic  data,  and  they  could  not  be 
used  to  obtain  EMI  response  in  time  domain  for  different  wave  forms  and  sensors. 

Originally,  in  [55]  and  [90],  to  determine  the  amplitudes  of  responding  sources,  a  numerical 
procedure  based  on  the  method  auxiliary  source  (MAS)  and  hybrid  MAS  thin  skin  approximation 
(MAS/TSA)  was  proposed  and  it  is  generalized  here.  In  this  procedure,  by  utilizing  a  full  3-D  EMI 
solver,  the  modal  responding  coefficients,  or  strengths  of  a  reduced  source  set  (RSS)  [55,  90]  are 
determined  by  employing  a  physically  complete  numerical  simulation  of  the  object’s  response  to 
each  fundamental  excitation  mode.  The  full  MAS  model-based  approach  has  an  advantage  over 
the  data-derived-based  approach,  because  it  is  a  well-posed  EMI  problem  and  it  is  not  dependent 
on  measurements.  In  this  approach  the  RSS  can  be  obtained  very  accurately  for  any  excitation 
mode  and  in  an  ultra-wide  band  frequency  range.  Thus,  the  technique  allows  users  to  calculate 
the  EM  field  in  both  the  frequency  and  time  domain  and  to  control  the  number  of  input  spheroidal 
modes.  Here  the  TD  EMI  response  for  a  given  UXO  is  calculated  directly  from  ED  RSS  sources 
just  by  using  the  convolution  theorem  [92-94] . 
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The  section  is  organized  as  follows:  In  Sec.  4.5.b,  the  generalized  standardized  excitation  ap¬ 
proach  is  presented,  Sec.  4.5.c  describes  the  frequency-dependent  magnetic  susceptibility  model, 
and  Sec.  4.5.d  shows  several  experimental  and  numerical  results,  demonstrating  the  applicability 
of  the  GSEA  in  both  frequency  and  time  domain,  and  for  frequency-dependent  magnetic  soil. 


4.5.b  Generalized  Standardized  Excitation  Approach  (GSEA) 

Recently,  the  SEA  has  been  developed  and  applied  to  UXO  discrimination  [55,  88,  90].  All  those 
studies  assume  that  an  object  is  placed  in  free  space.  Here,  the  SEA  is  generalized  to  take  into 
account  conducting  and  magnetically  susceptible  host  media  effects.  To  illustrate  the  GSEA,  let  us 
assume  that  an  object  is  placed  in  a  background  with  magnetic  permeability  /ii  and  conductivity 
Oi,  Eig.  4.5.1.  The  object  is  illuminated  by  an  arbitrarily  oriented,  time-varying  primary  magnetic 
field.  We  surround  the  object  with  a  fictitious  spheroid,  which  is  introduced  only  as  a  computational 
aide  in  the  decomposition  of  the  primary  magnetic  field  into  fundamental  spheroidal  modes.  We 
choose  spheroids  because  they  can  assume  the  general  proportions  of  elongated  objects  of  interest, 
such  as  UXO,  which  are  also  typically  bodies  of  revolution  (BOR).  Oblate  spheroids  can  also 
be  used  for  flattened  shapes.  In  general,  the  fictitious  surface  could  be  a  smooth  closed  surface, 
as  applicable  for  a  related  standardized  source  set  approximation  described  in  [55,  90].  On  the 
fictitious  spheroid  given  hy  ^  =^o  (Eig.  4.5.1),  the  primary  magnetic  field  can  be  expressed  as: 

OO  OO  1 

ff'”'  =  E  E  E  (4.5.1) 

m=0n=m p=0 

The  bpmn  are  coefficients  needed  to  express  the  primary  field  and  is  the  pmn  mode  of  the 
primary  magnetic  field  component  when  bpmn  =1  •  The  normal  component  of  the  primary  magnetic 
field  on  the  fictitious  spheroid,  can  be  written  as 

f-f  P  ^  N  1 

Hi' —  E  E  E  (4,5,2) 

m=0n=m p=0 

where  are  the  standard  spheroidal  coordinates;  d  is  the  spheroid’s  interfocal  distance,  E™ 

are  associated  Eegendre  functions  of  the  first  kind,  and  Tpm  (0)  represents  the  azimuthal  functions 
given  by 

=  I  P  =  “  (4,5,3) 

pniKYj  y  sm(m<^),  p=l. 

By  the  orthogonality  of  the  associated  Eegendre  functions,  the  spheroidal  expansion  coefficients 
bpmn  can  be  derived  as 

1  iTt 

7“  =  «»■  <P)TrMd<pdn,  (4,5,4) 

where  7=2form  =  p  =  0  and  7=1  otherwise.  The  integration  in  (4.5.4)  is  evaluated  by  numerical 
integration.  This  completes  the  decomposition  of  the  primary  field  (t],  0)- 

After  the  primary  magnetic  field  is  decomposed  into  the  pmn  spheroidal  modes,  then  the  com¬ 
plete  solution  for  the  target  to  each  field  is  obtained.  Since  the  object  is  placed  in  a  conduct¬ 
ing  and  magnetically  susceptible  background,  the  magnetic  field  in  the  entire  computational  space 
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(Regions  1  and  2,  Fig.  4.5.1)  satisfies  Helmholtz’s  wave  equation  and  it  ean  be  represented  with 
magnetic  dipoles  as 

liaij)  =  I  Bairy)-Paif^)ds\  (4.5.5) 

s 

Here  Ga  is  the  Dyadic  Green’s  Function,  Pa{P)  are  amplitudes  of  magnetic  dipoles,  and  a 
=  1,2  corresponds  to  region  1  or  2.  Here 

=  /=  VV\  p-jkaR 

Ga  (r;r)  =  U  +  7^  {f;r)  ,Ga  {n?')  =  - - -  (4.5.6) 

ka  =  \/(0^allo{(0£o  +  jOa) 

where  7  is  the  unit  dyad,  ka  is  the  wave-number  in  a  region,  Oa  and  /daare  the  conductivity 
and  relative  permeability  of  a  region  respectively,  and  R=  |r  — r'l  is  the  distance  between  source 
and  observation  points.  We  assume  that  the  relative  electric  permittivities  of  both  regions  are 
1;  Ga{f;f')  is  the  fundamental  solution  for  the  wave  equation  of  the  Hertzian  magnetic  vector 
potential,  whereas  the  exp  { jot}  time  convention  has  been  implied  and  suppressed. 

On  the  surface  of  the  object,  total  magnetic  fields  satisfy  the  following  boundary  conditions: 

+HP^)]  =  [fix  7/2]  (4.5.7) 

fi-Ati(//^"  +  //n  =  n-PiHi  (4.5.8) 

where  fiis  a  unit  normal  vector  on  the  real  surface,  H  is  the  primary  magnetic  field,  H2  is  the  total 
magnetic  field  in  region  2,  and  /i«  is  a  region’s  relative  magnetic  permeability. 

To  determine  amplitudes  of  the  magnetic  dipoles  we  have  to  solve  the  entire  boundary  value 
EMI  problem  for  each  mode  of  the  primary  field  just  once  and  then  store  the  amplitudes  of  the 
responding  sources.  To  do  so,  one  extends  each  mode  to  the  physical  surface  within  the  enclosing 
spheroid  and  applies  the  3-D  MAS-TSA  method  [55,  90].  Finally,  the  target’s  response  for  each 
primary  magnetic  field  component  bpmnHpmn  is  expressed  similar  to  (4.5.5)  as: 

_  ^  =  _ 

=  (>p™£G,(r,fi)P!'7  (4.5.9) 

k=l 

where  Fi^™”is  the  amplitude  of  the  k^^  auxiliary  magnetic  dipole,  located  at  the  target’s  domain 
[24-25],  corresponding  to  the  //p„„response,  and  N  is  the  number  of  auxiliary  magnetic  sources. 
Using  P\  ^  for  each  fundamental  mode,  the  total  response  at  any  point  outside  the  scatterer  can  be 
represented  as: 

00  CX3  1  TV  

=  E  E  E  v»E5i(r,77r- 

m=0n=mp=0  k=\ 

Thus,  after  pre-computation  of  the  coefficients  for  any  given  object,  the  EMI  scattering 
problem,  for  any  particular  3-D  configuration  involving  it,  breaks  down  merely  to  determine  the 
spheroidal  modal  expansion  coefficients  bp,nn- 

In  equation  (4.5.10)  a  substantial  number  of  responding  sources  P\  ^  are  required  to  represent 
the  scattered  magnetic  field  outside  the  object  (including  physical  surface).  Similar  to  (4.5.10)  we 
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can  re-express  that  field  quite  aeeurately  in  terms  of  a  reduced  number  of  sourees  from  the  fietitious 
spheroid.  The  amplitudes  of  this  redueed  number  of  sourees  for  eaeh  input  pmn  spheroidal 

mode  ean  be  determined  by  solving  a  linear  system  of  equations  for  the  normal  eomponent  of  the 
seattered  magnetie  field  as  it  is  shown  in  [55,  90],  and  finally  the  eomplete  seeondary  magnetie 
field  ean  be  represented  as 


H'Uf) 


oo  oo  1  ^red _ 


m=0n=m  p=0  ;=1 


(4.5.11) 


Note  that,  while  we  may  ultimately  be  able  to  express  the  seattered  field  using  a  small  number 
of  sourees,  this  is  a  fundamentally  different  strategy  from  what  has  been  applied  heretofore  in 
the  simple  independent  dipole  models.  In  the  latter,  eaeh  souree  responds  only  to  the  primary 
field  striking  it  loeally.  However,  here  the  pf"*”,  i  =  1,2, . .  .,Nfed  responding  sourees  aet  together, 
not  in  response  to  loeal  stimuli  but  to  express  the  response  of  the  entire  objeet  to  the  distributed 
exeitation  of  the  pmn  mode.  In  addition,  using  the  MAS-MAS/TSA  numerieal  eode  the  ^1’"“  can 
be  generated  and  stored  for  any  number  of  pmn  spheroidal  modes  in  the  infinite  series  (4.5.11). 
Onee  this  is  done,  then  the  truneation  eriterion  ean  be  determined  from  the  input  primary  magnetie 
field  easily. 

Overall,  the  entire  GSEA  approaeh  ean  be  deseribed  briefly  as  follows: 

1.  For  a  given  UXO  amplitudes  of  responding  magnetie  dipoles  pf'”",  i  =  l,2,...A^^gj  rings 
are  determined  and  sorted  in  the  universal  library  for  any  number  of  basie  spheroidal  mode 
exeitations  pmn  =  1,2,.... 

2.  Onee  step  1  is  done,  then  for  a  given  sensor  the  primary  field  is  deeomposed  into  spheroidal 
modes,  the  spheroidal  modal  deeomposition  bpmn  ooeffleients  are  ealeulated,  and  the  neees- 
sary  number  of  spheroidal  modes  is  determined. 

3.  Use  redueed  set  of  sourees  (RSS)  to  ealeulate  EMI  response  for  eaeh  pmn-th  basie 
exeitation.  Seale  eaeh  pmn-EMl  response  on  the  bpmn  eoeffieients  and  ealeulate  the  target’s 
eomplete  response  by  just  superposing  responses  of  eaeh  basie  exeitation. 

The  GSEA,  whieh  is  based  on  the  MAS-TSA  and  introdueed  here,  can  produce  a  target’s  ultra- 
wideband  frequeney  response.  Thus,  the  proposed  GSEA  ean  be  used  direetly  to  obtain  an  objeet’s 
TD  EMI  responses  via  eonvolution  theorem  without  reealeulating  amplitudes  of  the  responding 
sourees.  This  makes  the  GSEA  a  unified  model  to  treat  both  ED  and  TD  data,  and  is  attraetive  from 
apraetieal  point  of  view,  as  many  state-of-the-art  EMI  sensors  (EM-63,  EM-61,  Zonge  NanoTEM) 
are  operating  in  TD. 

Eet  us  briefly  deseribe  the  important  formulas  required  to  eompute  the  TD  EMI  response  for 
a  general  eurrent  waveform  I{t)  flowing  in  a  transmitter  loop.  By  using  the  eonvolution  theorem 
[29]  the  indueed  voltage  in  the  reeeiver  eoil  ean  be  expressed  as 


A'{t-x)l'{x)  dx-A{0)l' 


(4.5.12) 


where  A  (t)  represents  an  objeet’s  impulse  response  and  the  prime  means  the  derivative  with  respeet 
to  time  t.  Equation  (4.5.12)  represents  the  TD  response  of  an  objeet  to  a  general  exeitation  eurrent 
/(t)souree. 
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4.5.C  Frequency-Dependent  Magnetic  Susceptibility 

There  are  three  types  of  magnetic  susceptibility  that  generate  soils’  EM  responses:  (1)  induced, 
(2)  permanent/remanent,  and  (3)  viscous  remanent  (VRM).  The  soil  magnetic  properties  are  deter¬ 
mined  by  the  presence  of  iron  and  iron-oxide  particles.  Permanent  remanent  magnetization  is  the 
magnetization  that  exists  in  the  absence  of  any  applied  field.  Induced  magnetization  is  the  magneti¬ 
zation  that  arises  in  the  presence  of  an  external  magnetic  field,  and  viscous  remanent  magnetization 
(VRM)  is  a  phenomenon  that  occurs  when  magnetization  of  an  object  placed  in  an  external  mag¬ 
netic  field  changes  in  a  time  relative  to  the  applied  field.  This  means  that  the  object’s  susceptibility 
is  a  complex  frequency-dependent  [45] 

Z(®)  =/(«)+ (4.5.13) 


where  mgs  the  angular  frequency,  j  is  the  unit  complex  number,  and  x'{0})  and  the  real 

and  imaginary  parts  for  the  frequency-dependent  magnetic  susceptibility.  There  are  various  refer¬ 
ences  [31  and  references  therein]  that  discuss  many  aspects  of  soil  magnetic  properties.  The  most 
common  frequency-dependent  complex  magnetic  susceptibility,  assuming  that  magnetic  relaxation 
time  constants  are  uniformly  distributed  between  times  X\  and  T2,  is  modeled  as  follows  [49]: 


Z(«) 


Xo{  1 


1 


■In 


j(OX2  + 1 


ln(T2/Tl)  j(OTl  +  l 


(4.5.14) 


where  ;t^o  is  the  DC  value  of  the  susceptibility. 

For  magnetically  susceptible  ground  in  the  present  of  a  metallic  object,  the  field  that  is  mea¬ 
sured  by  the  sensors  contains  two  parts  and  can  be  written  as 

F'”“(m)  =  W{co)  +h"^\co)  (4.5.15) 

where  and  are  respectively  magnetic  fields  produced  by  the  magnetically  sus¬ 

ceptible  soil  and  the  object.  The  magnetic  field  contains  all  interactions  between  the  object 
and  the  susceptible  host  medium.  For  determining  contribution  of  each  magnetic  field 

and  H"^\co)  in  the  measured  field,  a  full  EMI  problem  must  be  solved.  In  this  section  the  soil  is 
considered  to  be  a  uniform  half-space,  and  interaction  between  the  soil  and  the  object  is  taken  into 
account  by  using  the  image  source  method  [49]. 


4.5. d  Numerical  Results 

This  section  presents  some  numerical  and  experimental  data  that  demonstrate  the  GSEA  as  a  uni¬ 
fied  model  for  any  FD  or  TD  sensor  configuration,  and  its  applicability  to  complex  data  sets:  novel 
waveforms,  multi-axis,  vector,  or  tensor  magnetic,  or  electromagnetic  induction  data,  or  any  com¬ 
bination  of  magnetic  and  EMI  data;  and  ability  to  take  into  account  the  influence  of  conductive  and 
magnetically  susceptible  geological  soils  on  metal  detectors. 

4.5. d.(l)  The  universal  RSS  As  it  was  discussed  above,  the  reduced  source  set  (RSS)  in  Eq. 
(4.5.11)  depends  only  on  the  target’s  geometry  and  electromagnetic  properties.  To  validate  such  a 
unique  characteristic  of  the  RSS,  here  comparisons  between  RSS  modeled  and  experimental  data 
are  given  for  an  actual  UXO  (81  mm)  in  both  FD  and  TD.  The  data  were  collected  by  two  EMI 
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systems  that  are  widely  used  in  the  UXO  diserimination  eommunity:  (1)  a  wideband  frequeney- 
domain  sensor  (GEMS)  developed  by  Geophex  Ltd.  [95]  and  (2)  a  time-domain  instrument  (EM63) 
developed  by  Geonies  [96].  The  data  were  eolleeted  on  the  U.S.  Army  Engineer  Researeh  and 
Development  Center  test  stand  site. 


Figure  4.5.2:  Frequency-domain  EMI  response  for  81  mm  UXO;  (a)  GEM-3  excitation;  (b)  nose  up  vertical, 
(c)  nose  up  45  inclination,  (d)  horizontal. 


Figure  4.5.2  shows  the  eomparisons  between  GSEA  and  aetual  data  in  the  frequeney  domain. 
In  this  case  data  were  collected  for  the  UXO  oriented  in  three  different  directions  relative  to  the 
GEM-3  sensor’s  head:  (1)  vertical  tail  up,  (2)  45®  degree  nose  up,  and  (3)  transverse.  The  GEM-3 
frequency  range  is  from  30  Hz  to  50  kHz.  The  comparisons  between  measured  and  actual  exper¬ 
imental  data  are  in  very  good  agreement  for  all  orientations.  Note  that  the  RSS  produces  results 
in  an  ultra-wideband  frequency  range,  from  magneto- static  (0  Hz)  to  EMI  frequency  limit.  This 
allows  users  to  accurately  compute  the  scattered  field  at  any  required  frequencies  by  simple  inter¬ 
polation  and  to  obtain  EMI  responses  readily  in  TD  via  inverse  Fourier  transforms.  To  illustrate 
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this  capability,  here  the  first  impulse  responses  in  TD  are  evaluated  by  applying  the  digital  filter 
teehnique  [92]  to  the  inverse  sine  transform  [93]  as  follows: 

2 

A{t)  =  —  I  ImB{G))sm(Ot  dG)  (4.5.16) 

7t  Jo 

where  ImB{(o)  represents  an  imaginary  part  of  the  magnetie  flux  that  is  ealeulated  via  RSS.  Then 
the  indueed  voltage  is  ealeulated  using  the  time  eonvolution  teehnique  (4.5.12). 


Time  [m  s] 


Figure  4.5.3:  Time-domain  EMI  response  for  81 -mm  UXO;  (a)  EM-3  excitation;  (b)  nose  up  vertical,  (c) 
nose  up  45  inclination,  (d)  horizontal. 

In  the  EM63  instrument,  the  eurrent  waveform  eonsists  of  an  exponential  eurrent  inerease  fol¬ 
lowed  by  a  linear  ramp  off.  The  eurrent  has  the  three  pulses  per  measurement.  For  eomparisons 
between  RSS  and  TD  data,  here  the  same  81-mm  UXO  is  ehosen.  The  objeet  was  exeited  in  three 
orientations  1 :  (1)  vertieal  tail  up,  (2)  45®  degree  nose  up,  and  (3)  transverse  orientations.  For  all 
three  exeitations,  the  vertieal  distanee  between  the  sensor’s  transmitter  loop  and  the  eenter  of  the 
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cylinder  is  h  =  60.00  cm.  The  TD  indueed  voltage  is  calculated  by  an  inverse  Fourier  transform  of 
the  frequency-domain  magnetic  field  flux  using  equations  (4.5.12)  and  (4.5.16).  The  comparisons 
are  depieted  in  Fig.  4.5.3.  They  show  very  good  convergence  between  the  measured  and  modeled 
TD  data.  Therefore,  the  universal  RSS  should  be  able  to  produce  EMI  responses  for  a  given  target 
for  any  sensor  in  FD  and  TD  domains. 


Figure  4.5.4:  EMI  response  for  a  magnetically  susceptible  soil  versus  the  sensor  height  and  frequency,  (a  ) 
Inphase  part;  (b)  quadrature  part. 


4.5.d.(2)  EMI  response  for  soil  with  frequency-dependent  susceptibility  In  this  section,  first  FD 
EMI  response  is  studied  for  a  magnetically  susceptible  half-space.  The  half-space  is  illuminated 
by  a  ED  sensor.  In  these  simulations  for  the  sensor  model  the  following  parameters  are  used: 
current  /o  =1  A,  and  a  100-cm  x  100-cm  transmitter  loop.  The  soil’s  frequency-dependent  sus¬ 
ceptibility  is  assumed  to  be  the  same  as  in  equation  (4.5.14)  with  a  realistic  Xo  =  0.005  D.C  value 
of  susceptibility  [33],  and  X\  =  10^^  [sec]  and  X2  =  10^^  [sec]  time  constants.  Eigure  4.5.4  shows 
soil’s  responses  in-phase  (right)  and  quadrature  (left)  parts  as  a  function  of  frequency  and  sensor 
height.  The  results  illustrate  that  the  soil’s  EMI  responses  strongly  depend  on  both  the  frequency 
and  sensor  height.  As  the  sensor  approaches  the  soil,  soil’s  response  increases  and  it  stays  almost 
constant  for  antenna  heights  between  1  cm  and  10  cm.  Note  that  all  parameters  in  Eigure  4  are 
in  Eogarithmical  scale.  The  soil  response’s  quadrature  part  approaches  to  maximum  between  10 
kHz  and  100  kHz,  whereas  at  low  frequencies  the  in-phase  part  is  dominant.  At  highest  frequencies 
(more  than  lOOkHz)  both  parts  of  soil’s  response  decrease.  Thus  the  soil  with  frequency-dependent 
magnetic  susceptibility  produces  significant  EMI  responses  almost  entire  UWB  frequency  range. 

Einally,  to  illustrate  the  soil  effect  on  a  buried  object’s  EMI  responses,  several  numerical  exper¬ 
iments  were  done  in  both  ED  and  TD.  The  81 -mm  UXO  was  buried  under  a  magnetically  suscep¬ 
tible  half-space.  The  half-space  electromagnetic  parameters  are  exactly  the  same  {Xq  =  0.005,  Xi 
=  10^^  [sec]  and  X2  =  10^^  [sec])  as  in  the  previous  paragraph.  The  entire  structure  is  illuminated 
by  (1)  the  GEM-3  sensor  with  current  Ig  =1  A,  10-  and  20-cm  radii  coils,  and  (2)  the  EM-63  with 
/o  =1  A,  100-cm  X  100-cm  transmitter  loop.  The  sensors  are  placed  10  cm  above  the  half-space. 
The  UXO  is  oriented  45®  nose  up  and  its  center  is  at  32-cm  depth  for  the  GEM-3  excitation  and 
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50  cm  for  the  EM-63  sensor.  The  results  are  depicted  in  Fig.  4.5.5.  These  results  suggest  that 
a  magnetieally  suseeptible  half-spaee  could  significantly  affect  both  frequency  and  time  domain 
EMI  responses. 


Figure  4.5.5:  EMI  response  for  81  mm  UXO;  (a)  Frequency  domain,  (b)  Time  domain. 

These  ealeulations  are  based  on  the  funetional  form  for  ground  suseeptibility  given  by  4.5.14, 
whieh  the  Steve  Billings,  et  al.  at  UBC/Sky  state  has  “applied  to  every  soil  they’ve  ever  eneoun- 
tered”  in  the  TD.  Also,  the  value  of  5  x  10-3  is  on  the  order  of  magnitude  of  the  single  frequeney 
suseeptibility  values  obtained  by  Janet  Simms  in  numerous  measurements  at  the  UXO  standardized 
test  sites  [33-35].  However  the  question  arises  whether  this  is  the  appropriate  value  to  insert  for 
Xo-  In  partieular,  we  have  never  seen  ED  soil  effeets  of  the  magnitude  and  type  shown  in  Fig.  4.5.5, 
left.  This  may  mean  that  funetional  form  4.5.14  pertains  sueeessfully  to  TD,  but  that  another  larger 
effeet  appears  in  the  frequeney  domain  eorresponding  to  the  instantaneous  soil  response.  The  latter 
would  not  appear  (or  at  least  would  only  appear  indireetly)  in  the  TD  measurements.  The  method 
used  here  to  perform  “absolute”  ealibrations  of  instrument  data  for  eonversion  into  aetual  magnetie 
field  units  allows  us  to  infer  aetual  broadband  x  values  from  sueh  data.  Applying  this  to  both  the 
new  GEM-3^'^  as  well  as  emerging  TD  instruments  should  allow  us  to  eomplete  the  eharaeteriza- 
tion  of  soils  of  interest,  modifying  the  funetional  form  4.5.14  as  neeessary  to  inelude  all  relevant 
ED  and  TD  effeets. 

4.5.e  CONCLUSION 

The  generalized  standardized  exeitation  approaeh,  whieh  is  a  fast,  universal,  and  rigorous  forward 
modeling  system,  has  been  developed  and  demonstrated.  The  GSEA  is  applieable  to  any  ED  or 
TD  sensor  eonfiguration,  and  to  any  data  set:  novel  waveforms,  multi-axis,  veetor  and  tensor,  or 
magnetie  or  eleetromagnetie  induetion  data,  or  any  eombination  of  magnetie  and  EMI  data. 

The  proposed  system  has  been  tested  against  aetual  data  in  both  the  frequeney  (GEM-3)  and 
time  (EM-63)  domains.  Exeellent  agreements  between  the  GSEA  and  experimental  data  have  been 
demonstrated  here.  The  GSEA  teehnique  eould  be  used  for  building  a  new  type  of  EMI  sensor  as 
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well  as  for  optimal  survey  designing. 

By  using  an  appropriate  dyadic  Green’s  function  the  GSEA  takes  into  account  a  host  medium. 
EMI  response  from  a  magnetically  susceptible  half-space  is  analyzed  versus  sensor  height  and 
frequency.  Numerical  tests  show  that  a  frequency-dependent,  magnetically  susceptible  half-space 
significantly  affects  the  81 -mm  UXO’s  EMI  responses  in  both  ED  and  TD.  This  study  suggests 
that  a  magnetically  susceptible  host  medium  must  be  taken  into  account  in  UXO  discrimination 
problems. 

4.6  NSMS  Modeling  of  GEM-3^+  and  HAP  Method 

The  normalized  surface  magnetic  source  (NSMS)  model  [97-100]  is  a  rigorous  physically  com¬ 
plete  model  of  any  metallic  target’s  interaction  with  an  arbitrary  primary  field  in  the  electromag¬ 
netic  induction  (EMI)  regime.  Based  on  a  number  of  well-defined  measurements  for  a  given  target, 
the  system  solves  for  its  fundamental,  characteristic  response  parameters.  After  target  character¬ 
ization,  the  response  of  that  target  to  any  combination  primary  field  modes  can  be  quickly  and 
accurately  calculated. 

The  NSMS  model  can  be  considered  as  a  generalized  surface  dipole  model  [68,  82-87,  95,  96, 
101,  102],  which  in  a  special  limiting  case  coincides  with  the  infinitesimal  dipole  approach.  In 
the  model,  the  response  of  an  object  to  the  exciting  field  of  a  sensor  is  taken  to  originate  from  a 
layer  of  equivalent  magnetic  dipoles  distributed  over  a  spheroidal  surface  surrounding  the  target. 
Such  a  distribution  can  be  generated  by  spreading  positive  charge  density  Q  over  the  exterior  of 
a  thin  smooth  surface  and  an  identical  negative  charge  on  its  inner  surface,  resulting  in  a  double 
layer  of  charge  separated  by  an  infinitesimal  distance.  This  configuration  introduces  the  proper 
discontinuities  in  the  tangential  components  of  the  magnetic  flux  density  vector  B  but  does  not 
affect  its  normal  component.  (A  single  layer  of  charge,  on  the  other  hand,  would  in  no  way  affect 
the  transition  of  the  tangential  components  of  the  flux  density  but  could  model  a  discontinuity 
in  the  normal  component  of  B.)  Since  there  are  no  free  magnetic  charges  in  nature,  the  normal 
component  of  the  magnetic  flux  density  must  always  be  continuous  across  the  boundary  between 
two  media;  the  secondary  magnetic  field  outside  an  object  can  thus  be  accounted  for  by  a  double 
magnetic  charge  layer  of  density  Q  distributed  on  an  auxiliary  surrounding  surface;  this  density  is 
proportional  to  the  component  of  the  primary  (exciting)  magnetic  field  normal  to  the  surface. 

The  amplitudes  of  the  auxiliary  sources  are  determined  directly,  and  with  great  speed  and  accu¬ 
racy,  by  minimizing  the  difference  between  measured  and  modeled  data  for  a  known  object-sensor 
configuration.  The  charge  density  is  a  property  of  the  target,  and  its  integral  Q  over  the  spheroidal 
surface  constitutes  a  global  magnetic  “capacitance”  of  sorts  that  for  a  given  object  is  invariant  in 
the  sense  that  different  computational  constructs  -  e.g.,  surrounding  surface,  measurement  grid, 
object  location  or  orientation  -  and  different  primary  fields  produce  the  same  value.  Once  the 
source  strengths  are  known,  the  NSMS  model  can  be  used  in  discrimination  processing  that  either 
compares  the  measured  fields  to  those  stored  in  a  library  of  known  objects  -  i.e.,  a  “pattern  match¬ 
ing”  technique  -  or  uses  the  integrated  NSMS  as  a  discriminant  -  i.e.,  “genuine”  inversion.  In 
“pattern  matching”  discrimination  processing,  an  optimization  or  search  algorithm  is  used  to  de¬ 
termine  which  known  UXO  has  the  catalogued  ultrawideband  Q  distribution  that  best  reproduces 
the  signal  received  by  the  survey  sensor.  Used  in  this  way,  the  NSMS  system  is  a  faster  forward 
model  than  other  physically  complete  forward  models  such  as  the  standardized  excitation  approach 
(SEA)  or  the  generalized  SEA  [88,  103-106]. 

We  use  the  total  NSMS,  Q,  of  a  target  as  a  discriminant  in  Sec.  5.  This  Q  may  also  be 
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parametrized  by  a  Pasion  Oldenburg  type  model  and/or  eombined  with  other  parameters,  for  ex¬ 
ample  from  the  dipole  model,  in  order  to  further  eharaeterize  a  target  and  elassify  it  based  on  its 
EMI  response  [107].  When  extraeting  Q  from  training  data  in  order  to  ereate  a  given  target’s  “fin¬ 
gerprint,”  the  loeation  of  the  target  is  generally  know  a  priori.  However,  when  the  target’s  position 
and  orientation  are  not  known,  then  inverse  algorithm  has  to  eombine  inverting  not  only  for  whieh 
target  is  present,  but  where  that  target  is.  This  turns  out  to  be  a  eombination  of  a  linear  and  non¬ 
linear  proeess  and  is  often  time  eonsuming.  To  this  end,  we  have  developed  the  HAP  method  {H 
(magnetie  field),  A  (magnetie  veetor  potential),  and  (magnetie  sealar  potential))  whieh  under 
eertain  assumption  ean  find  the  loeation  and  orientation  of  a  dipolar  souree  without  an  expensive 
inversion  [107-110]. 

While  it’s  easy  to  show  that  simple,  infinitesimal  dipoles  are  frequently  defieient  for  modeling 
a  target’s  general  response  under  diverse  eireumstanees,  the  HAP  approaeh  assumes  that  under  any 
partieular  eireumstanees  a  dipole  representation  is  likely  to  furnish  a  good  first  order  representation 
of  the  response  at  hand.  In  order  to  find  an  estimate  of  the  target’s  loeation  and  orientation  in  a 
rapid,  robust  manner,  we  temporarily  assume  that  the  sourees  of  the  anomalies  being  surveyed  are 
eaused  by  simple  dipole  sourees.  Using  the  magnetie  dipole  approximation  in  the  EMI  regime,  the 
seattered  magnetie  field  H,  magnetie  veetor  potential  A,  and  magnetie  sealar  potential  y/  ean  be 
written  as: 


H 


4kR^  [  R2  ^ 


RxP 

^°4^’ 


RP 

4nR^ 


(4.6.1) 


where  R  =  f  —  f'  with  F  and  f'  being  the  position  of  the  observation  and  of  the  souree,  respeetively. 

Erom  (4.6.1),  it  ean  be  shown  also  that  yr  =  {H  ■  R) /2  and  A  =  iiq{H  x  R)  .  Now,  taking  RxA, 
we  ean  solve  for  R  as 


R 


2yrH-{AxH)/po 


(4.6.2) 


Similarly,  taking  the  eross  produet  of  A  in  (4.6.1)  from  the  left  side,  the  magnetie  dipole  moment 
P  ean  be  expressed  as 

P  =  47tR{Ry/+{AxH)/iio).  (4.6.3) 

Thus,  under  magnetoquasistatie  (MQS)  assumptions  of  a  dipole  souree,  (4.6.2)  and  (4.6.3)  repre¬ 
sent  an  analytieal  method  for  finding  the  loeation,  R  and  orientation,  P  of  the  dipole. 

Currently,  {//,  A,  yr}  are  found  from  EMI  data  through  the  Normalized  Surfaee  Magnetie  Souree 
Method  (NSMS).  An  equivalent,  flat,  two  dimensional  NSMS  surfaee  loeated  in  between  the  mea¬ 
surement  loeations  and  the  target  ean  aeeurately  reproduee  the  response  from  the  target(s)  being 
interrogated  [100].  This  NSMS  surfaee  eonsists  of  a  layer  of  magnetie  eharges  or  dipoles  from 
whieh  {H,A,  yr}  and  henee  {R,P}  can  be  calculated.  Note  that  no  information  regarding  the  tar- 
get(s)  is  required  to  calculate  this  NSMS  surface  layer. 

The  NSMS  model  must  be  adapted  to  each  unique  instrument  one  time.  This  has  been  accom¬ 
plished  for  the  GEM- 3^”''  and  the  results  of  matching  GEM-3^'''  data  to  the  predicted  response 
from  the  NSMS  model  are  shown  in  Eigs.  4.6. 1-4. 6. 3. 
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Figure  4.6.1:  Comparison  of  Hx  from  the  GEM-3^+  and  NSMS  model  predicted  response.  Data  at  90Hz 
over  a  steel  sphere,  inphase  part. 


Benjamin  Barrowes,  Kevin  O’Neill 


-90- 


engineering 

AT  DARTMOUTH 


MM- 15  37,  GEM-3  ^+-  Final  Report 


4  MATERIALS  AND  METHODS:  GEM-3^+  SOFTWARE 
4.6  NSMS  Modeling  of  GEM-3  and  HAP  Method 


Figure  4.6.2:  Comparison  of  Hy  from  the  GEM-3^+  and  NSMS  model  predicted  response.  Data  at  90Hz 
over  a  steel  sphere,  inphase  part. 


Figure  4.6.3:  Comparison  of  from  the  GEM-3^+  and  NSMS  model  predicted  response.  Data  at  90Hz 
over  a  steel  sphere,  inphase  part.  3  levels  with  8 1  points  on  each  level. 
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4.7  Electromagnetic  induction  from  highly  permeable  and  conductive  ellipsoids  under  ar¬ 
bitrary  excitation  -  Application  to  the  detection  of  unexploded  ordnances 

Symmetrical  objects  such  as  spheres  or  spheroids  may  be  useful  as  models  of  some  subsets  of 
possible  targets;  and  for  testing  instruments  they  may  be  particularly  useful  when  one  can  examine 
signals  for  expected  symmetries.  However,  particularly  when  all  vector  components  are  being 
recorded  and  examined,  a  truly  3-D  object  solution  offers  the  most  general  reference  response. 
Heretofore  there  has  been  no  broadband  3-D  analytical  solution  available,  except  by  relatively 
laborious  numerical  methods.  The  approximate  solution  produced  here  is  rapidly  evaluated  under 
diverse  conditions.  It  represents  a  milestone  in  benchmark  cases  that  can  be  used  hereafter  in  a  wide 
variety  of  contexts.  Specifically  in  our  discrimination  processing,  with  this  new  solution  we  can  for 
the  first  time  model  the  responses  of  unseen  objects  quickly  without  having  to  assume  restrictive 
symmetries.  We  now  have  a  3-D  solution  to  guide  analysis,  statistical  learning  machines,  etc,  for 
inferring  overall  geometrical  properties  such  as  volume  and  aspect  ratios,  from  signals  across  the 
breadth  of  the  EMI  band. 

The  secondary  field  produced  by  three-dimensional  highly  permeable  and  conductive  objects 
is  computed  in  the  electromagnetic  induction  regime,  with  the  purpose  of  modeling  unexploded 
ordnances  and  surrounding  clutter.  The  analytical  formulation  is  based  on  the  ellipsoidal  coordi¬ 
nate  system,  able  to  model  real  three  dimensional  geometries  as  opposed  to  bodies  of  revolutions 
like  within  a  spheroidal  approach.  At  the  frequencies  of  interest  (tens  of  Hertz  to  hundreds  of 
kilo-Hertz),  the  conduction  currents  in  the  soil  are  negligible  and  the  fields  are  computed  in  the 
magneto-quasistatic  regime  based  on  the  Laplace  equation.  Inside  the  objects,  where  the  wave 
equation  governs  the  field  distribution,  the  currents  are  assumed  to  have  a  small  penetration  depth, 
allowing  for  the  analytical  simplification  of  the  field  components,  which  become  decoupled  at  the 
surface.  This  approximation,  valid  across  the  entire  frequency  spectrum  because  of  the  high  per¬ 
meability  and  conductivity,  avoids  the  necessity  of  using  ellipsoidal  wave  functions  and  results  in 
a  considerable  saving  of  computational  time.  Numerical  results  compare  favorably  with  numerical 
and  experimental  data  .  Finally,  the  optimization  approach  used  to  match  our  numerical  predictions 
with  experimental  data  demonstrates  the  possibility  of  remotely  inferring  the  material  properties 
of  objects. 

4.7.a  Introduction 

One  theoretical  approach  models  the  UXO  as  a  (possibly)  magnetic  and  (possibly)  conductive 
spheroid  [2,  111-113].  The  model  is  therefore  mathematically  more  involved  but  offers  impor¬ 
tant  flexibility  in  terms  of  geometries  (from  elongated  to  flat  objects  within  the  same  formalism), 
material  parameters,  and  distance  to  sensor. 

In  this  section,  we  extend  the  second  approach  to  ellipsoidal  geometries  with  the  underlying 
motivation  of  providing  an  improved  physical  model  of  the  clutter  surrounding  UXO,  directly  use¬ 
ful  for  example  in  statistical  and  optimization  methods  such  as  in  [114-118].  Clutter  is  indeed 
most  often  composed  of  random  three  dimensional  (3D)  pieces  of  steel  which  are  jamming  the 
measured  signals  and  are  in  part  responsible  for  the  high  false  positive  alerts.  Their  proper  mod¬ 
eling  is  therefore  essential  and  the  ellipsoidal  coordinate  system  provides  a  general  framework  to 
do  so,  allowing  for  the  modeling  of  real  3D  geometries  while  retaining  the  speed  advantage  of 
analytical  calculations,  fundamental  for  example  in  inversion  algorithms  [1 19][120][121-123].  In 
addition,  we  use  in  this  paper  the  small  penetration  approximation  (SPA),  justified  by  the  major 
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steel  component  of  clutter  objects.  SPA  was  originally  developed  as  a  high  frequency  approxima¬ 
tion  for  spheroidal  geometries  [113]  under  the  assumption  that  the  currents  induced  on  the  surface 
of  the  object  penetrate  only  slightly  into  its  volume.  Remarkably,  however,  SPA  has  been  shown 
to  yield  accurate  results  over  the  entire  EMI  frequency  range  for  highly  permeable  and  highly  con¬ 
ductive  objects,  such  as  those  made  of  steel.  This  success  has  prompted  various  implementations 
of  the  SPA,  both  analytical  [113]  and  numerical  [124,  125].  In  this  paper,  we  generalize  the  SPA 
to  ellipsoidal  geometries  and  confirm  the  good  agreement  obtained  over  the  entire  frequency  range 
by  comparing  our  predictions  with  both  well  established  numerical  results  and  with  measurement 
data. 


The  configuration  under  study  is  depicted  in  Fig.  4.7.1:  a  homogeneous  ellipsoid  of  semi-axes 
(a  >  b  >  c)  along  the  directions  (x,y,£)  (see  the  Appendix  for  mathematical  details),  of  perme¬ 
ability  /ii  =  lOO/lo  and  conductivity  Oi  =  10^  [S/m]  is  embedded  in  a  non-magnetic  and  non- 
conductive  soil  which,  at  the  frequencies  of  interest,  is  essentially  transparent  and  modeled  as 
free-space.  A  sensor  located  above  the  object  emits  a  primary  field  H"  and  measures  the  secondary 
field  H  produced  by  the  object.  We  subsequently  use  two  models  for  the  primary  field:  a  sim¬ 
ple  uniform  field,  valid  for  small  or  deeply  buried  objects,  and  the  field  produced  by  the  GEM-3 
instrument  [1],  more  realistic  for  large  or  close  objects. 

In  the  EMI  regime,  the  magnetic  field  outside  the  ellipsoid  is  irrotational  and  can  therefore  be 
written  in  terms  of  a  potential  function  which  obeys  the  Eaplace  equation.  Denoting  the  primary 
and  secondary  potentials  as  C/"  and  U^,  respectively,  the  total  magnetic  field  in  the  region  exterior 
to  the  ellipsoid  is  written  as 

H2  =  H‘’ =  -yU^ .  (4.7.1) 


In  the  interior  region  to  the  object,  the  magnetic  field  H\  obeys  the  wave  equation 

Vx  VxFi-kf77i  =0  (4.7.2) 

where  k\  =  ittXTi/ii .  This  limiting  form  of  k\  that  applies  in  this  parameter  space  means  that  (4.7.2) 
leads  in  fact  to  a  diffusion  equation  representing  magnetic  diffusion,  though  it  retains  the  form  of 
a  wave  equation.  The  two  quantities  H\  and  H2  are  related  by  the  usual  boundary  conditions. 

For  details  on  the  ellipsoidal  coordinate  system,  the  separation  of  the  wave  equation  in  ellip¬ 
soidal  coordinates,  and  the  derivation  of  the  secondary  field  in  the  ellipsoidal  coordinate  system 
as  a  function  of  the  type  of  incidence  (uniform  and  GEM-3)  within  the  SPA  approximation,  see 
[126]. 

The  ellipsoidal  coordinate  system  is  defined  based  on  the  definition  of  the  master  ellipsoid  of 
semi-axes  dimensions  (a,h,c)  along  the  (i,y,£)  directions  respectively  (a>  b>  c): 


222 
y  z 

- h  - 1 - =  1 

^2^  u2^  r-2  ^  • 


(4.7.3) 


b^ 

Upon  defining  the  variable  k^  =  and  —  b^  [127],  three  families  of  orthogonal  sur¬ 

faces  can  be  obtained  [128],  used  to  defined  the  ellipsoidal  coordinate  system  (p,  /i,  v) : 


X 

^^p2 

~2 


■h^ 


+ 


y 


p2  —  k^ 

,2 


p2  -h'^  ^2  _  1^: 


=  1, 


=  1, 


(4.7.4a) 


(4.7.4b) 


V" 


V" 


k^ 


=  1, 


(4.7.4c) 
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Figure  4.7.1:  Configuration  of  the  problem:  an  ellipsoid  of  material  parameters  (/Ti  =  100/Xo,  CJi  = 
10^  [S/m])  is  surrounded  by  free-space  {1X2  =  iM),  02  ~  0)  and  is  excited  by  a  GEM-3  instrument  [1]  modeled 
as  two  current  loops  (figure  not  to  scale).  The  semi-axes  of  the  ellipsoid  are  denoted  by  a,  b,  and  c  (with 
a  >  b  >  c)  along  the  x,  y  and  £  directions,  respectively  (see  the  Appendix  for  mathematical  details  on  the 
ellipsoidal  coordinate  system). 


and  where  >k^  >  >h^  >\^  >  0. 

The  magnetic  field  is  obtained  from  the  gradient 


yu 


^  \  dU  .  \  dU  ^  I  dU 

P~i — 5 - ^  P~i — "5 - ^  ^7 — 5 — 

hp  op  hp  op  hy  ov 


(4.7.5) 


and  is  written  as 


h"  =  -yu‘’  = 


Y^bl!’h{p,p\v) 

n,p 


for  the  primary  field. 


Y,b[^h{p,py) 

n,p 

(4.7.6) 
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Introducing  the  field  expansions  and  separating  the  terms  funetions  of  and  yields 


ih—  [l'{a)E{a)+I{a)E'{a)]E{ii,v) 


n,p 


Ml 


hp  d 
hphy  d jLi 


hphy  d\^ 


I(a)E(a,il,  v')  -  ^/v/(a)E(a.fJ,  v") 


I* 

n,p 


ip) 

n 


iki^E{p',iiX 

Ml 


+ 


hr 


h^Kdll  y,,yp)E(ay,v)  + ^^f^E(a,)l".v) 


+ 


hr 


h^hy  dv  \h 


hr 


'-^fy]E{a,py)  +  '-j^fyE{a,py 


(4.7.7) 


and  it  is  understood  that  this  expression  is  evaluated  at  the  surfaee  of  the  ellipsoid  p  =  a.  It  should 
be  noted  that  most  of  the  terms  in  Eq.  (4.7.7)  are  straightforward  to  compute:  the  derivatives  and 
seeond  derivatives  of  the  Lame  funetions  with  respect  to  p  and  v  do  not  present  any  extreme 
diffieulty. 


4.7.b  Secondary  fields  due  to  uniform  incidence 


The  first  situation  we  shall  eonsider  is  the  one  of  a  uniform  ineidenee.  The  latter  is  not  merely  a 
theoretieal  idealization  sinee  objeets  that  are  far  from  the  GEM-3  sensor,  or  those  that  are  small 
(in  the  order  of  a  few  eentimeters),  do  not  see  a  very  inhomogeneous  primary  exeitation.  The 
assumption  of  uniform  ineidenee  in  these  eases  is  therefore  exeellent. 

Eor  the  sake  of  later  comparison  with  known  solutions,  we  first  suppose  the  primary  field  to  be 
in  the  [xz)  plane,  although  the  general  ease  ean  be  treated  just  as  easily.  A  uniform  field  derives 
from  a  simple  linear  potential  in  x  and  z,  whieh  ean  also  be  expanded  in  the  ellipsoidal  eoordinate 
system  as 

U°  =  HQX?,me+HQZ?,me  =Y,b\^’hl!’\p,py  (4.7.8) 

n,p 

where  0  =  0  for  axial  exeitation  and  6  =  nj2  for  transverse  exeitation,  and  where  Hq  is  the  am¬ 
plitude  of  the  uniform  field.  Realizing  that  x  and  z  ean  be  direetly  expressed  in  terms  of  Lame 
functions  of  the  first  kind  with  n  =  1  we  write 


,  „  HosinO  ^(1) ,  ,  HqcosO  ^(3) ,  , 

^  =^^£r^(p,M,v)+-^L==Er^p,M,v) 


hk 


kVk^  —  h?- 


Sp) 


so  that  the  bn  eoeffieients  of  Eq.  (4.7.7)  are  directly  given  by 


Hq  sin  0 
hk 


0,  b 


(3) 


HocosO  (p) 

kyyy?' 


0. 


(4.7.9) 


(4.7.10) 


The  expansion  of  the  seeondary  field  also  deserves  attention  before  proeeeding  with  the  numer- 
ieal  solution  of  the  b\{’^  .  Indeed,  there  exist  an  analytieal  solution  to  the  external  potential  due  to  a 
homogeneous  ellipsoid  [129]: 


y{.x,y,z)  =  Kb 


3^ 


2 


—  h^ 


P  d5 

s^-k^_  yJs^-hP-y/s^-k'^ 


(4.7.11) 
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where  Kq  is  a  eonstant.  We  see  that  the  first  term  of  Eq.  (4.7.1 1)  is  direetly  proportional  to  Kq  (i.e. 
n  =  0),  while  the  sueeessive  three  terms  are  direetly  proportional  to  Ki,  Li,  and  M\  (i.e.  n  =  1), 
respeetively.  The  potential  of  Eq.  (4.7.11)  ean  therefore  be  expressed  with  exaetly,  and  only,  four 
ellipsoidal  term,  direetly  related  to  the  demagnetization  faetors  diseussed  in  [130,  131].  In  addition, 
sinee  the  n  =  0  term  is  disearded  in  our  solution  (the  ellipsoid  being  uneharged),  the  expansion  of 

the  seeondary  potential  is  redueed  to  n  =  1  only  {p  =  1,2,3),  and  =  0.  Moreover,  in  the 

(2) 

speeifie  ease  of  a  uniform  primary  field  in  the  (xz)  plane,  B\  =0  and  we  are  left  with  only  two 
unknowns  to  determine. 

We  first  apply  our  method  to  the  eomputation  of  the  seeondary  field  produeed  by  a  series  of 
spheroidal  geometries,  under  axial,  transverse,  and  oblique  ineidenees.  The  main  purpose  here  is  to 
being  able  to  eompare  the  results  with  those  obtained  in  our  previous  work  [2] .  The  eorrespondenee 
between  the  elongation  of  the  spheroids  and  the  size  of  the  ellipsoids  is  given  in  Table  2.  The  axial 


Spheroidal 

Ellipsoidal  dimensions 

Measured  field 

elongation 

a>  b>  c  (along  x,  y,  z) 

Fig.  4.7.2(top) 

Fig.  4.7.2  (bottom) 

e  =  6 

a  =  1.5  cm,  b  =  c  =  0.25  cm 

Hz 

e  =  3 

a  =1.5  cm,  ft  =  c  =  0.5  cm 

H, 

Hz 

e  =  1 

a  =  ft  =  c  =  1.5  cm 

H 

c  or  Hz 

e  =  1/3 

a  =  ft  =  1 .5  cm,  c  =  0.5  cm 

H, 

e  =  1/6 

a  =  ft  =  1 .5  cm,  c  =  0.25  cm 

H, 

Table  2:  Correspondence  between  the  elongations  of  spheroids  and  the  size  of  the  ellipsoids  as  shown  in 
Fig.  4.7.1  and  defined  by  Eq.  (4.7.3). 


and  transverse  eases  are  shown  in  Fig.  4.7.2  for  spheroids  loeated  10  em  below  the  sensor  and 
with  elongations  e  G  {6,3, 1,1/3, 1/6}  {e  >  1  eorresponding  to  prolate  objeets  while  e  <  1  to 
oblate  ones).  The  overall  responses  of  the  objeets  as  funetion  of  frequeney  follow  the  now  well- 
understood  trend:  at  low  frequeneies  the  seeondary  magnetie  field  is  aligned  with  the  exeiting  field 
while  at  high  frequeneies,  the  indueed  eurrents  are  eonfined  to  the  surfaee  of  the  ellipsoid  and  are 
out  of  phase  by  n  due  to  Eenz’  law.  The  imaginary  part  of  the  field  therefore  vanishes  at  both  ends 
of  the  frequeney  speetrum  and  peaks  at  an  intermediate  frequeney  eorresponding  to  the  resonanee 
of  the  objeet. 

In  the  axial  ease,  the  magnetie  field  sees  more  objeet  surfaee  with  prolate  geometries,  produe- 
ing  a  stronger  response  than  oblate  geometries,  while  this  tendeney  is  reversed  under  transverse 
exeitation  for  the  same  reason.  The  position  of  the  peak  in  the  imaginary  part  is  direetly  related 
to  the  resonanee  of  the  objeet,  also  elosely  related  to  the  size  of  the  objeet  along  the  direetion  of 
the  primary  field.  In  the  transverse  ease,  elongations  of  e  G  (1/6, 1/3, 1}  present  the  same  size 
along  the  primary  field  and  therefore,  a  similar  loeation  of  the  resonant  peak.  Prolate  objeets,  on 
the  other  hand,  have  a  varying  size  along  this  direetion  and  thus,  different  resonant  peaks.  Similar 
arguments  ean  be  applied  to  explain  the  results  obtained  under  the  axial  exeitation. 

The  first  eonelusion  when  examining  Fig.  4.7.2  is  that  the  agreement  over  the  entire  frequeney 
range  is  reasonably  good  in  most  of  the  eases,  thus  proving  that  the  SPA  is  an  adequate  approxima¬ 
tion  for  highly  magnetie  and  eonduetive  objeets,  as  already  shown  in  [113].  The  prolate  eases  are 
seen  to  usually  agree  better  with  the  referenee  markers  [2],  and  even  more  so  under  the  transverse 
exeitation.  The  worst  agreement  is  obtained  for  oblate  geometries  under  transverse  exeitation,  with 
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increased  disagreement  as  the  objeets  flatten,  in  agreement  with  the  conclusions  reported  in  [125]. 

Although  not  shown  for  the  sake  of  clarity,  simulations  have  been  extended  up  to  extreme  elon¬ 
gations  of  10  and  1/10  and  eonfirm  these  findings.  A  possible  explanation  is  that  under  these 
eireumstanees,  a  very  large  portion  of  the  objeet  surfaee  has  a  very  high  radius  of  eurvature  so  that 
a  dominant  radial  deeay  approximation  of  the  eurrents  may  beeome  less  aeeurate. 

Fig.  4.7.3  presents  the  results  for  a  eonfiguration  similar  to  the  one  eorresponding  to  Fig.  4.7.2 
but  where  the  objeet  is  rotated  at  45  degrees.  Similar  eonelusions  to  the  previous  ease  ean  be 
drawn:  prolate  spheroids  are  well  modeled  whereas  oblate  ones  present  the  highest  mismateh,  the 
latter  being  more  pronouneed  for  more  extreme  geometries.  The  disagreement  for  eG  {1/3, 1/6} 
is  here  seen  to  be  more  evident  at  low  frequeneies,  while  at  high  frequeneies  (where  the  SPA  is 
more  and  more  valid)  the  agreement  is  usually  reasonably  good. 

Finally,  Fig.  4.7.4  shows  the  seeondary  field  produeed  by  a  real  ellipsoid  of  dimensions  1  emx  2  emx  3  em 
along  the  x,  y,  and  £  direetions,  respeetively,  under  an  axial  exeitation.  The  depth  and  material  pa¬ 
rameters  are  otherwise  identieal  to  the  previous  oases.  The  referenoe  solution  is  obtained  from 
the  method  of  auxiliary  sourees  (MAS)  presented  in  [132], [3],  and  has  been  trunoated  at  1  MHz 
beeause  of  the  fine  mesh  and  long  computation  time  required  at  higher  frequencies.  At  low  fre¬ 
quencies,  however,  the  MAS  is  very  aeeurate  and  shows  that  our  method  suffers  from  some  errors 
which  we  attribute  to  the  choice  of  the  looation  of  the  matching  points.  Nonetheless,  the  differenoe 
between  the  seeondary  field  produced  by  an  object  of  dimensions  1  emx 2  emx  3  cm  (the  present 
ellipsoid)  and  1  emx  1  emx  3  cm  (a  spheroid  of  elongation  e  =  3)  is  quite  significant  and  allows 
us  to  conclude  that  away  from  the  lowest  frequeneies,  the  agreement  between  the  two  methods  is 
reasonably  good. 


4.7.C  Secondary  fields  due  to  GEM-3  incidence 


Although  the  assumption  of  uniform  primary  field  is  suffieient  in  many  eases,  it  is  still  beneficial  to 
model  the  primary  field  generated  in  real  field  measurements,  typically  as  obtained  from  the  GEM- 
3  instrument.  The  GEM-3  is  a  broadband  eleetromagnetie  sensor  manufaetured  by  Geophex  [1], 
whose  transmitting  sensor  head  eonsists  of  two  eurrent  loops  of  radii  of  20  cm  and  11.074  cm. 
These  two  current  loops  are  arranged  as  bueking  eoils  with  the  smaller  loop  plaeed  eoneentrioally 
in  the  interior  of  the  larger  loop  (schematieally  illustrated  in  Fig.  4.7.1).  The  two  currents  are 
optimized  so  that  the  primary  magnetie  field  at  the  center  of  the  instrument  is  minimal  and  the 
measurements  elosely  eorrespond  to  the  secondary  field  H\ 

The  analytieal  expression  of  the  axial  and  transverse  magnetie  field  eomponents  produeed  by  a 
eireular  eurrent  loop  in  a  eylindrieal  eoordinate  system,  based  on  the  eomplete  elliptieal  integrals 
of  the  first  and  seeond  kind,  ean  be  found  in  [74,  75,  133]  and  reads: 


Mo/o _ jz-zo) _ 

2;r  r[(r-hro)2-h(z-Zo)2]i/2 


-K{kc)  + 


r'^  +  rl  +  jz-Zo)^ 

(r-ro)2-h(z-zo)2 


E{kc) 


(4.7.12a) 


Bz{r,^,z) 


_ Mo/o _ 

2;r[(r  4- ro)2 -I- (z  -  Zo)2]  1/2 


K{kc) 


r^-rl  +  {z-Zof 
(r-ro)2  +  (z-Zo)^ 


E{kc)  , 


(4.7.12b) 


where  ro  is  the  radius  of  the  eurrent  loop  loeated  in  the  z  =  zo  plane,  4  is  the  eurrent  flowing  in 
the  ^  direetion,  and 


kc  = 


4ar 


(r-ha)2-|-  (z-Zo)^' 


(4.7.13) 
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Axial  excitation 


Transverse  excitation 


Figure  4.7.2:  Real  and  imaginary  parts  of  the  secondary  field  produced  by  spheroids  of  elongations 
shown  in  the  legend,  buried  10  cm  below  the  sensor,  and  under  axial  (upper  panels)  and  transverse  (lower 
panels)  excitations.  The  material  parameters  are  jUi  =  lOOjUo  and  0\  =  10^  [S/m].  The  sphere  case  (denoted 
by  ’e=l  (n)’)  has  been  normalized  by  1/3  for  display  convenience.  The  markers  are  obtained  from  the 
method  presented  in  [2] . 


Benjamin  Barrowes,  Kevin  O’Neill 


-98- 


engineering 

AT  DARTMOUTH 


MM- 15  37,  GEM-3  ^+-  Final  Report 


4  ELLIPSOIDAL  SPA 


Figure  4.7.3:  Similar  to  Fig.  4.7.2  for  an  object  buried  at  30  cm  and  rotated  by  45  degrees.  Upper  panel:  Hx- 
Lower  panel:  Note  that  for  this  specific  case  only,  the  ranges  have  been  set  to  v  G  [0,0. 7]h  for  prolate 

spheroids  and  /i  G  [0, 0.7]  {k  —  h)  +  h  for  oblate  spheroids. 
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Figure  4.7.4:  Real  and  imaginary  parts  of  the  secondary  field  produced  by  an  ellipsoid  of  dimensions  1  cm 
in  X,  2  cm  in  y,  and  3  cm  in  z-  The  material  parameters  are  /ii  =  IOO/Iq  and  ai  =  10^  [S/m],  and  the  ellipsoid 
is  buried  10  cm  below  the  sensor.  The  reference  markers  have  been  obtained  using  the  method  of  auxiliary 
sources  [3]. 


K  and  E  are  the  eomplete  elliptieal  integral  funetions  of  the  first  and  seeond  kind,  respeetively. 
Note  that  this  model  agrees  with  the  alternative  methods  of  using  a  numerieal  integration  or  the 
applieation  of  the  Biot-Savart’s  Law,  and  has  the  added  benefit  of  being  eomputationally  less  ex¬ 
pensive. 

The  first  step  in  modeling  a  GEM-3  ineidenee  is  to  determine  the  expansion  eoeffieients  of 
the  primary  field  in  Eq.  (4.7.6),  simply  done  using  a  point  matehing  teehnique.  Henee  an  ellipsoid 
is  first  defined  in  spaee,  yielding  a  series  of  (x,y,z)  eomponents.  At  those  loeations,  the  magnetie 
field  is  eomputed  from  Eq.  (4.7.12),  and  the  field  eomponents  are  transformed  into 
Eq.  (4.7.6)  is  then  used  aeross  many  points  to  eompute  the  matrix  of  an  over-determined  system, 
from  whieh  the  are  solved  for  in  the  least  square  sense.  In  doing  so,  we  find  that  the  eoeffieients 
are  more  stable  when  the  matehing  is  done  on  Hy  for  prolate-type  objects  and  on  Hp  for  oblate-type 
objects.  The  matching  using  Hp,  reasonably  good  in  both  cases,  can  also  be  used  for  the  purpose 
of  comparison. 

Due  to  the  non-uniform  nature  of  the  primary  field,  terms  of  higher  degrees  than  n  =  I  need  to 
be  included  in  the  ellipsoidal  expansion  and  therefore,  all  the  coefficients  need  to  be  solved 
for  simultaneously.  The  lack  of  direct  orthogonality  between  the  various  terms  of  Eq.  (4.7.7)  results 
in  a  somehow  unstable  numerical  solution:  the  influence  of  some  modes  {n,p)  can  leak  into  other 

modes  and  degrade  the  accuracy  of  the  corresponding  .  This  is  typically  the  case  when  some 
are  very  small,  and  need  to  be  solved  as  zero  numerically:  although  the  results  are  numerically 
small  as  expected,  the  non-zero  B^f^  can  be  slightly  modified,  introducing  errors  in  the  secondary 
field  values.  In  the  previous  examples,  this  phenomenon  has  been  bypassed  by  using  a  physical 
argument  to  force  some  coefficients  to  zero  and  to  solve  for  the  non-zero  ones  only.  In  the  non- 
uniform  incidence  case,  however,  such  reasoning  cannot  in  general  be  applied  and  all  the  B^^'^  need 
to  be  solved  for  simultaneously. 

In  order  to  minimize  the  numerical  leakage  problem,  we  resort  to  solving  for  the  b\{’^  in  two 
steps.  Eirst  we  solve  for  all  the  B^j’^  in  a  usual  manner.  Second,  the  B^^^  are  used  to  compute 
the  field  at  the  sensor  location.  The  contribution  of  each  mode  is  compared  and  those  modes 
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Figure  4.7.5:  Real  and  imaginary  parts  of  the  secondary  field  produced  by  a  spheroid  of  elongation  e  =  3 
(left  panel)  and  e  =  1/3  (right  panel)  buried  10  cm  below  the  GEM-3  instrument.  The  material  parameters 
are  /ii  =  IOO/Iq  and  ai  =  10^  [S/m].  The  primary  field  is  modeled  using  Eqs.  (4.7.12).  The  reference 
markers  (marked  ‘ref’  in  the  legend)  are  obtained  from  the  method  presented  in  [2]. 


contributing  less  than  a  threshold  (1/1000  of  the  maximum  eontribution  in  our  ease)  are  disearded. 
The  are  then  solved  anew  only  for  the  eontributing  modes.  Note  that  this  seleetion  works  best 
if  it  is  performed  at  the  first  frequeney,  and  kept  identieal  for  subsequent  frequeneies. 

The  results  for  a  GEM-3  axial  exultation  of  a  prolate  {e  =  3)  and  an  oblate  {e  =\  /3)  spheroid 
are  shown  in  Fig.  4.7.5.  The  neeessity  of  ineluding  more  modes  in  the  ellipsoidal  expansion  is 
apparent,  although  the  number  of  modes  is  still  signifieantly  lower  than  under  a  spheroidal  eoordi- 
nate  approaeh.  The  agreement  obtained  here  with  up  to  n  =  3  for  the  prolate  ease  and  up  to  n  =  2 
for  the  oblate  ease  is  seen  to  be  good.  The  transverse  exultation  exhibits  a  similar  good  agreement 
and  is  not  shown  for  the  sake  of  brevity. 

Finally,  measurement  data  are  taken  as  a  last  validating  element  of  our  analytieal  approaeh. 
Measurements  were  performed  using  the  GEM-3  instrument  along  the  three  axes  of  a  steel  ellipsoid 
of  size  15  emx7.5  emx3.75  em  in  the  f,  y,  and  z  direetions,  respeetively  (yielding  an  aspeet  ratio 
of  1:2:4),  The  magnetie  permeability  and  eleetrie  eonduetivity  of  the  steel  are  a  priori  unknown 
and  ean  be  optimized  for  when  matehing  the  experimental  data  to  the  theoretieal  ones. 

The  eomparisons  between  measurements  and  analytieal  results  are  illustrated  in  Fig.  4.7.6  for 
the  data  eolleeted  along  the  three  axes.  The  eaption  provides  information  on  the  sealing  faetors  7 
optimized  when  matehing  the  data  as  well  as  on  the  values  of  the  permeability  and  the  eonduetivity 
for  eaeh  ease.  It  ean  be  seen  that  overall,  the  numerieal  approaeh  is  very  sueeessful  at  predieting 
the  EMI  response  of  the  ellipsoidal  objeet  aeross  the  entire  frequeney  range.  This  good  agreement 
represents  an  important  validation  of  the  method  presented  in  this  paper. 

A  more  speeifie  examination  of  the  results  reveals  both  good  and  less  good  aspeets,  however. 
First,  the  sealing  faetor,  whieh  eonverts  the  units  of  GEM-3  measurements  into  [A/m],  is  relatively 
eonstant  in  the  three  independent  optimizations  results  and  yields  7~  4.5  x  10^.  Along  with  the 
overall  good  agreement  in  the  magnetie  field,  this  is  an  important  result  whieh  eonfirms  the  ae- 
euraey  of  our  measurement  teehnique.  Seeond,  the  first  and  seeond  eases  in  Fig.  4.7.6  produee 
very  similar  results  for  the  permeability  and  the  eonduetivity,  whieh  is  another  important  positive 
result.  The  numerieal  values  eonfirm  the  expeetation  that  steel  is  a  highly  permeable  and  eondue- 
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Case  A  Case  B 


Figure  4.7.6:  Simulated  (solid  lines)  and  experimentally  measured  (circles)  real  and  imaginary  parts  of  the 
magnetic  field  along  the  three  axes  of  an  ellipsoid:  case  A  is  along  the  longest  dimension  (15  cm),  case  B  is 
along  the  second  longest  dimension  (7.5  cm),  case  C  is  along  the  shortest  dimension  (3.75  cm).  The  location 
of  the  resonant  frequency  in  the  three  cases  follow  this  geometrical  progression,  as  expected.  The  optimized 
parameters  in  the  respective  cases  are:  (/ii  =  140/Xo>  =  4.5  x  10®  [S/m],  7  =  4.7  x  10^)  (Hi  =  150/Xo, 

ai  =  5  X  10®  [S/m],  7  =  4.5  x  10^)  (^ti  =  20no,  ai  =  1.2  x  10®  [S/m],  7  =  4.2  x  10^).  The  dashed  line  in 
case  C  has  been  obtained  with  (/Ti  =  100/Io^  CJi  =  1  x  10^  [S/m],  7  =  4.2  x  10^).  The  numerical  data  have 
been  obtained  with  tree  modes  in  all  cases  (N  =  3). 
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live  material,  which  is  a  necessary  condition  for  our  SPA  approximation  to  be  trustable  across  the 
entire  frequency  range.  More  importantly,  these  results  show  that  we  are  able,  to  some  extent,  to 
infer  the  material  properties  of  an  object  in  a  totally  remote  and  non-invasive  way.  Third,  it  is  seen 
that  the  last  situation  in  Fig.  4.7.6  stands  out  from  the  other  two,  with  different  estimates  of  the 
permeability  and  of  the  conductivity.  We  believe  that  this  case  is  tantamount  to  cases  of  very  flat¬ 
tened  oblate  spheroids  studied  previously  (e  =  1  / 6),  where  the  currents  flowing  around  the  object 
see  a  very  large  radius  of  curvature  at  all  points  around  the  circumference,  producing  less  accurate 
results  within  an  SPA  approximation.  Moreover,  the  cost-function  associated  with  this  third  situa¬ 
tion  presents  many  local  minima  so  that  the  optimized  values  jli  =  20/lo  and  ai  =  1.2  x  10^  [S/m] 
have  to  be  understood  as  being  one  possible  solution.  In  fact,  the  dashed  lines  correspond  to  the 
values  of  the  magnetic  field  obtained  with  jii  =  100/io  and  (7i  =  1  x  10^  [S/m]  (the  initial  guess 
in  the  optimization  procedure),  and  yield  a  reasonably  good  match  with  the  measured  data  as  well. 
Consequently,  we  can  consider  this  last  case  as  less  reliable  than  the  other  two,  but  yet  constitutes 
an  important  situation  to  consider  in  future  work. 

4.7.d  Conclusion 

The  motivation  of  modeling  clutter  surrounding  unexploded  ordnances  (UXOs)  as  three-dimensional 
steel  objects  in  the  electromagnetic  induction  regime  has  prompted  the  study  of  the  ellipsoidal  co¬ 
ordinate  system  presented  in  this  paper.  The  hybrid  approach,  combining  a  static  formulation  based 
on  the  Laplace  equation  in  the  region  exterior  to  the  object  and  a  full-wave  approach  based  on  the 
wave  equation  in  the  region  interior  to  the  object,  has  been  formulated  based  on  the  Lame  functions 
and  the  small  penetration  approximation  (SPA).  The  latter,  originally  derived  as  a  high  frequency 
approximation,  is  shown  to  be  valid  across  the  entire  frequency  spectrum  of  the  GEM-3  instru¬ 
ment  for  highly  permeable  and  highly  conductive  objects,  in  agreement  with  previous  conclusions. 
Importantly,  the  SPA  bypasses  the  necessity  of  using  ellipsoidal  wave  functions  inside  the  object 
and  results  in  a  considerable  reduction  of  computational  time,  yielding  an  algorithm  that  produces 
broadband  secondary  magnetic  field  responses  in  seconds  under  uniform  excitation,  and  in  tens  of 
seconds  under  GEM- 3  incidences. 

Numerical  results  have  been  favorably  compared  with  three  independent  results:  those  obtained 
with  a  well-established  method  for  the  study  of  spheroidal  geometries  (treated  here  as  a  special 
case  of  an  ellipsoidal  geometry),  those  obtained  with  the  method  of  auxiliary  sources  for  metallic 
objects  (where  the  ellipsoidal  geometry  is  a  particular  case  of  a  general  three-dimensional  shape), 
and  those  obtained  by  GEM-3  measurements  of  a  manufactured  ellipsoid.  The  good  agreement  of 
all  these  results  prompts  us  to  conclude  that  the  method  presented  in  this  paper  is  both  fast  and 
accurate,  and  can  therefore  be  used  in  subsequent  inversion  algorithms  for  the  classification  of 
buried  UXOs  in  clutter  contaminated  soils  as  well  as  for  the  estimation  of  material  parameters  in  a 
non-invasive  way. 
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5  Results  and  Discussion  -  GEM-3^+ 

The  GEM-3^+  can  be  used  in  several  different  measurement  modes,  namely,  dynamic  mode, 
static  mode,  and  hybrid  mode.  These  correspond  to  what  type  of  information  the  user  wants  to 
know.  Dynamic  mode  may  be  used  mostly  for  the  detection  of  metallic  objects,  while  static  mode 
could  give  higher  fidelity  lower  noise  data  for  discrimination.  The  hybrid  mode,  which  allows  for 
measurements  taken  while  the  sensor  head  is  in  motion,  and  relies  on  post  processing  to  extract 
individual  data  points,  may  turn  out  to  be  the  most  commonly  used  mode. 

In  static  mode,  the  GEM-3^”''  would  be  placed  on  a  certain  location  and  held  immobile  while 
one  data  shot  is  taken.  Because  the  sensor  is  not  moving,  the  data  may  be  integrated  for  a  very 
accurate  low  noise  data  point.  Several  of  these  in  a  rough  grid  would  aid  inversion  efforts  toward 
discrimination.  The  main  drawback  to  this  mode  is  that  the  data  acquisition  sequence  would  require 
a  longer  time  than  data  acquired  in  a  dynamic  or  hybrid  mode. 

Acquiring  data  in  the  the  dynamic  mode  would  allow  the  user  to  hold  and  rotate  the  sensor 
head  on  the  handle  much  more  freely  without  having  to  pause  between  the  data  shots.  At  a  10  Hz 
measurement  rate,  substantial  amounts  of  data  could  be  acquired  in  a  short  period  of  time.  One 
concern  when  acquiring  data  in  dynamic  mode  is  the  effects  of  sensor  head  movement  on  data 
quality.  Another  concern  would  be  data  stability  of  each  10  Hz  measurement  without  integrating 
to  reduce  noise.  If  the  sensor  head  is  not  moved  too  quickly,  movement  issues  can  be  kept  to  a 
minimum.  Eurthermore,  when  a  post  processing  step  is  added  (resulting  in  measuring  in  the  hy¬ 
brid  mode)  data  acquired  while  the  sensor  head  was  moving  too  rapidly  could  be  flagged  and  or 
discarded.  To  address  stability  concerns  of  the  10  Hz  measurements,  we  plotted  the  standard  devi¬ 
ation  of  a  few  seconds  of  data  acquired  while  the  instrument  was  not  moving  (see  Eig.  5.1.11).  The 
GEM-3^'*"  has  exceptional  stability  over  short  time  periods  even  at  the  10  Hz  measurement  rate 
and  especially  in  the  mid  frequency  band.  Note  that  we  normally  only  regard  signals  greater  than 
100  ppm  as  significant.  During  post  processing,  measurements  taken  sufficiently  close  together 
(for  example,  when  the  instrument  was  held  immobile  during  dynamic  acquisition  mode)  could  be 
averaged  together  leading  to  higher  quality,  higher  SNR  data.  All  of  these  data  collection  modes 
rely  explicitly  on  a  very  accurate  positioning  system  such  as  the  beacon  system  described  above 

5.1  Initial  Data  and  Reasonableness  Checks 

We  have  taken  some  initial  data  shots  of  a  spherical  object.  The  main  purpose  of  this  initial  data  set 
was  for  a  reasonableness  check.  That  is,  do  the  separate  components  of  the  received  secondary  field 
have  the  right  sign,  reasonable  amplitude,  range,  etc.  Eigure  5.1.1  shows  the  path  of  a  sphere  as  it 
is  moved  along  a  line  below  the  sensor  head  and  is  offset  in  the  positive  y  direction.  Eigure  5.1.2 
shows  the  result  if  of  extracting  the  quadrature  data  recorded  at  450Hz  and  plot  it  as  the  sphere  is 
moved  along  this  path.  Reasonableness  checks  on  this  data  include 

•  Hy  is  always  negative 

•  Hx  starts  negative  then  switches  to  positive 

•  H,  is  always  positive 

•  \Hy\>  \Hx\  due  to  the  proximity  of  the  sphere  to  the  x  —  y  plane 

•  results  are  symmetrical  around  the  yz  plane 
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Figure  5.1.1:  Initial  reasonableness  check  of  the  GEM-  3^+.  Object  is  moved  along  a  path  below  and  offset 
to  the  GEM-3  ^+. 


Figures  5.1.3  and  5.1.4  show  shots  from  along  a  line  erossing  the  eenter  of  the  3  GEM-3^''' 
Rx  coils,  on  a  45°  angle  relative  to  the  Hx  and  Hy  coils.  Spectra  are  shown  for  three  positions: 
the  one  (10  cm)  before  crossing  the  origin  (”pre”),  over  the  center,  and  after  (’’post”)  crossing  the 
origin.  Symmetry  &  appropriate  pattern  (magnitude)  similar  to  those  noted  above  are  evident. 

Figures  5.1.5  and  5.1.6  show  the  components  from  Rx  coils  1  &  2  (horizontal:  Hx  and  Hy). 
Values  are  pretty  small,  not  far  above  the  noise/background  variation  floor  (the  sphere  was  quite 
small).  The  responses  pass  both  quantitative  and  quantitative  reasonableness  checks  based  on  the 
amplitudes  and  signs  of  the  3  components. 


Benjamin  Barrowes,  Kevin  O’Neill 


-105- 


engineering 

AT  DARTMOUTH 


MM- 15  37,  GEM-3  ^+-  Final  Report 


5  RESULTS  AND  DISCUSSION  -  GEM-3^+ 
5.1  Initial  Data  and  Reasonableness  Checks 


Q  at  450  Hz 

Figure  5.1.2:  Initial  reasonableness  check  of  the  GEM-3^+.  Quadrature  response  at  450Hz. 
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Figure  5.1.3:  GEM-3^+  Z  component  of  received  field  with  sphere  directly  over  the  center  of  the  instru¬ 
ment. 
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Z  component 


Figure  5.1.4:  GEM-3 Z  component  of  received  field  with  sphere  to  the  sides  of  the  instrument. 


Components  1  &  2:  Inphase 


Figure  5.1.5:  GEM-3^+  X  and  Y  component  of  received  field  af  3  positions.  Inphase. 
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Components  1  &  2:  Quadrature 


Figure  5.1.6:  GEM-3^+  X  and  Y  component  of  received  field  at  3  positions.  Quadrature. 
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We  also  acquired  data  over  Sin  spheres  of  brass,  aluminum,  and  steel  in  a  9x9  grid.  Figures  5.1.7 
and  5.1.8  show  data  from  different  transects  of  the  received  field  for  the  and  Hx  components 
respectively  for  the  steel  sphere.  These  data  traces  should  be  on  top  of  each  other  if  the  sphere  was 
exactly  aligned.  Since  the  data  are  slightly  off  but  in  a  consistent  direction,  we  conclude  that  we 
are  in  fact  seeing  the  results  of  a  slight  misalignment  from  center. 


Figure  5.1.7:  component  of  data  from  different  transects  of  the  received  field  over  a  9x9  grid. 

Figures  5.1.9  and  5.1.10  show  the  and  the  Hx,  Hy  components  of  the  received  magnetic 
field  respectively  over  a  steel  sphere  located  southeast  of  the  sensor.  As  expected  the  transverse 
components  of  the  received  magnetic  field  are  negatives  of  each  other.  Also,  all  three  components 
of  the  received  magnetic  field  are  available  across  the  entire  band  from  30  Hz  100  kHz. 

Figure  5.1.11  shows  the  standard  deviation  of  the  measured  H  from  the  combining 

10  seconds  of  lOHz  data.  The  values  are  seen  to  be  quite  stable  over  the  short  time  periods  of 
collecting  one  data  show.  The  GEM-3  does  exhibit  some  slight  drift  during  the  warm  up  phase 
(Ril5  minutes)  especially,  similar  to  other  GEM’s. 
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Figure  5.1.8:  component  of  data  from  different  transects  of  the  received  field  over  a  9x9  grid. 


Figure  5.1.9:  component  of  the  received  magnetic  field  over  a  steel  sphere  located  fhe  soufheasf  of  fhe 

sensor. 
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5  RESULTS  AND  DISCUSSION  -  GEM-3^+ 
5.1  Initial  Data  and  Reasonableness  Checks 


Figure  5.1.10:  Hx,  Hy  components  of  the  received  magnetic  field  over  a  steel  sphere  located  the  southeast  of 
the  sensor. 


Figure  5.1.11:  Standard  deviation  of  the  raw  lOHz  data  of  the  3  field  components  of  the  received  magnetic 
field. 
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5  RESULTS  AND  DISCUSSION  -  GEM-3^+ 
5.2  Blind  Test  Suite  #1  Training  Data 


Figure  5.2.1:  GEM-3^+  on  the  lab  grid  for  blind  test  suite  #1. 


5.2  Blind  Test  Suite  #1  Training  Data 

Blind  tests  were  performed  for  the  GEM-3  in  that  laboratory  data  was  aequired  at  CRREL  and 
handed  to  our  eolleague  Eridon  Shubitidze  at  Dartmouth  College.  There  are  two  stages  to  blind 
testing  with  the  NSMS  or  SEA  models.  The  first  stage  involves  aequiring  eharaeterization  data  of 
a  target  in  order  to  set  up  a  library  of  responses  from  possible  targets  in  the  test.  The  seeond  stage 
is  inverting  for  position,  orientation,  and  whieh  target  is  in  eaeh  test.  These  were  single  target  tests. 

The  eharaeterization  phase  for  this  round  of  blind  tests  was  eomplete  in  April,  2008.  Data 
was  taken  on  a  3x3  (20em  spaeing,  blue  dots  in  Eig.  5.2.2  for  baekground  data)  or  a  9x9  (5om 
spaeing)  grid  at  multiple  elevations  and  orientations  of  eaeh  library  target.  Eigure  5.2.1  shows 
the  GEM-3^+  on  the  lab  measurement  grid,  and  Eig.  5.2.2  shows  the  grid  numbering  system 
used  for  this  data.  Eigure  5.2.3  shows  a  60mm  mortar  in  the  3  test  eonfigurations  used  to  aequire 
eharaeterization  data.  Generally,  the  three  9x9  measurement  grids  were  taken  at: 

•  horizontal,  nose  north  (or  -i-y),  abbrev  HNN. 

•  nose  north  45degrees  down  from  horizontal,  abbrev  NND 

•  nose  north  45degrees  up  from  horizontal,  abbrev  NNU 

Depth  measurements  are  to  the  nearest  point  with  the  eenter  of  the  target  eentered  at  (0,0),  and 
baekground  grids  of  3x3  were  aequired  in  between,  before,  and  after  eaeh  9x9  measurement  grid. 
The  suite  of  targets  we  aequired  eharaeterization  data  for  this  test  are  pietured  in  Eig.  5.2.4. 
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5.2  Blind  Test  Suite  #1  Training  Data 


North 


Figure  5.2.2:  Grid  layout  used  for  blind  test  suite  #1. 


Figure  5.2.3:  60mm  mortar  in  measurement  setup  for  blind  test  suite  #1. 
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5  RESULTS  AND  DISCUSSION  -  GEM-3^+ 
5.3  Blind  Test  Suite  #1  Results 


Figure  5.2.4:  Possible  targets  in  library  for  test  suite  #1. 


5.3  Blind  Test  Suite  #1  Results 

A  library  was  constructed  based  upon  the  lab  characterization  data  described  above.  The  the  HAP 
method  and  the  NSMS  model  (see  Sec.  4.6)  were  applied  to  subsequent  blind  test  data  (also  on 
a  9x9  grid).  The  total  normalized  magnetic  charge,  Q  for  each  blind  test  target  was  compared 
to  the  library  Q  and  ranked  according  to  best  match.  Figures  5.3. 1-5. 3. 8  show  Q  for  the  blind 
test  data  compared  to  that  of  the  best  match  from  the  library.  In  summary.  Fig.  5.3.9  shows  the 
geometry  and  tabulates  the  results  from  this  blind  test.  All  8  targets  were  identified  successfully 
with  positions  and  orientations  very  close  to  the  actual.  This  blind  test  is  encouraging  in  that  it 
suggests  the  NSMS  model,  combined  with  the  HAP  method  for  fast  target  pinpointing,  is  very 
effective  at  correctly  discriminating  between  targets  with  vector  data  such  as  that  provided  by  the 
GEM-3^+. 
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Figure  5.3.1:  Total  NSMS  of  blind  test  1. 
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Figure  5.3.2:  Total  NSMS  of  blind  test  2. 
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5.3  Blind  Test  Suite  #1  Results 


Figure  5.3.3:  Total  NSMS  of  blind  test  3. 


Figure  5.3.4:  Total  NSMS  of  blind  test  4. 
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5.3  Blind  Test  Suite  #1  Results 


Frequency  [Hz] 


Figure  5.3.5:  Total  NSMS  of  blind  test  5. 


Figure  5.3.6:  Total  NSMS  of  blind  test  6. 
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Figure  5.3.7:  Total  NSMS  of  blind  test  7. 
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Figure  5.3.8:  Total  NSMS  of  blind  test  8. 
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Blind 
test  # 

Estimated  (Actual) 

6(0) 

IP) 

X,  [cm] 

y,  [cm] 

z,  [cm] 

1 

93.60C  45) 

90.0  (  90) 

3.0  (  0) 

0.01(  1.5) 

-29.68(-28) 

2 

Se.l/C  90) 

20.9  (135) 

-5.5  (  -5) 

5.03(  5) 

-34.46(-31) 

3 

105. 36(  0) 

76.4  (arb) 

-15.6  (-20) 

-15.8  (-15) 

-43. 46 (-42) 

4 

91.08(arb) 

2.25(arb) 

18.2  (  20) 

-0.5  (  0) 

-25.19  (-22) 

5 

81.05  (135) 

10.11(135) 

20.2  (  15) 

-0.7  (  5) 

-26.22  (-26) 

6 

86.  34(180) 

67. 05 (arb) 

5.3  (  10) 

-12.64(-20) 

-20.25 (-18) 

7 

65.56(  90) 

16.63(  45) 

12.9  (  18) 

4.13(  12) 

-25. 81  (-20) 

S 

66.06(  90) 

14.01(  90) 

3.35  (  5) 

15.45(  10) 

-43.10(-40) 

Figure  5.3.9:  Geometry  diagram  and  positions  for  blind  tests  derived  from  the  HAP  method  compared  to 
truth. 
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Figure  5.4.1:  Contents  of  test  plot  #2  blind  test  suite  #2. 


5.4  Blind  Test  Suite  #2  -  CRREL  Test  Plots  with  Beacon  Positioning 

On  October  23,  2008,  CRREL  personnel  acquired  data  over  15  target  clusters  in  the  CRREL  test 
plots.  Each  target  cluster  consisted  of  either  one  or  2  emplaced  targets.  The  test  plot  emplacements 
consisted  of  5  target  clusters  in  3  separate  20’x20’x8’  concrete  plots.  The  items  in  each  cluster 
are  shown  in  Eigs.  5.4. 1-5. 4. 3.  In  these  figures,  a  red  box  around  an  alternate  target  indicates  a 
target  that  used  to  be  in  that  position,  but  was  thereafter  replaced  with  the  target  listed.  Target 
designations  such  as  LG-22  and  so  forth  are  CRREL  designations  and  refer  to  different  pieces  of 
scrap  with  specifications  available  upon  request.  Test  plot  #1  was  not  used  for  this  suite  of  tests. 

The  GEM-3^+  including  the  functioning  beacon  positioning  system  was  used  to  acquire 
data  over  all  15  target  clusters.  Appendix  C.3  describe  the  data  processing  of  this  data  while 
in  Sec.  5. 4. a,  the  results  of  a  multi-object  NSMS  inversion  [98,  99,  134]  is  presented.  References 
to  “lithos”  and  data  paths  in  general  are  CRREL/Dartmouth  specific.  All  relevant  files  are  available 
upon  request. 
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Figure  5.4.2:  Contents  of  test  plot  #3  blind  test  suite  #2. 
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Figure  5.4.3:  Contents  of  test  plot  #4  blind  test  suite  #2. 
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5.4  Blind  Test  Suite  #2 


5.4.a  Blind  Test  Suite  #2  Results 

The  results  for  the  multiobjeet  NSMS  inversion  are  presented  on  the  next  two  pages. 

The  inverted  results  for  this  multi-object  NSMS  method  are  mixed  but  encouraging.  The  it¬ 
erative  multi-object  approach  here  is  a  slow,  somewhat  “brute  force”ish  method.  Also  remember 
that  this  data  has  background  subtraction  issues  as  outlined  in  the  last  sections.  Even  with  these 
difficulties,  the  technique  here  was  able  to  correctly  identify  individual  UXO  in  about  half  of  the 
cases  presented  (the  green  highlighted  boxes).  For  many  other  cases,  the  method  correctly  distin¬ 
guished  between  UXO  and  clutter  (yellow  boxes),  but  only  rarely  produced  a  full  misclassification 
(red  boxes). 
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GEM  3D  plus  processing  results  for  multi  objects  : 
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6  Conclusions  and  Implications 

SERDP  project  MM-1537  entitled  “Handheld  Frequency  Domain  Vector  EMI  Sensing  for  UXO 
Discrimination”  is  complete.  The  SERDP  project  MM1537  has  resulted  in  a  unique  frequency 
domain  instrument  that  adds  transverse  receivers  to  the  standard  GEM-3  platform  while  also  in¬ 
corporating  several  advances  in  hardware  and  software.  The  most  notable  other  hardware  advance 
aside  from  the  transverse  receiver  coils  (the  3^  part  of  GEM-3^'''),  a  positioning  system  based 
on  the  primary  field  of  the  instrument  itself  has  been  designed  and  installed  which  yields  subcen¬ 
timeter  accuracy  out  about  2.5  meters  from  the  beam  containing  two  triaxial  receivers. 

Alongside  the  hardware  advances  made  during  this  project,  the  accompanying  software  and 
models  are  also  of  note.  The  normalized  surface  magnetic  charge  (NSMS)  model  was  advanced 
and  adapted  for  the  GEM-3^’'^.  As  well  the  Standardized  Excitations  Approach  (SEA)  was  sim¬ 
ilarly  adapted  and  applied  to  ]gemd  data.  These  physics  based  rigorous  models  boast  arbitrary 
fidelity  and  can  take  into  account  any  primary  field  and  predict  the  response  from  any  physical  tar¬ 
get.  A  method  to  extract  an  absolute  scale  which  translates  GEM-3  ppm  into  physical  quantities 
(Amps/meter)  was  also  derived  and  published. 

The  GEM-3^”'"  is  now  a  fully  functioning  EMI  instrument  capable  of  vector  sensing  of  mag¬ 
netic  anomalies  while  being  well  located  within  a  limited  range.  The  transverse  receivers  add 
critical  data  diversity  and  information  to  rigorous  models  and  help  in  inversion  routines.  For  pre¬ 
cise,  queued  interrogation  of  anomalies,  the  GEM-3^”''  provides  diverse,  accurate,  frequency  do¬ 
main  data  of  the  secondary  EMI  field  suitable  for  inversion  and  discrimination  with  high  fidelity, 
rigorous  models. 
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C  OTHER  SUPPORTING  MATERIAES 


A  Appendices  -  Supporting  Data 
B  List  of  Scientific/Technical  Publications 

Publications  produced  in  whole  or  part  with  funds  from  this  projeet  during  FY06  are  [5,  98,  105, 
no.  111,  121,  126,  134-139], 

C  other  Supporting  Materials 

C.l  Data  Acquisition  Sequence  from  the  GEM-3 


Benjamin  Barrowes,  Kevin  O’Neill 


-Appendix- 

-127- 


ENGINEERING 
AT  DARTMOUTH 


MM- 15  37,  GEM-3  ^+-  Final  Report 


C  OTHER  SUPPORTING  MATERIAES 
C.l  Data  Acquisition  Sequence  from  the  GEM-3^"'" 


Running  Gem-3D  with  GPS  unit. 


_  More  than  a 
meter  awav  " 


You  must  follow  these  instruction  to  the  letter  or  you  will  not  have  good  data. 

Turn  on: 

1. )  Gem-3D  -  Battery  or  Power  Supply  -  11.6  volts  DC 

2. )  GPS  -  Battery  or  Power  Supply  -  12  volts  DC 

3. )  Computer 

4. )  Start:  WinGem  v4 

\VmGEMv4 

Press  the  Transmit  button  -  Then  the  Null  Em  Button  -  zero  it. 


3b' start!  ,  zr3  E  E3  -fi»  |  E  WinCEM  -  G30-1 


Easy  Button  Mar 


|Ln:?  FIDiO  100:00;lI.t55 
12:52  AM 


5.)  Start:  GEMPOS 


GEMPOS 
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C  OTHER  SUPPORTING  MATERIAES 
C.l  Data  Acquisition  Sequence  from  the  GEM-  3^+ 


6.)  Press  the  Transmit  button  -  Then  the  Null  Em  Button  -  zero  it. 


7.)  Then  minimize  reduce  screen  by  hitting  the  minus  sign. 
Do  not  shut  off,  to  silence  beeping. 


.=jn|xj 


8.)  Go  back  to  WiuGem  -  you  are  ready  to  test  -  Press  the  Store  Button 
You  can  use  automatic  file  names  or  create  your  own. 
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0  bbb. survey. GEM. IDOl.c 


File  name:  |t_^m3d.sufV^.GEM.csv 

Save  as  type:  |  Survey  Files  P.Survey.GEM.IDOl  .csv) 


“3 

"3 


Ready 

1 - 

It-  :-KBPS:ERltSn:  123  lExt.Pow.  ]EM:  000  |Lni?  FID:  0  |00;04;I5.164 
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9. ) 

10. ) 

11.) 

12.) 


Press  Store  again  when  test  is  finished. 

Then  minimize  screen  by  hitting  the  minus  sign. 

Do  not  shut  off. 

You  must  run  Merge  after  each  test  or  data  will  be  lost. 
Start:  Merge 

merge 


^xj' 
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C  OTHER  SUPPORTING  MATERIAES 
C.l  Data  Acquisition  Sequence  from  the  GEM-  3^+ 


13.)  Press  Select  Survey  -  fine  the  test  that  you  have  just  finished  and  open  it 


14.)  Then  press  Merge  LPS-3d  Survey  button.  When  Merge  is  finished  it  will  close 
automatically. 


“^Untitled  -  MERGER 
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Serrsor  Offsets  From  GPS  Center  Point 

P . - 3 


Select  Survey 


j  r 


{  010CT08  124F55^:.r,.-:- 


Survey  Type 

LPS-3D 

Standard 

Merge  IPS -30  Survey 

Ready 


15.) 


Check  data  with 


Windows 

Explorer 


and  Microsoft 
Office  Excell 
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16.)  data  should  look  like  this. 


Not  this. 


•  Ble  Edit  !^w  Insert  Fsrmet  loots  Qata  \i0ndow  Help 

i  j  tji  d  ^  .jiy  *  -j  a-  y  ■ 


Type  a  ouesbon  for 


;  Draw  *  ]  AutoShapes  ♦N  Nl  = 


I  ,^tem=Loc8lPo8 


^R^.GEAUOQ  1  Jka.csv 


Z] 

1  A  1  B 

m 

rsn 

rnn 

m 

1  1 

.System=Loca!Pos 

2 

.45  Total  Samples  In  File 

“3^ 

.Info  Date;  2008 

Time.  19J 

DataLabel  ID01 

jH 

Xl 

;NoteT 

.Note2: 

1 

1 

i 

L  1 

[ 

6 

7 

X  Y 

z 

alphV 

BETA 

Q 

LPSTtm 

LPSmsO 

1  •  1 

PowerLnII  30Hz 

1 2.38E+09  3"90E+02 

8 

* 

.  1 

[  •  1 

I2.38E+09 

3.89E+02 

9 

•  • 

• 

• 

• 

2.35E+09 

3  91E+02 

10 

2  36E+09 

3.89E+02 

11 

2.34E+09 

3  90E+02 

12 

2.36E+09 

3.90E+02 

13 

• 

* 

• 

2.38E+09 

3.89E+02 

14 

2.38E+09 

3  90E+02 

15 

*  1 

2.3SE+09 

3  90E+02 

16 

• 

* 

^  1 

• 

2.34E+09 

3  90E+02 

17 

1  • 

2  37E+09 

3  90E+02 

18 

] - - 

] - - 

tz? 

2  37E+09 

3  90E+02 

I  P  TT 


90Hz  i1_210Hz 
3  89E+02 
3.89E+02 
3  89E+02 
3.89E+02 
3.89E+02 
3.89E+02 
3.89E+02 
3.89E+02 
3  89E+02 
3.89E+02 
3  89E+02 
3.89E+02 


Q1_210Hzi1_450Hz 
limilltlllll  3*87E+02 
umKHIt#  3.87E+02 
unttmuil  3  87E+02 
m#mil  3.87E+02 
Itlll/UltUlfll  3.87E+02 
limtllUlUtlt  3  87E-4’02 
tHUIIIHIIHIf  3  87E-r02 
3.87E+02 
3  87E+02 
3.87E+02 
3  87E+02 
3.e7E+02 


17.)  You  are  ready  to  start  your  next  test. 
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C.2  Beacon  Positioning  System  Construction  Details 
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C.2  Beacon  Positioning  System  Const];yg^i^g(fj^^s 

Progress  Report  (October  2006) 


Progress  Report 
on 

Handheld  Frequency-Domain  Vector  EMI  Sensor  for  UXO  Discrimination 
Contract  Number:  W913E5-06-C-001 1 


April  2007 


from 

Geophex,  Ltd. 
605  Mercury  Street 
Raleigh,  NC  27603 


1 .  Progress  Summary 

This  report  covers  our  progress  on  the  Beacon  navigator.  CRREL  issued  a  contract 
modification  early  2007  so  that  this  task  is  completed  in  FY07.  As  of  this  reporting,  we  have 
completed  the  physical  construction  of  the  coils  and  housing.  Early  this  month,  we  started 
electronic  wiring  and  assembly. 

2.  Navigator  Geometry  and  Construction 

The  Beacon  Navigator  is  to  be  used  for  accurately  locating  the  GEM-3D  active  EMI  sensor 
Geophex  built  for  CRREL.  Figure  1  shows  how  it  works:  the  GEM  transmitter  acts  as  an  active 
beacon  that  can  be  located  by  a  pair  of  receiving  coils  in  the  vicinity,  each  of  which  consists  of 
three  orthogonal  coils.  This  method  accurately  produces  the  sensor's  position,  as  well  as  the 
sensor's  orientation  as  defined  by  two  tilt  angles.  Tests  have  shown  that  the  method  produces 
accuracy  of  better  than  a  centimeter  when  the  receiving  coils  are  at  a  range  of  several  meters. 


Two  3-Axis 
Receiving  Coiis 

(XO,  YO,  ZO) 

(X2,  Y2,  Z2)  ^ 

Sensor 
Transmitter 

Figure  1.  A  pair  of  3D  coils  measure  the  magnetic  field  generated  by  the  EMI  sensor.  The 
six  coil  outputs  are  sufficient  to  determine  both  the  sensor  location  and  tilt  angles. 

We  built  in  2006  a  prototype  Beacon  Navigator  housed  in  a  10cm  x  10cm  x  L5m  PVC  tube 
that  contains  a  3D  coil  at  each  end.  The  output  of  the  prototype  navigator  is  routed  by  hardwire 
(~10m  long)  to  the  GEM-3  sensor.  The  Navigator  computes  the  sensor  location  using  the 
voltage  output  from  the  coils  and  sends  the  results  to  the  GEM-3  that  merges  the  location  data 
with  the  sensor  data.  Under  this  CRREL  project,  we  will  refine  the  existing  navigator  by 
improving  and  simplilying  its  performance,  and  producing  an  advanced  navigator.  We  will  then 


{X1,Y1,Z1) 
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Progress  Report  (October  2006) 

integrate  the  Navigator  with  the  GEM-3D,  using  the  Bluetooth  wireless  communication  so  that 
it  will  be  completely  detached  from  the  sensor  operator. 

Calibration:  Since  the  relative  geometry  of  the  two  3D  navigator  coils  is  fixed,  only  a  few 
ground-truth  points  are  sufficient  for  a  full  calibration.  One-time  calibration  is  sufficient  for  its 
lifetime  usage. 

3.  Progress  Status 

3.1  Navigator  Sensing  Coils 

Table  1  shows  the  preconstmction  design  figures  of  the  three  nested  coils.  The  z-axis  coil  (the 
middle  of  the  nest),  for  instance,  has  a  mean  winding  cross-section  of  6.8  cm  x  6.8  cm  and  has 
2,000  turns  of  #34-guage  wire.  The  turns  for  other  nested  coils  are  adjusted  to  have  a  similar 
area-turns  product  to  equalize  their  voltage  outputs.  Table  1  also  shows  predicted  inductance 
(in  Henry)  and  resistance  (in  ohm)  for  each  coil. 


As  designed 


-  Z  Middle  Section 

2000  turns 

0.873  AWG#34 

9.25 

area  X  turn  (m*m) 

0.068  a  (m) 

0.334  inductance  (H) 

64 

x-sec  (mm**2) 

0.068  b  (m) 

474.9  resistance  (ohm) 

0.703 

time  constant  (msec) 

-  X  Outermost 

1784  turns 

0.873  AWG#34 

9.25 

area  x  turn  (m*m) 

0.072  a  (m) 

0.290  inductance  (H) 

57.088 

x-sec  (mm**2) 

0.072  b  (m) 

448.5  resistance  (ohm) 

0.647 

time  constant  (msec) 

-  Y  Innermost 

2260  turns 

0.873  AWG#34 

9.26 

area  x  turn  (m*m) 

0.064  a  (m) 

0.387  inductance  (H) 

0.766 

time  constant  (msec) 

0.064  b  (m) 

505.1  resistance  (ohm) 

Table  1.  Design  figures  of  the  three  nested  coils 

Figure  2  shows  actual  coils  constructed  based  on  the  design  figures.  Each  of  the  two  sensing 
coils  is  wound  on  a  light-density,  3 -inch  plastic  cube  that  contains  a  pair  of  recessed  grooves 
for  coil  on  each  axis.  The  coils  are  to  be  mounted  at  the  ends  of  a  3”  x  3”  x  5  ft  square  PVC 
tube,  as  shown  in  Figure  3. 

Table  2  shows  the  measured  inductance,  resistance,  and  resonance  frequency  of  actual  coils 
that  were  built  based  on  the  design  figures.  The  predicted  and  actual  are  close  within 
manufacturing  accuracy. 
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_ _ Progress  Report  (October  2006) 


Figure  2.  Two  3-D  nested  sensing  coils  with  preamps  attached 


Figure  3.  Square  PVC  tube  housing 


Actual 

Coil  1 

Coll  2 

Rx-Z 

half  half  total 

0.116  0.117  0.323 

240  241  481 

half  half  total 

0.117  0.117  0.334 

235  235  470 

Ind  (H) 
res  (Ohm) 

14.7 

14.6 

res  freq  (kHz) 

Rx-X 

0.102  0.101  0.288 

230  229  460 

0.101  0.101  0.287 

226  226  452 

Ind  (H) 
res  (Ohm) 

16.3 

15.7 

res  freq  (kHz) 

Rx-Y 

137  137  0.393 

263  263  537 

0.138  0.138  0.394 

253  253  507 

Ind  (H) 
res  (Ohm) 

13.8 

13.1 

res  freq  (kHz) 

Table  2.  Measured  inductance,  resistance,  and  resonance  frequency  of  the  actual  coils 
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Progress  Report  (October  2006) 


3.2  Sensing  Electronics 

The  final  electronic  design  for  the  navigator  is  still  in  progress.  Two  primary  design  issues  are: 

•  Wireless  communication  with  the  GEM-3D  and 

•  Automatic  gain  control  scheme  for  each  preamp  to  prevent  input  saturation  when  the 
GEM-3  is  close  to  the  navigator. 

The  second  issue  is  a  newly-discovered:  when  the  GEM-3D  is  at  less  than,  say,  1  m  from  the 
navigator,  the  front-end  voltage  becomes  high  to  saturate  the  preamp.  We  are  currently 
working  on  an  automatic  gain-control  scheme  to  deal  with  this  issue. 

Figure  4  shows  the  main  DSP  section  that  will  process  the  6-channel  time-series  data. 


Figure  4.  DSP  electronics  for  the  navigator 

4.  Plan  for  the  Next  Reporting  Period 

We  will  continue  the  electronic  assembly  of  the  navigator.  We  expect  to  fire  it  up  within  the 
next  month  and  to  begin  the  test  and  calibration  phase. 

5.  Estimated  Percentage  of  Contract  Completion  for  Year  2  Task:  40% 
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C.3  GEM-3  Data  Processing  -  Static  and  Dynamic 
C.3.a  Introduction 

This  short  document  explains  the  internal  strueture  of  the  program  for  data  read-in  and  analysis  of 
GEM-3^+  files. 

Two  oases  are  eonsidered: 

1.  The  “Beaeon  sliding  on  a  grid”  data  taken  in  the  lab  by  Ben  Barrowes  on  10/7/08.  This  set  is 
stored  on  lithos  in  direetory 

DATA/ gemSdplus/ 3D_1 0-7-2 0 0 8_Ben/ 

and  is  henoeforth  oalled  Ben. 

2.  The  “Beaeon  free  form”  data  taken  over  the  test  plots  by  Kevin  O’Neill  on  10/23/08.  This 
set  is  stored  in 

DATA/GEM-3D+_OCT-23-08_TESTPLOTS_EXPORT/ 

and  is  oalled  Kevin  from  now  on. 

From  now  on  when  I  refer  to  a  “file”  I  really  mean  a  oolleotion  of  three  .  csv  files:  IDOl, 
ID02,  and  ID03,  oorresponding  to  the  GEM-3^'''’s  three  reoeivers. 

C.3.b  Reading  the  data 

1.  All  filenames  start  with  an  identifying  number.  The  ranges  (of  usable  files)  are  [21  :  49] 
(29  files)  for  Ben  and  [1,3, 4, 8, 9, 12  :  18,21,22]  (14  files)  for  Kevin.  The  variable  Ncase 
ranges  from  1  to  29  (14). 

2.  Some  files  have  a  strange  feature:  IDOl  has  one  more  line  than  either  ID02  or  ID03.  This 
happens  both  in  Ben’s  and  in  Kevin’s  data. 

Kevin’s  writeup  mentions  a  “slight  timing  asynohronization”  in  the  data  he  studied  ( ( 2  - 1 ) , 
or  3_*  in  our  present  notation),  whioh  happens  to  be  one  of  those  eases.  He  suggests  two 
solutions:  interpolating  or  just  ignoring  the  small  differenees. 

I  have  taken  the  seeond  route  and  simply  deleted  the  first  row  when  the  set  is  problematie. 
At  some  point  I  tried  interpolating  (sinee  eaeh  datafile  ineludes  the  time  at  whieh  eaeh  point 
was  measured)  but  did  not  see  mueh  of  a  differenee.  I  ean  redo  this  if  neeessary. 

3.  The  previous  version  of  the  data-reading  proeedure  made  some  use  of  UNIX  system  ealls 
to  determine  the  filenames  and  read  in  the  values.  I  fixed  that  feature,  whieh  made  it  both 
system-dependent  and  slow,  and  replaeed  it  with  Matlab-only  funetions. 

As  it  stands,  the  proeedure  ean  read  .  cs  v-like  files  that  need  not  be  all  numerie.  The  only  ad 
hoes  here  respond  to  Kevin’s  having  deleted  most  of  the  headers  and  the  “EPSTime”  eolumn 
that  presumably  was  preventing  him  from  using  csvread. 

4.  All  three  files  eorresponding  to  a  given  “file”  are  read  at  onee  and  stored  in  a  3D  array. 
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5.  Each  receiver  provides  its  own  position/tilt  information.  I  have  kept  this  as  both  a  2D  matrix 
(e.g.,  Xpos,  with  leading  eapital,  is  Npt  x  3,  where  Npt  is  the  number  of  measurements) 
and  as  an  averaged  array  (xpos,  all  lowerease). 

We  do  note  that  question  arises  as  to  the  aeeuraey  of  the  measurement  loeations  given  by  the 
beaeon  system.  A  look  at  Xpo  s  and  the  other  two  sueh  variables  reveals  two  eolumns  that  are 
identieal  and  distinet  from  the  third.  The  differenees  are  most  notable  for  the  x-eoordinate, 
where  they  reaeh  3  em  and  are  not  random  but  periodie  and  roughly  ;r/2  out  of  phase  with 
the  aetual  x-values.  (See  Fig.  C.3. 2  for  a  plot.) 

6.  We  also  delete  rogue  points  in  whieh  the  loeation  information  was  patently  wrong,  for  ex¬ 
ample  featuring  large  diseontinuities.  In  Ben’s  ease,  it  is  straightforward  to  find  those  points, 
beeause  they  had  a  yaw  angle  of  exaetly  zero.  In  Kevin’s  ease  we  had  to  find  them  by  hand. 


C.3.C  Background  subtraction 

C.3.c.(l)  Initial  approach  We  start  by  noting  that  we  negleeted  the  angular  information  in  the 
baekground  subtraetion.  In  all  eases  some  effort  was  put  into  holding  the  sensor  as  level  and 
well-aligned  with  the  beaeon  eoordinate  system  as  possible. 

Ben  took  a  whole  baekground  file,  Ncase  =  2,  whieh  makes  the  subtraetion  straightforward: 
just  read  in  that  file,  average  over  all  points,  and  subtraet  the  average  from  eaeh  data  file.  We  do 
not  have  to  be  too  observant  of  the  height  in  this  ease,  as  there  is  no  soil. 

Kevin’s  files  are  a  little  different,  sinee  he  takes  the  baekground  values  from  the  data  them¬ 
selves.  (Also,  his  “ealibration”  file  is  over  a  shotput.)  In  his  words. 

The  program  next  asks  you  to  enter  some  data  point  number  around  whieh  it  will  aver¬ 
age  (±2  points)  to  obtain  a  baekground  to  subtraet  from  all  data  treated  subsequently. 

It’s  often  hard  to  eome  up  with  a  good  ehoiee,  whieh  has  implieations  for  deeper  pro- 
eessing.  Here  we’re  just  doing  some  exploratory  visualization  so  it’s  not  so  important. 

Now,  however,  is  the  time  to  perform  the  “deeper  proeessing”  he  mentions.  We  eannot  use  Ben’s 
baekground  file,  sinee  that  does  not  inelude  soil  response. 

On  the  other  hand,  Kevin  ehooses  his  baekground  point  thus: 

In  any  ease,  in  the  favorable  ease  shown  here,  point  #25  oeeurs  after  the  surveying  has 
gotten  under  way  and  is  around  the  beginning  of  the  first  sweep;  and  there  isn’t  mueh 
going  on  in  either  of  the  two  response  measures.  Thus  we  enter  this  and  the  program 
averages  over  points  23-27,  plotting  the  resulting  baekground  profile  so  you  ean  see  if 
anything  waeky  is  going  on,  i^  more  strange  than  usual. 

This  gives  rise  to  two  related  possibilities  for  finding  the  baekground: 

1 .  Deeide  upon  a  fixed  number  of  baekground  points,  sort  the  measurements  in  deereasing  order 
of  distanee  from  the  eenter,  and  use  the  (averaged)  furthest  measurements  as  baekground. 

2.  Carry  out  a  multi-file  survey  and  use  the  very  furthest  measurements  as  baekground,  after 
averaging  them. 


Now,  how  do  we  determine  an  “optimal”  number  of  baekground  points  that  will  not  skew  the 
measurements?  We  must  look  at  the  way  the  (frequeney  by  frequeney)  average  for  the  furthest 
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points  change  as  we  inerease  or  deerease  Nbk;  also,  we  want  to  see  how  it  eompares  to  the  field 
values  themselves. 

It  turns  out  the  baekground  values  so  ehosen  exeeed  the  eombined  baekground  -P  UXO  mea¬ 
surements  in  a  few  eases,  and  moreover  the  response  is  not  flat.  This  may  well  be  due  to  the  height 
dependenee  of  the  soil  response,  whieh  is  nominally  On  the  other  hand,  the  response  of  a 
dipole  is  approximately  J  whieh  lends  eonfidenee  to  our  taking  the  points  presumably  furthest 
from  the  target. 

The  user  ean  eontrol  (dietate)  Nbk,  the  number  of  points  to  take  as  baekground.  The  program 
finds  a  “eenter”  point  defined  by  =  XcX  +  ycf  +  zqz,  where  Xc  and  are  the  averages  of  x  and  y 
and  zo  is  the  lowest  possible  value  of  z,  eenters  the  positions  there,  and  measures  the  radial  distanee 
from  that  origin  to  every  measurement  point.  It  then  eounts  the  Nbk  points  furthest  from  the  origin, 
averages  their  field  values,  and  uses  this  as  a  baekground.  The  funetion  takes  Nbk  =  20  as  default, 
though  that  value  ean  be  ehanged. 


C.3.c.(2)  Two-height  approach  The  method  deseribed  above  has  some  flaws.  First,  it  mashes 
together  “baekground”  data  eorresponding  to  different  heights,  whieh  ean  bias  the  results.  Also, 
the  distanee  from  the  “eenter,”  whieh  we  estimate  as  the  average  of  the  x-  and  y-loeations,  may 
not — and  in  multi-object  cases  will  not — eoineide  with  the  eenter  of  anything,  and  thus  it  may  turn 
out  that  the  baekground  points  are  in  faet  right  on  top  of  a  target. 

We  then  had  to  find  a  baekground-subtraetion  teehnique  that  would  try  to  overeome  these  prob¬ 
lems.  In  the  end  we  implemented  an  ad  hoe  two-height  approaeh. 

Kevin’s  measurements  consist  of  two  zigzags,  one  elose  to  the  ground  making  five  passes  over 
the  measurement  area  and  the  other  retraeing  the  original  in  the  opposite  direetion  and  some  10  em 
above.  The  measurements  were  taken  by  hand,  so  the  separations  between  the  two  zigzags  are 
not  neeessarily  elear-eut,  and  in  some  eases  there  are  “higher”  points  that  are  aetually  eloser  to 
the  ground  than  some  in  the  “lower”  zigzag.  In  our  proeedure,  however,  we  assume  two  distinet, 
eonseeutive  heights. 

We  start  by  finding  by  hand  a  single  point,  halfway  in  the  transition  region  between  the  two 
zigzags,  and  divide  eaeh  data  set  into  the  “before”  and  “after”  intervals.  We  then  find  a  five- 
eonseeutive-point  interval  for  eaeh  of  the  zigzags  that  represents  the  best  eompromise  between  the 
following  desirable  (and  often  eonflieting)  features: 

1.  The  interval  is  far  from  the  eenter  of  the  measurement  trajeetory. 

2.  The  points  lie  at  roughly  the  same  height,  one  representative  of  the  eorresponding  zigzag 

3.  The  five  points  lie  along  a  fairly  smooth  path  so  the  mean  of  the  field  values  is  a  reasonable 
average. 

4.  The  resulting  fields  will  have  these  eharaeteristies: 

•  The  z-field  should  be  positive  in  as  many  points  as  possible. 

•  The  X-  and  y-fields  should  have  roughly  as  many  positive  values  as  negative,  and  the 
zero  erossings  should  be  roughly  at  the  eenter  of  eaeh  sweep. 

We  eoneentrated  on  the  quadrature  parts  of  the  fields,  since  they  are  presumably  less  affeeted 
by  ground  effeets  and  did  not  have  the  large  anomalous  de  shift  that  the  inphase  parts  had. 
Moreover,  the  quadrature  response  is  positive  at  all  frequencies,  so  any  negative  amplitudes 
will  be  due  to  geometry. 
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5.  The  above  rules  should  hold  in  as  large  a  range  of  frequeneies  as  possible. 

It  is  this  method  that  we  used  in  the  end  to  prepare  the  data  for  further  proeessing. 
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Figure  C.3.1:  Sensor  position  in  the  example  data  run.  The  convention  in  this  and  the  following  figures  is: 
the  sensor  “warms  up”  as  it  takes  data:  i.e.,  the  trajectory  is  initially  blue  and  transitions  into  red  during  the 
run.  The  background  points  are  black  dots. 

C.3.d  An  example  of  the  first  approach 

As  an  example  we  inelude  our  analysis  of  Kevin’s  file  Ncase  =  2,  corresponding  to  location  2-1. 

For  this  example  we  use  the  first  approach  to  background  subtraction.  We  choose  Nbk  =  10 
in  this  case;  the  points  happen  to  lie  at  the  very  end  of  the  data  run,  and  appear  as  black  dots  on 
Fig.  C.3.1.  The  black  star  in  that  plot  is  the  center  from  which  the  background  points  are  most 
distant. 

Figures  C.3. 3  and  C.3.4  respectively  show  the  inphase  and  quadrature  parts  of  the  background- 
subtracted  secondary  field  at  /  =  450  Hz.  On  the  left-hand- side  panel  we  have  the  raw  values.  The 
dashed  lines  represent  the  background  value,  computed  by  averaging  the  field  values  component 
by  component  at  the  Nbk  points,  which  appear  in  the  figure  as  dots.  On  the  right-hand- side  panels 
we  have  the  values  obtained  after  subtracting  the  background. 

The  quadrature  parts  are  further  depicted  as  “profiles”  on  Fig.  C.3. 5. 
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Figure  C.3. 2:  Discrepancy  in  the  a:  position  between  the  first  and  third  data  files  of  a  set. 
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Inphase  fields  for  f  =  450;  z+1 3000  Inphase  normalized  fields  for  f  =  450 


Figure  C.3. 3:  Background- subtracted  field  at  /  =  450  Hz,  inphase  part 
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Quadrature  fields  for  f  =  450 


Figure  C.3.4:  Background-subtracted 


Quadrature  normalized  fields  for  f  =  450 


field  at  /  =  450  Hz,  quadrature  part 
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Figure  C.3. 5:  Background-subtracted  3D  fields  at  /  =  450  Hz. 
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Figure  C.3. 6:  Quiver  plot  showing  the  projection  on  the  x-y  plane  of  the  quadrature  part  of  the  secondary 
field  at  /  =  450  Hz.  Joining  the  tails  of  the  arrows  can  give  a  approximate  estimate  of  the  locations  of  the 
targets. 
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Figure  C.3.7:  Sensor  position  in  the  second  example  data  run. 


C.3.e  Example  of  the  second  approach 

To  illustrate  the  second,  final  background-subtraction  method  we  take  another  example,  location 


2-4. 


Figure.  C.3.7  again  plots  the  locations,  with  labels  and  hollow  markers  identifying  every  10th 
point.  The  background  points  can  be  seen  as  solid  black  dots,  to  the  left  and  right  of  the  figure. 

Figures  C.3.8  and  C.3.9  respectively  show  the  inphase  and  quadrature  parts  of  the  background- 
subtracted  secondary  field  at  3930  Hz.  The  description  is  as  before:  the  left-hand-side  panel  shows 
raw  values,  the  dashed  lines  represent  the  (two-level)  background  level,  the  dots  depict  the  points 
used  to  find  the  background.  On  the  right-hand- side  panels  we  have  “normalized”  fields. 

The  field  profiles  of  Fig.  C.3.10  look  two-peaked,  and  the  impression  that  this  plot  may  have 
two  objects  is  strengthened  by  looking  at  the  quiver  plot  in  Fig.  C.3.1 1. 
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Figure  C.3.10:  Background-subtracted  3D  fields  at  /  =  3930  Hz. 
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Figure  C.3.11:  Quiver  plot  showing  the  magnitudes  of  the  x-  and  ^-components  of  the  quadrature  part  of 
the  secondary  field  at  /  =  3930  Hz.  Joining  the  tails  of  the  arrows  can  give  a  approximate  estimate  of  the 
locations  of  the  targets. 
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